Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 508,957:1 +C 8.6% coding (57/1890 nt) estP ← esterase EstP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947508,9571.C 8.6% 262.3 / 20.1 70coding (57/1890 nt)estPesterase EstP
Reads supporting (aligned to +/- strand):  ref base . (33/31);  new base C (2/4);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.36e-01

GACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGCTCCCAAGTCCATTCGTTTCCCCGT  >  NC_002947/508808‑509106
                                                                                                                                                      |                                                                                                                                                     
gACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTgg                                                                                                                                                        <  1:1074192/150‑1 (MQ=255)
gACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTgg                                                                                                                                                        <  1:91023/150‑1 (MQ=255)
gACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTgg                                                                                                                                                        <  2:1337018/150‑1 (MQ=255)
      ggggCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCa                                                                                                                                                  >  1:1026597/1‑149 (MQ=255)
      ggggCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCa                                                                                                                                                  <  2:1026597/149‑1 (MQ=255)
      ggggCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAg                                                                                                                                                 >  2:314561/1‑150 (MQ=255)
              gTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAggccagggc                                                                                                                                         <  1:1010012/150‑1 (MQ=255)
               ttGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGGCCTGCa‑gccaggcgt                                                                                                                                        <  2:1587362/150‑4 (MQ=255)
                tGCCGTCGGCGTCACGGTTGGTGAACCGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAggccagggcca                                                                                                                                       <  1:4701/150‑1 (MQ=255)
                tGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGGCCTGCa‑gccaggcgtt                                                                                                                                       <  2:1071289/150‑5 (MQ=255)
                  ccGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGAGCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGt                                                                                                                                     <  1:1908188/150‑1 (MQ=255)
                   cGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGATGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTg                                                                                                                                    <  1:487247/150‑1 (MQ=255)
                     tCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGaa                                                                                                                                  <  2:1714035/150‑1 (MQ=255)
                       ggCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGc                                                                                                                                >  1:1159560/1‑150 (MQ=255)
                        gCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCg                                                                                                                               >  2:524610/1‑150 (MQ=255)
                         cGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAAATGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGa                                                                                                                              <  1:245428/150‑1 (MQ=255)
                          gTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAg                                                                                                                             <  2:1599138/150‑1 (MQ=255)
                           tCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAgg                                                                                                                            <  1:735173/150‑1 (MQ=255)
                                gTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTaaa                                                                                                                       <  1:1786231/150‑1 (MQ=255)
                                gTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTaaa                                                                                                                       <  1:314561/150‑1 (MQ=255)
                                     tGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCa                                                                                                                  >  1:1794814/1‑150 (MQ=255)
                                      gAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCaa                                                                                                                 >  2:382815/1‑150 (MQ=255)
                                          cgcgcgCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAAc                                                                                                             <  1:96125/150‑1 (MQ=255)
                                           gcgcgcCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACg                                                                                                            >  1:418548/1‑150 (MQ=255)
                                           gcgcgcCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACg                                                                                                            <  1:951007/150‑1 (MQ=255)
                                             gcgcCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACggg                                                                                                          <  1:657140/150‑1 (MQ=255)
                                               gcCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGc                                                                                                        >  1:538112/1‑150 (MQ=255)
                                                            ccAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGGCCTGCa‑gccaggc                                                                                                                                          <  1:137563/103‑2 (MQ=255)
                                                            ccAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGGCCTGCa‑gccaggc                                                                                                                                          >  2:137563/1‑102 (MQ=255)
                                                            ccAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCggg                                                                                           <  1:1871166/150‑1 (MQ=255)
                                                              aacaaGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGGGGGTAGAGCGCAATAACGGGGCTTTTCGCATCGGGGa                                                                                         >  1:1066133/1‑150 (MQ=255)
                                                                 aaGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCACCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACa                                                                                      <  1:834918/150‑1 (MQ=255)
                                                                     tCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCt                                                                                  >  1:293184/1‑150 (MQ=255)
                                                                        ggAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCtttt                                                                               >  1:1621636/1‑150 (MQ=255)
                                                                               gCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGtt                                                                        <  2:1060815/150‑1 (MQ=255)
                                                                                cccGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGGCCTGCa‑gccaggcgtt                                                                                                                                       <  1:108171/86‑5 (MQ=255)
                                                                                cccGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGGCCTGCa‑gccaggcgtt                                                                                                                                       >  2:108171/1‑82 (MQ=255)
                                                                                  cGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTc                                                                                                   >  1:1651990/1‑120 (MQ=255)
                                                                                  cGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTc                                                                                                   <  2:1651990/120‑1 (MQ=255)
                                                                                    gCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTAt                                                                   >  2:227943/1‑150 (MQ=255)
                                                                                      aTCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCa                                                                 >  1:1801649/1‑150 (MQ=255)
                                                                                       tCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAg                                                                >  2:1404048/1‑150 (MQ=255)
                                                                                        cGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAgg                                                               <  1:828640/150‑1 (MQ=255)
                                                                                          gCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGc                                                             >  2:1075482/1‑150 (MQ=255)
                                                                                            gAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCaa                                                           <  2:538112/150‑1 (MQ=255)
                                                                                              ggCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCCACCCCTGG‑CTGCCGGCCAGGGCCAGGGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGc                                                         <  2:688110/150‑1 (MQ=255)
                                                                                              ggCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCATGGCAAGc                                                         <  1:382815/150‑1 (MQ=255)
                                                                                                 tGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCa                                                      >  1:289674/1‑150 (MQ=255)
                                                                                                   tCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATc                                                    <  1:518576/150‑1 (MQ=255)
                                                                                                   tCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATc                                                    <  1:54747/150‑1 (MQ=255)
                                                                                                   tCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGAGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATc                                                    >  1:1015348/1‑150 (MQ=255)
                                                                                                      ccAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGAc                                                 >  1:109770/1‑150 (MQ=255)
                                                                                                       cAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGAcc                                                >  2:262440/1‑150 (MQ=255)
                                                                                                                tCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTGTTCGCATCGGGGAACATCCTTTTCTTTGTGTTTATCAGGGCGAGCCCATCGACCTTAGCAAAt                                       >  1:1531783/1‑150 (MQ=255)
                                                                                                                 cAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGttt                                                                       <  1:895100/117‑1 (MQ=255)
                                                                                                                 cAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGttt                                                                       >  2:895100/1‑117 (MQ=255)
                                                                                                                 cAGGGTTGAATAGGGCGAAGGGCCGGCCAACGACAGT‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGGCCTTAGCAAAtc                                      >  1:912586/1‑150 (MQ=255)
                                                                                                                    ggTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAAtcgtc                                   <  1:1075482/150‑1 (MQ=255)
                                                                                                                      ttGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATcgtcgt                                 >  1:354352/1‑150 (MQ=255)
                                                                                                                      ttGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCAAGGTCAACACGGATCAACTGGGCTTTTCGCATCTGGGAACATCCTTTTCTTTGTTTTTCTCAGTGCAAGCCCATCGACCTTAGCAAATcgtcgt                                 <  1:135240/150‑1 (MQ=255)
                                                                                                                         aaTAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCt                              >  1:1434781/1‑150 (MQ=255)
                                                                                                                          aTAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCtt                             >  1:260917/1‑150 (MQ=255)
                                                                                                                             gggCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGc                          >  2:653425/1‑150 (MQ=255)
                                                                                                                              ggCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAAtc                                      <  1:979667/137‑1 (MQ=255)
                                                                                                                              ggCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGACTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAAtc                                      >  2:979667/1‑137 (MQ=255)
                                                                                                                                 gAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCATTTTCTTCGGTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGCTccc                      >  1:416450/1‑150 (MQ=255)
                                                                                                                                   aGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGACGGTAAAGCGCAATAACGGTGCTGTTCGCGTCGTGGAACATCCTTTTCTTTGGTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGCTCCCaa                    <  2:1159560/150‑1 (MQ=255)
                                                                                                                                     ggccggccAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGCTCCCAAGt                  <  2:289674/150‑1 (MQ=255)
                                                                                                                                           ccAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGCTCCCAAGTCCATTc            <  1:1666746/150‑1 (MQ=255)
                                                                                                                                             aaCGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGCTCCCAAGTCCATTCGt          >  2:276353/1‑150 (MQ=255)
                                                                                                                                               cGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGCTCCCAAGTCCATTCGttt        >  1:1648366/1‑150 (MQ=255)
                                                                                                                                                gCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCACTAACGGGGCCGTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGCTCCCAAGTCCATTCGTTTc       >  1:1279910/1‑150 (MQ=255)
                                                                                                                                                gCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCTCATCGACCTTAGCAAATCGTCGTGCTTCGCTCCCAAGTCCATTCGTTTc       >  2:275612/1‑150 (MQ=255)
                                                                                                                                                     g‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGCTCCCAAGTCCATTCGTTTCCCCgt  <  2:418548/150‑1 (MQ=255)
                                                                                                                                                     g‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGCTCCCAAGTCCATTCGTTTCCCCgt  >  1:1117687/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTTGAATAGGGCGAAGGGCCGGCCAACGCCTGG‑CTGCAGGCCAGGGCCAGTGAAGCGAGGGTAAAGCGCAATAACGGGGCTTTTCGCATCGGGGAACATCCTTTTCTTTGTTTTTATCAGGGCAAGCCCATCGACCTTAGCAAATCGTCGTGCTTCGCTCCCAAGTCCATTCGTTTCCCCGT  >  NC_002947/508808‑509106

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