Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,033,816 T→A 8.4% F145I (TTC→ATC)  PP_0891 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,033,8160TA 8.4% 288.8 / 14.6 97F145I (TTC→ATC) PP_0891hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (45/42);  new base A (3/5);  total (49/48)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

AAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCTTCACTGCGCGCTG  >  NC_002947/1033672‑1033964
                                                                                                                                                  |                                                                                                                                                     
aaGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAg                                                                                                                                                  <  1:1876541/150‑1 (MQ=255)
 aGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCAACTTGCATGGCCAGTCGCACAGCTACAATTTCCTGTTGCTGGCCACCCGCCTGGA‑AGCGCC‑GCGCATCCAgg                                                                                                                                                 <  1:2124062/150‑1 (MQ=255)
  ggCACAGATCGACCTGCAGCAGATCAGTGCCCTTGCCCTCGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACTCACCTGGA‑AGCGCC‑GCGCATCCAggc                                                                                                                                                <  2:1203605/150‑2 (MQ=255)
     acaGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAggt‑ggt                                                                                                                                            >  1:1176721/1‑150 (MQ=255)
     acaGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAggt                                                                                                                                                <  2:892460/147‑1 (MQ=255)
     acaGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAggt                                                                                                                                                >  1:892460/1‑147 (MQ=255)
      caGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGTCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑AGCGCC‑GCGCATCCAggcgggt                                                                                                                                            <  2:1943575/150‑6 (MQ=255)
       aGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGa                                                                                                                                          >  2:2388392/1‑150 (MQ=255)
        gATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGAc                                                                                                                                         <  2:1446558/150‑1 (MQ=255)
        gATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGAc                                                                                                                                         <  2:1417259/150‑1 (MQ=255)
            gACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGG‑CAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTg                                                                                                                                    <  2:185181/150‑1 (MQ=255)
                tGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGc                                                                                                                                 >  1:49169/1‑150 (MQ=255)
                tGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGc                                                                                                                                 >  2:67925/1‑150 (MQ=255)
                   aGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCctc                                                                                                                              >  2:285328/1‑150 (MQ=255)
                   aGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCctc                                                                                                                              <  2:2141330/150‑1 (MQ=255)
                    gCCTATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCctct                                                                                                                             >  2:1603562/1‑150 (MQ=255)
                      cGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTgg                                                                                                                           <  1:1275807/150‑1 (MQ=255)
                             tGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATATCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAggt                                                                                                                                                >  1:1810479/1‑123 (MQ=255)
                             tGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATATCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAggt                                                                                                                                                <  2:1810479/123‑1 (MQ=255)
                                    gCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGaa                                                                                                             <  1:1603562/150‑1 (MQ=255)
                                         cGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGAttt                                                                                                        <  1:2372426/150‑1 (MQ=255)
                                           gTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTg                                                                                                      >  2:2213438/1‑150 (MQ=255)
                                           gTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTg                                                                                                      <  2:2398109/150‑1 (MQ=255)
                                            tGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTgg                                                                                                     >  2:1617862/1‑150 (MQ=255)
                                             gCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGc                                                                                                    <  1:285328/150‑1 (MQ=255)
                                               acaGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCtt                                                                                                  >  2:2004150/1‑150 (MQ=255)
                                               acaGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑AGCGCC‑GCGCATCCAggcgg                                                                                                                                              >  2:1777603/1‑104 (MQ=255)
                                               acaGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑AGCGCC‑GCGCATCCAggc‑gg                                                                                                                                             <  1:1777603/107‑1 (MQ=255)
                                                 aGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTtgc                                                                                                <  1:1227307/150‑1 (MQ=255)
                                                   aTCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCtt                                                                                                                                     <  1:1509782/111‑1 (MQ=255)
                                                   aTCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCtt                                                                                                                                     >  2:1509782/1‑111 (MQ=255)
                                                    tagAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTtgctgc                                                                                             <  1:2072513/148‑1 (MQ=255)
                                                         ttGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCgggg                                                                                        >  2:519331/1‑150 (MQ=255)
                                                           gcgcCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCt                                                                                      <  2:1120171/150‑1 (MQ=255)
                                                                tGCCGCTGCCGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGGGAACCG‑‑GGGAGTGGGTGTTCCCTTCCAGGT‑GTTGACCATTGAGCCTCTGGTGCTGCCCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGaaa                                                                                 <  1:2388392/150‑1 (MQ=255)
                                                                  ccGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGTGTCTAAGGGTCGGCGCGCCGAATATTTTGGCTTGCTGCCGGGGCTGGAAAGc                                                                               <  1:1617862/150‑1 (MQ=255)
                                                                   cGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGcc                                                                              >  1:8520/1‑150 (MQ=255)
                                                                    gCTGACGGTCACCCTGCATGGCCAGTCGCAGAGCGACAATGTCCTGTTGCTGGCCACCCGCCTGGC‑TGGGCG‑GCACTTCTAGGT‑GGTGAACCTTGCGCTTCTGGTGCTACGAGCCGAAGATTTTGGGTTGATTGCGCGGCTGGCATGCCt                                                                             >  2:1141218/1‑150 (MQ=255)
                                                                    gCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGaaa                                                                                 <  1:296276/146‑1 (MQ=255)
                                                                    gCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGaaa                                                                                 >  2:296276/1‑146 (MQ=255)
                                                                    gCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCt                                                                             <  2:713373/150‑1 (MQ=255)
                                                                      tGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTgc                                                                           <  1:1162777/150‑1 (MQ=255)
                                                                          ggTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCaa                                                                       <  1:1562919/150‑1 (MQ=255)
                                                                          ggTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCaa                                                                       <  1:67925/150‑1 (MQ=255)
                                                                            tCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAGGCCTGCGCAAGt                                                                     >  2:417401/1‑150 (MQ=255)
                                                                              aCCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTc                                                                   >  2:345606/1‑150 (MQ=255)
                                                                              aCCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTc                                                                   >  2:2056908/1‑150 (MQ=255)
                                                                               cccTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGGTGCGGCCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCg                                                                  >  2:1244769/1‑150 (MQ=255)
                                                                                ccTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGc                                                                 >  2:777814/1‑150 (MQ=255)
                                                                                ccTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGc                                                                 <  1:2177248/150‑1 (MQ=255)
                                                                                ccTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGc                                                                 >  2:1826561/1‑150 (MQ=255)
                                                                                ccTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGc                                                                 >  2:1854301/1‑150 (MQ=255)
                                                                                 cTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGcc                                                                >  1:922297/1‑150 (MQ=255)
                                                                                    cATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAgg                                                             >  2:1974833/1‑150 (MQ=255)
                                                                                     aTGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGt                                                            <  1:1826561/150‑1 (MQ=255)
                                                                                         ccAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGTGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGaa                                                        >  1:160211/1‑150 (MQ=255)
                                                                                          cAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAg                                                       <  1:1832983/150‑1 (MQ=255)
                                                                                           aGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGt                                                      >  1:358560/1‑150 (MQ=255)
                                                                                               gCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCAt                                                  <  2:1176721/150‑1 (MQ=255)
                                                                                                 acaGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCa                                                >  1:1148348/1‑150 (MQ=255)
                                                                                                   aGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAAc                                              <  2:677993/150‑1 (MQ=255)
                                                                                                    gCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAAcc                                             <  1:206861/150‑1 (MQ=255)
                                                                                                     cTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCa                                                               >  2:1363061/1‑131 (MQ=255)
                                                                                                     cTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCa                                                               <  1:1363061/131‑1 (MQ=255)
                                                                                                     cTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAAccc                                            <  2:160211/150‑1 (MQ=255)
                                                                                                      tACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCg                                           >  1:1527196/1‑150 (MQ=255)
                                                                                                      tACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCg                                           >  2:2118378/1‑150 (MQ=255)
                                                                                                       aCAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGt                                          <  1:777814/150‑1 (MQ=255)
                                                                                                       aCAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGt                                          <  2:49169/150‑1 (MQ=255)
                                                                                                        cAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTc                                         >  1:690291/1‑150 (MQ=255)
                                                                                                           tGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGACTTTAAATCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGt                                      >  1:2090924/1‑150 (MQ=255)
                                                                                                           gacattGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCAt                                                  >  2:803015/6‑138 (MQ=255)
                                                                                                           gacattGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCAt                                                  <  1:803015/133‑1 (MQ=255)
                                                                                                             tCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAAc                                              <  2:467940/140‑1 (MQ=255)
                                                                                                             tCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAAc                                              >  1:467940/1‑140 (MQ=255)
                                                                                                               cTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCCTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCg                                  >  2:226752/1‑150 (MQ=255)
                                                                                                                         ccACCCGCCTGGA‑TGCGCG‑GCGCCTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCgg                        <  1:226752/150‑1 (MQ=255)
                                                                                                                          cACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTTCGGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCggg                       >  1:159480/1‑149 (MQ=255)
                                                                                                                          cACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGt                       >  1:2384137/1‑150 (MQ=255)
                                                                                                                          cACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCAGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGcc                                                                              >  1:733818/1‑95 (MQ=255)
                                                                                                                          cACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCAGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGcc                                                                              <  2:733818/95‑1 (MQ=255)
                                                                                                                          cACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTCGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTTCATCAACCCGTCGGTACCGGTTAGCGCGGt                       <  2:1409629/150‑1 (MQ=255)
                                                                                                                           aCCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTg                      >  1:253327/1‑150 (MQ=255)
                                                                                                                             ccGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTg                                                                                                      >  2:2096261/1‑68 (MQ=255)
                                                                                                                             ccGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTg                                                                                                      <  1:2096261/68‑1 (MQ=255)
                                                                                                                             ccGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCt                    >  2:1481331/1‑150 (MQ=255)
                                                                                                                              caccTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTg                   >  1:1709166/3‑150 (MQ=255)
                                                                                                                              cGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTg                   <  1:2004150/150‑1 (MQ=255)
                                                                                                                                 cTGGA‑AGCGCC‑GCGCATCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATc                >  2:1219275/1‑150 (MQ=255)
                                                                                                                                  tGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGATAAAAGTAATGATAGTTGTTCTCTGTGTCTTGCTGCCGGTGCTGGAAAGCCTGCGCAAGTTCGTCAGTTTGAAGTCCATCAACCCGTCGGTGCCGGTTAGCGCGGTGCTGATCt               >  1:1225886/1‑150 (MQ=255)
                                                                                                                                     a‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCTTCa            >  2:68324/1‑150 (MQ=255)
                                                                                                                                        gcgcC‑GCGCATCCAGGCGGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCTTCAc           >  1:1257203/1‑150 (MQ=255)
                                                                                                                                          gcg‑gcgCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTgc                                                                           >  1:854892/1‑83 (MQ=255)
                                                                                                                                          gcg‑gcgCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTgc                                                                           <  2:854892/83‑1 (MQ=255)
                                                                                                                                           gccgcgcATCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCTTCACTgc        <  1:1219275/147‑1 (MQ=255)
                                                                                                                                                gcTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCTTCACTGCGCGCt   <  1:560158/150‑1 (MQ=255)
                                                                                                                                                 cTTCCAGGT‑GGTGTCCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGGAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCTTCACTGCGCGCTg  <  2:1744126/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
AAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCTTCACTGCGCGCTG  >  NC_002947/1033672‑1033964

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