Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,388,379 C→G 6.0% 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 6.0% 287.8 / 6.5 93intergenic (+73/+71)pyrD/rmfdihydroorotate dehydrogenase/ribosome modulation factor
Reads supporting (aligned to +/- strand):  ref base C (47/40);  new base G (3/3);  total (50/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACCGCGAC  >  NC_002947/2388234‑2388526
                                                                                                                                                  |                                                                                                                                                   
cTCGGGCTTCATCTACAAAGGGCCGGAGGACATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAgg                                                                                                                                                 >  1:62809/1‑150 (MQ=255)
 tCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAggg                                                                                                                                                >  1:1010140/1‑150 (MQ=255)
 tCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAggg                                                                                                                                                >  2:59098/1‑150 (MQ=255)
      cTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCAAAGGCATAAAAAAGGGCCCCTTGAAGGGGccct                                                                                                                                           <  2:2443143/150‑2 (MQ=255)
      cTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGACCTCCGCGGTGGGGGAGACGGCGGGAACCCAGGCGGAGCTTGGCCGAGGCATAAAAAAGGGCCCCTTCAAGGGGcccc                                                                                                                                           >  1:534384/1‑150 (MQ=255)
         cATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTgg                                                                                                                                        <  2:702480/150‑1 (MQ=255)
         cATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTgg                                                                                                                                        <  1:275149/150‑1 (MQ=255)
          aTCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTggg                                                                                                                                       >  2:1581226/1‑150 (MQ=255)
          aTCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTggg                                                                                                                                       <  1:59098/150‑1 (MQ=255)
           tCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGGCCCCTTGAAGGGGcccttttt                                                                                                                                       >  2:2093376/1‑145 (MQ=255)
           tCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGc                                                                                                                                      >  1:1328613/1‑150 (MQ=255)
             tACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGGCCCCTTGAAGGGGccctttttta                                                                                                                                     <  1:2093376/150‑8 (MQ=255)
               cAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTc                                                                                                                                  <  1:1561545/150‑1 (MQ=255)
               aaaaGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTc                                                                                                                                  >  1:1016415/2‑150 (MQ=255)
                aaaGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTgg                                                                                                                                        >  2:2377038/1‑143 (MQ=255)
                aaaGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTgg                                                                                                                                        <  1:2377038/143‑1 (MQ=255)
                              gATTCGCGAAGCTGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGATTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGAt                                                                                                                   <  2:1328613/150‑1 (MQ=255)
                               aTTCGCGAAGCGTTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATg                                                                                                                  <  2:215355/150‑1 (MQ=255)
                                   gcgAAGCGGGGGATGCTATCGCGGCAATGCCGCGGTGACCAGACCGCCCCTGTGGGGGGGCGGGGGGGTGCAGCTGATGAGGTTTGCGAAGGGATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGGCGCCCGGATGGGGc                                                                                                              >  2:1036795/1‑150 (MQ=255)
                                         cGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATgcggcgcgttt                                                                                                        <  2:1808582/150‑6 (MQ=255)
                                          ggTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATg                                                                                                       <  2:1992951/150‑1 (MQ=255)
                                          ggTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATg                                                                                                       >  1:1992951/1‑150 (MQ=255)
                                            tGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGt                                                                                                     >  1:1008507/1‑150 (MQ=255)
                                              gATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCCCACGCCGCCCGGATGGGGCGTGCATGGTaa                                                                                                   >  1:1567568/1‑150 (MQ=255)
                                               aTGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTaaa                                                                                                  >  1:2282897/1‑150 (MQ=255)
                                                tGCTATCGCGGCAATGCCGCGGTGAGCAGCGCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAAt                                                                                                 >  1:1775373/1‑150 (MQ=255)
                                                 gCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCg                                                                                               <  2:354178/150‑1 (MQ=255)
                                                   tATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGa                                                                                              >  2:475864/1‑150 (MQ=255)
                                                    accgcgGCATTGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTaa                                                                                                   >  2:72426/3‑144 (MQ=255)
                                                    accgcgGCATTGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTaa                                                                                                   <  1:72426/142‑1 (MQ=255)
                                                    aTCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCAt                                                                                                        >  2:136070/1‑139 (MQ=255)
                                                    aTCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCAt                                                                                                        <  1:136070/139‑1 (MQ=255)
                                                     tCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGcccc                                                                                                                                           >  1:1875657/1‑103 (MQ=255)
                                                     tCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGcccc                                                                                                                                           <  2:1875657/103‑1 (MQ=255)
                                                                ccGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTc                                                                                 <  1:37854/150‑1 (MQ=255)
                                                                ccGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTc                                                                                 <  1:1051875/150‑1 (MQ=255)
                                                                  gcgGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGaa                                                                                             >  2:1671950/1‑136 (MQ=255)
                                                                  gcgGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGaa                                                                                             <  1:1671950/136‑1 (MQ=255)
                                                                   cgGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGtt                                                                              <  2:2272677/150‑1 (MQ=255)
                                                                     gTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCa                                                                            <  2:1010140/150‑1 (MQ=255)
                                                                        aCCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGcc                                                                         >  1:921034/1‑150 (MQ=255)
                                                                        aCCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGcc                                                                         >  1:1604168/1‑150 (MQ=255)
                                                                          cAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCaa                                                                       >  2:972213/1‑150 (MQ=255)
                                                                            gcccgcccCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACg                                                                     <  1:1581226/150‑1 (MQ=255)
                                                                             cccgcccCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACgg                                                                    >  2:714773/1‑150 (MQ=255)
                                                                              ccgcccCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGc                                                                   <  2:2282897/150‑1 (MQ=255)
                                                                              ccgcccCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGc                                                                   >  1:263756/1‑150 (MQ=255)
                                                                                gcccCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGt                                                                 >  2:2389767/1‑150 (MQ=255)
                                                                                  cccTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGa                                                               >  1:1813671/1‑150 (MQ=255)
                                                                                          gggggCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGtt                                                       <  2:1775373/150‑1 (MQ=255)
                                                                                          gggggCGGGTGGGTGCAGCTCGTGAGGTTTGCCAAGGCATAAAAAAGGTCCCCTGCAAGGGGCCCCTGGGCCTTCTCCCGCCGCCCAGATTGGTCGTGCATGGTAGTTATTTTTTTGTTTCCTCTTTCAGCCAACGTCGTGATTTTCGtt                                                       <  1:958446/150‑1 (MQ=255)
                                                                                            gggCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGa                                                     <  1:1500634/150‑1 (MQ=255)
                                                                                            gggCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGa                                                     <  1:2389767/150‑1 (MQ=255)
                                                                                              gCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGt                                                   <  2:533713/150‑1 (MQ=255)
                                                                                                  gtgggtgCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTAt                                               >  1:1184774/1‑150 (MQ=255)
                                                                                                      gtgCAGCTGATGATGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGt                                                                 <  1:726255/128‑1 (MQ=255)
                                                                                                      gtgCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGt                                                                 >  2:726255/1‑128 (MQ=255)
                                                                                                      gtgCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGAt                                           >  1:1769149/1‑150 (MQ=255)
                                                                                                            cTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGc                                     >  1:1091469/1‑150 (MQ=255)
                                                                                                             tgatgaGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGaa                                                                                             <  2:804747/93‑1 (MQ=255)
                                                                                                             tgatgaGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGaa                                                                                             >  1:804747/1‑93 (MQ=255)
                                                                                                              gatgaGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCgg                                   >  1:469143/1‑150 (MQ=255)
                                                                                                               atgaGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGt                                  >  2:1271263/1‑150 (MQ=255)
                                                                                                               atgaGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGt                                  <  2:1091469/150‑1 (MQ=255)
                                                                                                                 gaGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATTGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGc                                <  1:710365/150‑1 (MQ=255)
                                                                                                                    gTTTGCCAAGGCATAAGAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCa                                                                                                         >  2:1739660/1‑74 (MQ=255)
                                                                                                                    gTTTGCCAAGGCATAAGAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCa                                                                                                         <  1:1739660/74‑1 (MQ=255)
                                                                                                                    gTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAg                             <  1:1271263/150‑1 (MQ=255)
                                                                                                                     tttGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCa                              >  1:2317121/1‑148 (MQ=255)
                                                                                                                     tttGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCa                              <  2:2317121/148‑1 (MQ=255)
                                                                                                                         ccAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCAtt                        >  1:892879/1‑150 (MQ=255)
                                                                                                                          cAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGtt                                                                              >  1:2249020/1‑95 (MQ=255)
                                                                                                                          cAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGtt                                                                              <  2:2249020/95‑1 (MQ=255)
                                                                                                                               cATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGACTTCGCCCGCCGCCCGGATGGGGCGTGCGTGGTAAATCGAATTCAGTGTCCTCGTTCAGCCGACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTACCAGTCATTCCGTcc                  >  2:1018736/1‑150 (MQ=255)
                                                                                                                                aTAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGt                    <  2:1673669/147‑1 (MQ=255)
                                                                                                                                aTAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGt                    >  1:1673669/1‑147 (MQ=255)
                                                                                                                                aTAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGAGTCCTCGTTCCGCAAACGGCGTGGTTTTCGTTGAGTCTATGGATGCCAGCGGCGCCAGTAATTCCGTccc                 >  1:2147138/1‑150 (MQ=255)
                                                                                                                                 tAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCCGCGGTGCCAGTCATTCCGTCCCa                <  2:1184774/150‑1 (MQ=255)
                                                                                                                                 tAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCa                >  2:2220368/1‑150 (MQ=255)
                                                                                                                                 tAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCa                <  2:1604168/150‑1 (MQ=255)
                                                                                                                                  aaaaaaGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGACACCATCGATGCCAGCGGTGCCAGTCATTCCGTCCCAg               >  2:2017300/1‑150 (MQ=255)
                                                                                                                                     ccaggggCCCCTTGAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGAt                                           <  1:989957/116‑1 (MQ=255)
                                                                                                                                     ccaggggCCCCTTGAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGAt                                           >  2:989957/5‑120 (MQ=255)
                                                                                                                                      caggggCCCCTTGAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTa            >  1:730511/4‑150 (MQ=255)
                                                                                                                                             cccTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGt                                                        <  2:2006307/99‑1 (MQ=255)
                                                                                                                                             cccTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGt                                                        >  1:2006307/1‑99 (MQ=255)
                                                                                                                                             cccTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACcgc     <  1:972213/150‑1 (MQ=255)
                                                                                                                                             cccTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACcgc     >  1:1012725/1‑150 (MQ=255)
                                                                                                                                             cccTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACcgc     <  1:714773/150‑1 (MQ=255)
                                                                                                                                             cccTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCCCGCACAGCCAACGGCGTGATTTTAGTTGAATATTGGCATTCCTGCTGTGCCAGTCATTCCGTCCCAGTTATCACcgc     >  1:1115971/1‑150 (MQ=255)
                                                                                                                                               cTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACagcaa   >  2:1450312/1‑148 (MQ=255)
                                                                                                                                                ttCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGTAGAAGATATCTAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACCGCGAc  <  2:62809/150‑1 (MQ=255)
                                                                                                                                                ttCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACCGCGAc  >  2:1140198/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACCGCGAC  >  NC_002947/2388234‑2388526

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: