Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,724,771 T→A 5.1% Y21* (TAT→TAA) ‡ PP_3288 → universal stress protein family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,724,7710TA 5.1% 1062.1 / 31.4 314Y21* (TAT→TAA) ‡PP_3288universal stress protein family protein
Reads supporting (aligned to +/- strand):  ref base T (180/118);  new base A (2/14);  total (182/132)
Fisher's exact test for biased strand distribution p-value = 2.91e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 323 total aligned reads displayed.

GTAGGTGCCTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATG  >  NC_002947/3724622‑3724920
                                                                                                                                                     |                                                                                                                                                     
gTAGGTGCCTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTAt                                                                                                                                                       <  2:3161915/150‑1 (MQ=255)
    gTGCCTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGAcc                                                                                                                                                   <  2:5234302/150‑1 (MQ=255)
    gTGCCTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGAcc                                                                                                                                                   >  1:1859688/1‑150 (MQ=255)
    gTGCCTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGAcc                                                                                                                                                   >  1:831846/1‑150 (MQ=255)
    gTGCCTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGAcc                                                                                                                                                   <  2:2168597/150‑1 (MQ=255)
        cTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGgc                                                                                                                                               >  2:1695946/1‑150 (MQ=255)
        cTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCATAGACCgggg                                                                                                                                               <  1:7872412/150‑2 (MQ=255)
        cTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCATAGACCgggg                                                                                                                                               <  1:14529/150‑2 (MQ=255)
        cTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCATAGACCgggg                                                                                                                                               <  1:3419001/150‑2 (MQ=255)
         ttGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCATAGACCggggt                                                                                                                                              <  1:9020542/150‑3 (MQ=255)
              ttGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCtgctgc                                                                                                                                         <  2:8602486/150‑1 (MQ=255)
                gATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTgcgc                                                                                                                                       >  1:8448753/1‑150 (MQ=255)
                gATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTgcgc                                                                                                                                       >  2:2792428/1‑150 (MQ=255)
                  tAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTg                                                                                                                                     <  2:5110979/150‑1 (MQ=255)
                    aaCATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTggc                                                                                                                                   <  1:7556211/150‑1 (MQ=255)
                        tCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTggccaag                                                                                                                               <  2:9062947/150‑1 (MQ=255)
                         cAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTggccaagg                                                                                                                              >  2:5540781/1‑150 (MQ=255)
                         cAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTggccaagg                                                                                                                              <  1:63036/150‑1 (MQ=255)
                          aaGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTACTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTggccaaggc                                                                                                                             >  2:4268184/1‑150 (MQ=255)
                          aaGTTCGCAACGTGTAGCTGGCAGATGCTCAACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTggccaaggc                                                                                                                             >  1:3503651/1‑150 (MQ=255)
                                gCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAggg                                                                                                                       >  2:2938709/1‑150 (MQ=255)
                                 cAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGa                                                                                                                      >  2:716919/1‑150 (MQ=255)
                                        tAGCTGGCAGATGCTCTACCTTCGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGgcgc                                                                                                               <  2:1336569/150‑1 (MQ=255)
                                        tAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGACTGgcgc                                                                                                               >  2:9118969/1‑150 (MQ=255)
                                            tGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTggccaagg                                                                                                                              <  2:1151070/131‑1 (MQ=255)
                                              gCAGATGCTCTACCTACGTATCGCAACATCAGAGGTGTGCCATGGACACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTTATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGcaca                                                                                                         >  1:4513833/1‑150 (MQ=255)
                                                   tGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACAtcgtc                                                                                                    <  1:1249979/150‑1 (MQ=255)
                                                    gCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACAtcgtcg                                                                                                   >  1:7849916/1‑150 (MQ=255)
                                                        tacctacGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCttt                                                                                               >  1:4960662/1‑150 (MQ=255)
                                                        tacctacGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCttt                                                                                               >  2:8392024/1‑150 (MQ=255)
                                                        tacctacGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCttt                                                                                               <  2:7296407/150‑1 (MQ=255)
                                                         acctacGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGATTGGCGCTGCACATCGTCGCTTTc                                                                                              <  2:4677334/150‑1 (MQ=255)
                                                         acctacGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTc                                                                                              <  2:6241108/150‑1 (MQ=255)
                                                         acctacGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTc                                                                                              <  2:5429369/150‑1 (MQ=255)
                                                                  tCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTg                                                                                     >  2:3037236/1‑150 (MQ=255)
                                                                  tCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTg                                                                                     >  1:2760773/1‑150 (MQ=255)
                                                                  tCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTg                                                                                     >  1:7056991/1‑150 (MQ=255)
                                                                   cGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTgg                                                                                    <  1:7666888/150‑1 (MQ=255)
                                                                     cAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTTCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTTCACATCGTCGCTTTCGACTACGTGGaa                                                                                  <  2:8729857/150‑1 (MQ=255)
                                                                      aaCCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAg                                                                                 <  1:2539314/150‑1 (MQ=255)
                                                                         cTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTg                                                                                                                                     <  2:6600490/95‑1 (MQ=255)
                                                                         cTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGc                                                                              >  1:5666048/1‑150 (MQ=255)
                                                                         cTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGAACATCGTCGCTTTCGACTACGTGGAAGGGc                                                                              >  1:3565452/1‑150 (MQ=255)
                                                                          tCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCt                                                                             >  1:5553298/1‑150 (MQ=255)
                                                                            agagGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGCTCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTg                                                                           <  2:7802495/150‑1 (MQ=255)
                                                                            agagGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTg                                                                           <  2:6147226/150‑1 (MQ=255)
                                                                                   tgCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCAcggc                                                                    <  1:6922547/150‑1 (MQ=255)
                                                                                   tgCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCAcggc                                                                    >  2:8200115/1‑150 (MQ=255)
                                                                                     ccATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCAcggcgg                                                                  >  2:4569492/1‑150 (MQ=255)
                                                                                     ccATGGAAACTCCGCAAAGACTGCTGATGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCAcggcgg                                                                  >  1:4816712/1‑150 (MQ=255)
                                                                                      cATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCggg                                                                 >  2:520337/1‑150 (MQ=255)
                                                                                       aTGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGc                                                                >  2:7542350/1‑150 (MQ=255)
                                                                                       aTGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGc                                                                >  2:3521189/1‑150 (MQ=255)
                                                                                       aTGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGc                                                                >  2:2473144/1‑150 (MQ=255)
                                                                                        tGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTTATGCGCCGAACCCCGGTCTATGACCGGGCTGCAGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGcc                                                               >  2:6800164/1‑150 (MQ=255)
                                                                                        tGGAAACTACGCAAAGACTGCTGCTGATCGCCTCCCACCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGCGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGcc                                                               >  2:8907635/1‑150 (MQ=255)
                                                                                          gAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCAAGGCGGGCCTg                                                             >  2:7998233/1‑150 (MQ=255)
                                                                                            aaCTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGt                                                           <  2:4088058/150‑1 (MQ=255)
                                                                                               tCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGaa                                                        <  2:8673838/150‑1 (MQ=255)
                                                                                               tCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGaa                                                        <  1:2399594/150‑1 (MQ=255)
                                                                                                 cGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATg                                                      >  1:1344396/1‑150 (MQ=255)
                                                                                                    aaaGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGAcc                                                   >  1:6770455/1‑150 (MQ=255)
                                                                                                        aCTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGgc                                               <  1:1046504/150‑1 (MQ=255)
                                                                                                            ctgctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTg                                           <  1:1183798/150‑1 (MQ=255)
                                                                                                             tgctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTgg                                          >  1:3832310/1‑150 (MQ=255)
                                                                                                               ctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCa                                        <  2:120536/150‑1 (MQ=255)
                                                                                                                tgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAg                                       >  1:7508354/1‑150 (MQ=255)
                                                                                                                 gATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGt                                      >  1:2990078/1‑150 (MQ=255)
                                                                                                                   tCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCa                                    >  1:2660022/1‑150 (MQ=255)
                                                                                                                   tCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCa                                    >  1:2659881/1‑150 (MQ=255)
                                                                                                                   tCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCa                                    >  1:651322/1‑150 (MQ=255)
                                                                                                                     gCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATg                                  >  1:6405926/1‑150 (MQ=255)
                                                                                                                      ccTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATgc                                 >  1:7101018/1‑150 (MQ=255)
                                                                                                                      ccTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATgc                                 >  2:6498822/1‑150 (MQ=255)
                                                                                                                        tCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATgcgc                               >  2:7700492/1‑150 (MQ=255)
                                                                                                                          cccccTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGa                             <  2:3711558/150‑1 (MQ=255)
                                                                                                                          cccccTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGa                             >  1:2209999/1‑150 (MQ=255)
                                                                                                                          cccccTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGa                             >  1:6264363/1‑150 (MQ=255)
                                                                                                                           ccccTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAg                            >  1:523085/1‑150 (MQ=255)
                                                                                                                            cccTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAgg                           >  1:3017408/1‑150 (MQ=255)
                                                                                                                            cccTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAgg                           >  1:5351781/1‑150 (MQ=255)
                                                                                                                            cccTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAgg                           <  2:8261272/150‑1 (MQ=255)
                                                                                                                               tGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGtt                        <  1:5850170/150‑1 (MQ=255)
                                                                                                                                gATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTa                       >  2:4996440/1‑150 (MQ=255)
                                                                                                                                  tGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACg                     <  1:7420224/150‑1 (MQ=255)
                                                                                                                                     gcCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGACTGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCCCGGCGGGCCTGGTGACTGACCAGGCGCTGGCAGTCCTGCGCGAGGGTTACGTgc                  >  1:956919/1‑150 (MQ=255)
                                                                                                                                      ccGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTgca                 <  1:3312641/150‑1 (MQ=255)
                                                                                                                                      ccGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTgca                 >  2:1968216/1‑150 (MQ=255)
                                                                                                                                      ccGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTgca                 >  2:4760203/1‑150 (MQ=255)
                                                                                                                                       agcagcccGGTCATAGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTgcag                >  1:4144208/6‑150 (MQ=255)
                                                                                                                                        gAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTgcagc               <  1:1626193/150‑1 (MQ=255)
                                                                                                                                        gAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTgcagc               <  1:1728014/150‑1 (MQ=255)
                                                                                                                                           ccccGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTgcagcagc            >  1:8487345/1‑150 (MQ=255)
                                                                                                                                            cccGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTgcagcagca           >  2:6184340/1‑150 (MQ=255)
                                                                                                                                            cccGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTgcagcagca           >  1:8886876/1‑150 (MQ=255)
                                                                                                                                              cGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCAcc         <  1:6880228/150‑1 (MQ=255)
                                                                                                                                              cGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCAcc         >  1:7915610/1‑150 (MQ=255)
                                                                                                                                              cGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCAcc         >  1:8108735/1‑150 (MQ=255)
                                                                                                                                                gTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGt       >  1:2687469/1‑150 (MQ=255)
                                                                                                                                                     tGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATg  >  1:2048991/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTAGGTGCCTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATG  >  NC_002947/3724622‑3724920

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