Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,048,585 C→A 7.1% S199S (TCC→TCA)  PP_2661 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,048,5850CA 7.1% 259.6 / 6.9 84S199S (TCC→TCA) PP_2661hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (36/42);  new base A (2/4);  total (38/46)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.89e-01

GTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGG‑GCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAA  >  NC_002947/3048437‑3048734
                                                                                                                                                    |                                                                                                                                                      
gTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCt                                                                                                                                                       <  1:2322025/150‑1 (MQ=255)
gTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCt                                                                                                                                                       <  1:2209782/150‑1 (MQ=255)
 tttCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTg                                                                                                                                                      >  2:288021/1‑150 (MQ=255)
   tCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGaa                                                                                                                                                    <  2:1119285/150‑1 (MQ=255)
      gcgTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGcc                                                                                                                                                 >  2:1997444/1‑150 (MQ=255)
      gcgTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGcc                                                                                                                                                 >  2:1889037/1‑150 (MQ=255)
      gcgTGGTGGCGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCCCGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGcc                                                                                                                                                 >  2:747321/1‑150 (MQ=255)
        gtggtggtggtgTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAg                                                                                                                                               <  2:2177889/150‑1 (MQ=255)
        gtggtggtggtgTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAg                                                                                                                                               <  1:1336472/150‑1 (MQ=255)
         tggtggtggtgTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAgg                                                                                                                                              >  2:1800278/1‑150 (MQ=255)
         tggtggtggtgTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAgg                                                                                                                                              <  2:2057859/150‑1 (MQ=255)
         gggtggtggtgTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCAGGAAGCCAgt                                                                                                                                              <  1:990619/149‑2 (MQ=255)
          ggtggtggtgTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGga                                                                                                                                             <  1:1463036/150‑1 (MQ=255)
          ggtggtggtgTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCAGGAAGCCAgtg                                                                                                                                             <  2:1260454/150‑3 (MQ=255)
            tggtggtgTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGgaga                                                                                                                                           >  1:1343403/1‑150 (MQ=255)
            tggtggtgTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGgaga                                                                                                                                           >  2:375124/1‑150 (MQ=255)
             ggtggtgTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGAt                                                                                                                                          <  1:1199630/150‑1 (MQ=255)
                   gtgCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCt                                                                                                                                    <  2:1871218/150‑1 (MQ=255)
                    tgCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTc                                                                                                                                   <  1:826156/150‑1 (MQ=255)
                    tgCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGAGTTAGCAGTCACTGGCTTCCTGAAGCCAGGAGATTCCGCTc                                                                                                                                   >  2:2018391/1‑150 (MQ=255)
                      cccGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAg                                                                                                                                 >  2:350340/1‑150 (MQ=255)
                      cccGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGGGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCTTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAg                                                                                                                                 >  2:104852/1‑150 (MQ=255)
                       ccGAAGGATGAGCACTGGGCGCATAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGa                                                                                                                                >  2:1352270/1‑150 (MQ=255)
                       ccGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGa                                                                                                                                >  1:2086497/1‑150 (MQ=255)
                        cGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGAt                                                                                                                               <  1:288021/150‑1 (MQ=255)
                         gAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCGCGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATa                                                                                                                              >  1:1866385/1‑150 (MQ=255)
                         gAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATa                                                                                                                              >  1:1893637/1‑150 (MQ=255)
                         gAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATa                                                                                                                              >  1:2143068/1‑150 (MQ=255)
                            ggATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCCTCCTGAAGCCAGGAGATTCCGCTCAGATAATg                                                                                                                           <  2:1345088/150‑1 (MQ=255)
                             gATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGt                                                                                                                          <  1:375124/150‑1 (MQ=255)
                                    aCTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTggg                                                                                                                   <  2:1567316/150‑1 (MQ=255)
                                           gcAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCtt                                                                                                            <  2:2069378/150‑1 (MQ=255)
                                              gAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGc                                                                                                         <  1:106112/150‑1 (MQ=255)
                                               aaGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCa                                                                                                        <  1:1311405/150‑1 (MQ=255)
                                                 ggTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGc                                                                                                      >  1:2448362/1‑150 (MQ=255)
                                                  gTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGGTAGGGTTCGCAGGCACTGGCTTCCTGATGAAGAGTGATCGTGCGCGTATGATGTAATTGGGGCTGCTTTGCAGcc                                                                                                     >  1:1740068/1‑150 (MQ=255)
                                                    gAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCt                                                                                                   >  1:1671866/1‑150 (MQ=255)
                                                     aaGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCTTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCAGGAAGCCAgt                                                                                                                                              <  1:2124125/106‑2 (MQ=255)
                                                     aaGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCAGGAAGCCAgt                                                                                                                                              >  2:2124125/1‑105 (MQ=255)
                                                          cgcgGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCgg                                                                                             <  1:1234986/150‑1 (MQ=255)
                                                           gcgGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCggg                                                                                            <  1:1997444/150‑1 (MQ=255)
                                                            cgGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGc                                                                                           <  1:713402/150‑1 (MQ=255)
                                                              gtggtgGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTCATTGGGGCTGCTTTGCGGCCCTTCGCGGGctc                                                                                         <  2:1343403/150‑3 (MQ=255)
                                                                  tggATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGccc                                                                                     <  1:165336/150‑1 (MQ=255)
                                                                    gATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCcgc                                                                                   >  1:1755392/1‑150 (MQ=255)
                                                                          cAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAgg‑g                                                                            >  2:1159239/1‑150 (MQ=255)
                                                                          cAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAgg‑g                                                                            <  1:1889037/150‑1 (MQ=255)
                                                                           aGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGG‑Gc                                                                           <  2:2448362/150‑1 (MQ=255)
                                                                                   ctggcctggcAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGG‑GCCGCAAggc                                                                   <  2:413246/150‑1 (MQ=255)
                                                                                       cctggcAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGG‑GCCGCAAGGCGGcc                                                               <  2:1740068/150‑1 (MQ=255)
                                                                                        ctggcAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCGGCTCTTCGCGGGCTCGCCCGCGAAGG‑GCTGCaaa                                                                     >  1:725881/1‑142 (MQ=38)
                                                                                        ctggcAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCGGCCCTTCGCGGGCTCGCCCGCGAAGG‑GCTGCaaa                                                                     <  2:725881/143‑2 (MQ=38)
                                                                                             agcaACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGG‑GCCGCAAGGCGGCCCCGGtt                                                         >  1:1479696/1‑150 (MQ=255)
                                                                                                aaCCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGG‑GCCGCAAGGCGGcccccgggcat                                                      >  1:2434204/1‑143 (MQ=255)
                                                                                                 aCCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGG‑GCCGCAAGGCGGCCCCGGttcatt                                                     >  2:1432633/1‑150 (MQ=255)
                                                                                                 aCCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGG‑GCCGCAAGGCGGCCCCGGttcatt                                                     >  2:1682786/1‑150 (MQ=255)
                                                                                                  ccGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCGCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGG‑GCCGCAAGGCGGCCCCGGttcattc                                                    <  2:1893637/150‑1 (MQ=255)
                                                                                                   cGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGG‑GCCGCAAGGCGGCCCCGGttcattca                                                   >  2:1866196/1‑150 (MQ=255)
                                                                                                           tATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGG‑GCCGCAAGGCGGCCCCGGTTCATTCAATCGTAtt                                           <  1:497304/150‑1 (MQ=255)
                                                                                                                aaTGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGG‑GCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATAc                                      <  1:1215340/150‑1 (MQ=255)
                                                                                                                aaTGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGG‑GCCCCAAGGCGGCCCCGGTTCATTCAATCGTATTGATAc                                      <  1:2018391/150‑1 (MQ=255)
                                                                                                                aaTGCCCGCTTGCCAGCGTTCGCAGGCACGGGCGTCCTGAAGCCAGGAGATTCCGCTCAGATACTGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGG‑GCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATAc                                      <  2:1756825/150‑1 (MQ=255)
                                                                                                                   gCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGG‑GCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCg                                   >  1:175273/1‑150 (MQ=255)
                                                                                                                     ccGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGG‑GCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGcc                                 <  2:2086497/150‑1 (MQ=255)
                                                                                                                     ccGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGG‑GCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGcc                                 >  1:1951919/1‑150 (MQ=255)
                                                                                                                     ccGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCAAGATC‑AGACGGCTCTTCAAGTCTGCGCCGGTTCATTCAATCGTATTGATACCCGcc                                 >  1:1116280/1‑150 (MQ=255)
                                                                                                                      cGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGG‑GCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCt                                <  1:350340/150‑1 (MQ=255)
                                                                                                                       cagtGCCTGCGAACGCTGGCACTGGCTTCCTGGAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGG‑GCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTg                               >  2:1349677/4‑150 (MQ=255)
                                                                                                                       cagtGCCTGCGAACGCTGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGaa                                                                                <  1:988350/99‑1 (MQ=38)
                                                                                                                              aGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGG‑GCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTcc                        >  1:2433988/1‑150 (MQ=255)
                                                                                                                                  ttCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGc                                                                                           >  2:1047226/1‑80 (MQ=255)
                                                                                                                                  ttCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGc                                                                                           <  1:1047226/80‑1 (MQ=255)
                                                                                                                                     gCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGAGAGCCCGCGAAGG‑GCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGa                 <  1:1352270/150‑1 (MQ=255)
                                                                                                                                      cAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGGCCTTCGCGGGCGAGCCCGCGAAgg                                                                              <  2:1623951/89‑1 (MQ=255)
                                                                                                                                      cAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAgg                                                                              >  1:1623951/1‑89 (MQ=255)
                                                                                                                                      cAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGG‑GCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAg                <  2:175273/150‑1 (MQ=255)
                                                                                                                                       aGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGG‑GCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGc               >  2:747437/1‑150 (MQ=255)
                                                                                                                                         gCACTGGCTTCCTGAAGCCATGAGATGCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCCTTCCGATCTAGCCCGCGAAGG‑GCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCa             <  1:747321/150‑1 (MQ=255)
                                                                                                                                           ccTGGCTTCAGGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGG‑GCCGCAAGGCGGCCCCGGTTCATTCAAt                                                 <  2:252791/111‑1 (MQ=255)
                                                                                                                                           ccTGGCTTCAGGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGG‑GCCGCAAGGCGGCCCCGGTTCATTCAAt                                                 >  1:252791/2‑112 (MQ=255)
                                                                                                                                            cTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGG‑GCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATaa          <  2:1755392/150‑1 (MQ=255)
                                                                                                                                            cTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGG‑GCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATaa          <  2:1479696/150‑1 (MQ=255)
                                                                                                                                               gCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGG‑GCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGtc       >  1:1526459/1‑150 (MQ=255)
                                                                                                                                               gCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGG‑GCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGtc       >  1:1352818/1‑150 (MQ=255)
                                                                                                                                                   ccTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTCCGCGGGCGAGCCCGCGAAGG‑GCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGa   <  1:104852/150‑1 (MQ=255)
                                                                                                                                                    cTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGG‑GCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGaa  >  2:1906034/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGG‑GCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAA  >  NC_002947/3048437‑3048734

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: ‑