Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,966,573 T→G D621A (GAT→GCT)  ptsP ← PEP‑protein phosphotransferase enzyme I; GAF domain containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,966,5730TG98.5% 265.6 / ‑4.2 68D621A (GAT→GCT) ptsPPEP‑protein phosphotransferase enzyme I; GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (37/30);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

CACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTC  >  NC_000913/2966427‑2966718
                                                                                                                                                  |                                                                                                                                                 
cACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGtt                                                                                                                                               <  2:767929/151‑1 (MQ=255)
cACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGt                                                                                                                                                <  2:1059989/150‑1 (MQ=255)
   cGCACAAACGGAGATCGATTCCATGTCTTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGACGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGt                                                                                                                                            >  2:788422/1‑151 (MQ=255)
           cGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCGCCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACag                                                                                                                                    <  2:580490/151‑1 (MQ=255)
           cGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAAc                                                                                                                                      >  2:252923/1‑149 (MQ=255)
                   gATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGaaa                                                                                                                            <  1:311313/151‑1 (MQ=255)
                      tCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATc                                                                                                                          <  1:180418/150‑1 (MQ=255)
                        cATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGAc                                                                                                                       >  2:517363/1‑151 (MQ=255)
                          tGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGAccc                                                                                                                     >  2:641133/1‑151 (MQ=255)
                             aTTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCt                                                                                                                  <  2:573277/151‑1 (MQ=255)
                                tCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCttt                                                                                                                <  2:497536/150‑1 (MQ=255)
                                tCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTg                                                                                                               >  1:437735/1‑151 (MQ=255)
                                 ccGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTg                                                                                                               >  2:1109845/1‑150 (MQ=255)
                                  cGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGcc                                                                                                             <  2:821731/151‑1 (MQ=255)
                                    cTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTg                                                                                                               >  1:213306/1‑147 (MQ=255)
                                      tCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGcc                                                                                                             >  2:827317/1‑147 (MQ=255)
                                      tCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAg                                                                                                           <  2:385920/149‑1 (MQ=255)
                                       cccGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCGCCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATg                                                                                                        <  2:1126298/151‑1 (MQ=255)
                                       cccGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGAt                                                                                                         <  2:376614/150‑1 (MQ=255)
                                        ccGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGc                                                                                                       >  2:274532/1‑151 (MQ=255)
                                        ccGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGc                                                                                                       >  1:724223/1‑151 (MQ=255)
                                           ggCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgc                                                                                                    <  2:12930/151‑1 (MQ=255)
                                              gATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGgcagc                                                                                                 >  2:227136/1‑151 (MQ=255)
                                                  atcGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATa                                                                                              >  1:464491/1‑150 (MQ=255)
                                                    cGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATaa                                                                                             >  2:1194559/1‑149 (MQ=255)
                                                      ccAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATAc                                                                                         <  1:819023/151‑1 (MQ=255)
                                                       cAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATAcc                                                                                        >  2:223405/1‑151 (MQ=255)
                                                       cAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATAcc                                                                                        <  2:422113/151‑1 (MQ=255)
                                                       cAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATAcc                                                                                        >  1:459082/1‑151 (MQ=255)
                                                        agagCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCa                                                                                       <  2:1017215/151‑1 (MQ=255)
                                                          agCTCGTCACATTGCAGGCTGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAAAGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGCAATCGACCCGCATTGCCAGATGCGGCAGCATAAATACCAtt                                                                                     >  2:571175/1‑151 (MQ=255)
                                                           gCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTg                                                                                    <  2:558143/151‑1 (MQ=255)
                                                           gCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTg                                                                                    <  1:200741/151‑1 (MQ=255)
                                                           gCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTg                                                                                    <  1:1055916/151‑1 (MQ=255)
                                                              cGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGAc                                                                                  <  1:596484/150‑1 (MQ=255)
                                                                 aaCATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGCATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGc                                                                               >  2:749194/1‑150 (MQ=255)
                                                                       gCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCa                                                                        <  1:722302/151‑1 (MQ=255)
                                                                        cAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTcc                                                                         >  2:13463/1‑149 (MQ=255)
                                                                          ggATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCa                                                                     <  1:411482/151‑1 (MQ=255)
                                                                            aTGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCatg                                                                   <  2:314042/151‑1 (MQ=255)
                                                                                 gACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGcc                                                              >  1:885366/1‑151 (MQ=255)
                                                                                       cATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATAc                                                         <  2:1141118/150‑1 (MQ=255)
                                                                                          aaaTGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTgg                                                     >  2:173789/1‑151 (MQ=255)
                                                                                          aaaTGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTgg                                                     >  2:383688/1‑151 (MQ=255)
                                                                                              gTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGttt                                                  <  2:94965/150‑1 (MQ=255)
                                                                                                 cGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTggg                                               >  2:604778/1‑150 (MQ=255)
                                                                                                  gCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGa                                              >  2:155741/1‑150 (MQ=255)
                                                                                                   ccACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAAt                                            >  2:751145/1‑151 (MQ=255)
                                                                                                    cACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAAtt                                           >  2:939769/1‑151 (MQ=255)
                                                                                                    cACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAAtt                                           <  1:487331/151‑1 (MQ=255)
                                                                                                        cGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGt                                       >  1:632893/1‑151 (MQ=255)
                                                                                                            tATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAAcc                                   >  2:467812/1‑151 (MQ=255)
                                                                                                             aTTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCg                                  >  1:990717/1‑151 (MQ=255)
                                                                                                             aTTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCg                                  >  2:916274/1‑151 (MQ=255)
                                                                                                                    cGATCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCAtct                           <  2:1168235/151‑1 (MQ=255)
                                                                                                                      aTCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATctcc                         >  2:57975/1‑151 (MQ=255)
                                                                                                                      aTCAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATctcc                         <  1:179352/151‑1 (MQ=255)
                                                                                                                        cAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATctcctc                       >  2:514995/1‑151 (MQ=255)
                                                                                                                        cAACGGCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATctcct                        <  1:671258/150‑1 (MQ=255)
                                                                                                                             gCCAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACtt                  >  2:713228/1‑151 (MQ=255)
                                                                                                                               cAGAATGTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTc                 >  1:1090153/1‑150 (MQ=255)
                                                                                                                                     gTATTGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTcc           <  2:768363/150‑1 (MQ=255)
                                                                                                                                         tGAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCa       <  2:1181953/150‑1 (MQ=255)
                                                                                                                                          gAGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACg     >  1:760224/1‑151 (MQ=255)
                                                                                                                                           aGTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGt    >  2:750744/1‑151 (MQ=255)
                                                                                                                                            gTCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGtt   <  1:435901/151‑1 (MQ=255)
                                                                                                                                             tCAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTc  >  1:1967/1‑151 (MQ=255)
                                                                                                                                              cAGAGCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTc  >  2:460769/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CACCGCACAAACGGAGATCGATTCCATGTATTTCCGCTTCCCGGGCGATCATCGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTC  >  NC_000913/2966427‑2966718

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑