Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,388,379 C→G 7.5% intergenic (+73/+71) pyrD → / ← rmf dihydroorotate dehydrogenase/ribosome modulation factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,388,3790CG 7.5% 340.3 / 15.0 107intergenic (+73/+71)pyrD/rmfdihydroorotate dehydrogenase/ribosome modulation factor
Reads supporting (aligned to +/- strand):  ref base C (49/50);  new base G (5/3);  total (54/53)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

Only 100 of 108 total aligned reads displayed.

TATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACCGCGAC  >  NC_002947/2388230‑2388526
                                                                                                                                                     |                                                                                                                                                   
tataCTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTTCAGCTGATGAGGTTTGCAAAGTCATAAAAAAGGGCCCCTTc                                                                                                                                                     <  2:217161/150‑1 (MQ=255)
    cTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAgg                                                                                                                                                 <  2:2226962/150‑1 (MQ=255)
     tCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAggg                                                                                                                                                >  1:422528/1‑150 (MQ=255)
     tCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAggg                                                                                                                                                >  1:1359493/1‑150 (MQ=255)
       gggCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGc                                                                                                                                              >  1:982175/1‑150 (MQ=255)
         gCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTGAAGGGGccc                                                                                                                                            <  1:1680504/150‑1 (MQ=255)
          cTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAg                                                                                                                                                  >  1:1624388/1‑143 (MQ=255)
          cTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGcccc                                                                                                                                           >  2:1666338/1‑150 (MQ=255)
             cATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGCTGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTgg                                                                                                                                        <  2:2145208/150‑1 (MQ=255)
             cATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTgg                                                                                                                                        <  1:162114/150‑1 (MQ=255)
               tCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGc                                                                                                                                      >  1:2152107/1‑150 (MQ=255)
                cTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGcc                                                                                                                                     <  2:2237312/150‑1 (MQ=255)
                 tACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCt                                                                                                                                    >  1:51958/1‑150 (MQ=255)
                 tACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCt                                                                                                                                    >  1:1922842/1‑150 (MQ=255)
                 tACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCACCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCt                                                                                                                                    <  2:1799697/150‑1 (MQ=255)
                    aaaGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTcg                                                                                                                                 >  2:1073357/1‑150 (MQ=255)
                     aaGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTcgc                                                                                                                                <  1:1311353/150‑1 (MQ=255)
                     aaGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTcgc                                                                                                                                <  2:1140955/150‑1 (MQ=255)
                              cgcTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAATGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGccc                                                                                                                       <  1:2290997/150‑1 (MQ=255)
                               gcTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCg                                                                                                                      >  1:2126103/1‑150 (MQ=255)
                                  gATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGAt                                                                                                                   <  2:2368780/150‑1 (MQ=255)
                                     tCGCGAAGCGGTGGATGCTATCGCGGCAAGCCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATggg                                                                                                                <  1:506855/150‑1 (MQ=255)
                                       gcgAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGc                                                                                                              >  1:277001/1‑150 (MQ=255)
                                          aaGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTg                                                                                                           <  2:1922842/150‑1 (MQ=255)
                                            gCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCa                                                                                                         <  2:935101/150‑1 (MQ=255)
                                              ggTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATg                                                                                                       <  2:715233/150‑1 (MQ=255)
                                               gTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATgg                                                                                                      >  1:1628375/1‑150 (MQ=255)
                                                   aTGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTaaa                                                                                                  >  1:2129621/1‑150 (MQ=255)
                                                     gtcaccgcgGCATTGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTgg                                                                                                                                        >  2:1157301/6‑110 (MQ=255)
                                                     gtcaccgcgGCATTGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTgg                                                                                                                                        <  1:1157301/105‑1 (MQ=255)
                                                      tcaccgcgGCATTGCCGCGATGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTaaa                                                                                                  <  1:406461/143‑1 (MQ=255)
                                                       tATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGa                                                                                              >  1:1311625/1‑150 (MQ=255)
                                                        accgcgGCATTGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGaa                                                                                             >  1:1878107/3‑150 (MQ=255)
                                                        accgcgGCATTGCCGCGATGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATg                                                                                                       <  2:1956792/138‑1 (MQ=255)
                                                        accgcgGCATTGCCGCGATGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATg                                                                                                       >  1:1956792/3‑140 (MQ=255)
                                                        aTCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGaa                                                                                             <  1:253891/150‑1 (MQ=255)
                                                         tCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAAt                                                                                            >  2:5773/1‑150 (MQ=255)
                                                         tCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAAt                                                                                            >  2:28010/1‑150 (MQ=255)
                                                            cgGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCa                                                                                         <  1:2323285/150‑1 (MQ=255)
                                                            cgGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCa                                                                                         <  2:2152107/150‑1 (MQ=255)
                                                            cgGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCa                                                                                         <  1:28010/150‑1 (MQ=255)
                                                             ggCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAg                                                                                        >  2:1423286/1‑150 (MQ=255)
                                                                aaTGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTa                                                                                                    <  2:588063/135‑1 (MQ=255)
                                                                aaTGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTa                                                                                                    >  1:588063/1‑135 (MQ=255)
                                                                  tGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTcc                                                                                   <  1:1073357/150‑1 (MQ=255)
                                                                  tGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTcc                                                                                   <  1:5773/150‑1 (MQ=255)
                                                                   gCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGACATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCt                                                                                  >  1:483861/1‑150 (MQ=255)
                                                                      gcgGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGAt                                                                                                                   >  2:332519/1‑114 (MQ=255)
                                                                           gACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGc                                                                          <  2:1414030/150‑1 (MQ=255)
                                                                                  ccgcccCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGc                                                                   >  2:1594426/1‑150 (MQ=255)
                                                                                   cgcccCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCg                                                                  <  2:1628375/150‑1 (MQ=255)
                                                                                    gcccCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGt                                                                 >  2:230859/1‑150 (MQ=255)
                                                                                    gcccCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGACCCTGGGCCTTCGCCCGGGGCGCTGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGt                                                                 <  1:1499761/150‑1 (MQ=255)
                                                                                    cccacTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCACCGCAATAGAGAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGt                                                                 <  1:149305/146‑1 (MQ=255)
                                                                                     ccccTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTg                                                                <  1:1594426/150‑1 (MQ=255)
                                                                                      cccTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGccctt                                                                                                                                          >  1:1805731/1‑73 (MQ=255)
                                                                                      cccTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGccctt                                                                                                                                          <  2:1805731/75‑3 (MQ=255)
                                                                                               ggggCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTg                                                      >  1:90330/1‑150 (MQ=255)
                                                                                               gggaCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGCGCCCCTTCCAGGGGCCACTGTGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTCATTCGAATTCAGTGTCCTCGTTTATCCAACTGAGTGATTGTCGttc                                                      >  2:433053/1‑149 (MQ=255)
                                                                                                gggCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGa                                                     <  2:51958/150‑1 (MQ=255)
                                                                                                          gtgCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGt                                                                                                     <  1:2008124/92‑1 (MQ=255)
                                                                                                          gtgCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGt                                                                                                     >  2:2008124/1‑92 (MQ=255)
                                                                                                          gtgCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGAt                                           <  2:580280/150‑1 (MQ=255)
                                                                                                            gCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGc                                         <  2:483861/150‑1 (MQ=255)
                                                                                                              aGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCa                                       <  1:1423286/150‑1 (MQ=255)
                                                                                                              aGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGAGTGATTTTCGTTGAGTCTATGGATGCCa                                       <  1:1874444/150‑1 (MQ=255)
                                                                                                                  gatgaGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGAGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCgg                                   >  1:2105083/1‑150 (MQ=255)
                                                                                                                      aGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGcc                               >  1:556342/1‑150 (MQ=255)
                                                                                                                       ggTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGTATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCa                              <  2:918179/150‑1 (MQ=255)
                                                                                                                       ggTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCa                              <  2:277001/150‑1 (MQ=255)
                                                                                                                         tttGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTc                                                          <  2:749286/120‑1 (MQ=255)
                                                                                                                         tttGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTc                                                          >  1:749286/1‑120 (MQ=255)
                                                                                                                             ccAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCAtt                        >  1:1291039/1‑150 (MQ=255)
                                                                                                                             ccAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCAtt                        >  1:1457852/1‑150 (MQ=255)
                                                                                                                               aaGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTcc                      >  2:1795730/1‑150 (MQ=255)
                                                                                                                               aaGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGCGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCCGTCATTcc                      >  2:2299994/1‑150 (MQ=255)
                                                                                                                                    aTAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTccc                 >  1:1720966/1‑150 (MQ=255)
                                                                                                                                     tAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCa                >  1:2245632/1‑150 (MQ=255)
                                                                                                                                      aaaaaaGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAg               <  2:1878107/150‑1 (MQ=255)
                                                                                                                                       cacacgggCACCCTCATGGGGCCTCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGt              <  2:1720966/145‑1 (MQ=255)
                                                                                                                                        ccaggggCCCCTTGAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGtt             >  1:802590/5‑150 (MQ=255)
                                                                                                                                        ccaggggCCCCTTGAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGtt             >  2:733709/5‑150 (MQ=255)
                                                                                                                                         aaaGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGGATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTAt                                               >  1:1502543/1‑115 (MQ=255)
                                                                                                                                         aaaGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTAt                                               <  2:1502543/115‑1 (MQ=255)
                                                                                                                                            gggCCCCTTGAAGGGGCCCTTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTAt                                               >  1:1192409/1‑112 (MQ=255)
                                                                                                                                            gggCCCCTTGAAGGGGCCCTTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTGTCGTTGAGTCTAt                                               <  2:1192409/112‑1 (MQ=255)
                                                                                                                                            gggCCCCTTGAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGt                                                                 <  2:1629192/94‑1 (MQ=255)
                                                                                                                                            gggCCCCTTGAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGt                                                                 >  1:1629192/1‑94 (MQ=255)
                                                                                                                                            gggCCCCTTGAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCa         >  1:407208/1‑150 (MQ=255)
                                                                                                                                              gCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGt                                                                                                     <  2:2345181/56‑1 (MQ=255)
                                                                                                                                              gCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGt                                                                                                     >  1:2345181/1‑56 (MQ=255)
                                                                                                                                               ccccTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCa                                                                                         >  2:436692/1‑67 (MQ=255)
                                                                                                                                               ccccTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCa                                                                                         <  1:436692/67‑1 (MQ=255)
                                                                                                                                                cccTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACcgc     >  1:638367/1‑150 (MQ=255)
                                                                                                                                                cccTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACcgc     >  2:1849065/1‑150 (MQ=255)
                                                                                                                                                 ccTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACcgcg    >  2:1621660/1‑150 (MQ=255)
                                                                                                                                                  cTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCAc        >  1:258351/1‑145 (MQ=255)
                                                                                                                                                  cTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCAc        <  2:258351/145‑1 (MQ=255)
                                                                                                                                                   ttCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACCGCGAc  <  1:1310072/150‑1 (MQ=255)
                                                                                                                                                   ttCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGACGGGGCGTCACTGGTACGTGGAATTCAGTGTCCTATGGCTGCCAACGGCGTGATTTTCGTTGATTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACCGCGAc  >  1:45532/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACCGCGAC  >  NC_002947/2388230‑2388526

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: