Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,388,379 C→G 9.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 9.0% 271.3 / 15.5 89intergenic (+73/+71)pyrD/rmfdihydroorotate dehydrogenase/ribosome modulation factor
Reads supporting (aligned to +/- strand):  ref base C (41/40);  new base G (5/3);  total (46/43)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

ACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACCGC  >  NC_002947/2388233‑2388523
                                                                                                                                                   |                                                                                                                                                
aCTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAg                                                                                                                                               <  1:123448/150‑1 (MQ=255)
aCTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAg                                                                                                                                               <  1:1378806/150‑1 (MQ=255)
 cTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTGAAga                                                                                                                                              <  1:357651/150‑2 (MQ=255)
  tCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAggg                                                                                                                                             >  1:703785/1‑150 (MQ=255)
    gggCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGc                                                                                                                                           >  2:1628370/1‑150 (MQ=255)
       cTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGcccc                                                                                                                                        >  1:1809162/1‑150 (MQ=255)
        ttCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCACTGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCt                                                                                                                                       >  2:1246783/1‑150 (MQ=255)
          caccTACAAAGGGCCGGCGCTGATTCGCGCAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTgg                                                                                                                                     <  2:339982/147‑1 (MQ=255)
              tACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGGCCCCTTGAAGGGGccctttttta                                                                                                                                  >  2:1471880/1‑143 (MQ=255)
              tACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCt                                                                                                                                 >  2:998676/1‑150 (MQ=255)
              tACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCt                                                                                                                                 >  2:223500/1‑150 (MQ=255)
               aCAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCtt                                                                                                                                >  1:288584/1‑150 (MQ=255)
                   aGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTGAAGGGGccctttt                                                                                                                                     >  1:1137362/1‑137 (MQ=255)
                   aGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTGAAGGGGccctttt                                                                                                                                     <  2:1137362/141‑5 (MQ=255)
                      gCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTcgcccgc                                                                                                                         <  1:998676/150‑1 (MQ=255)
                        cGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCccgccg                                                                                                                       >  2:1882783/1‑150 (MQ=255)
                           cgcTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGccc                                                                                                                    <  1:1420616/150‑1 (MQ=255)
                                      gAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGt                                                                                                         >  1:1437436/1‑150 (MQ=255)
                                            gTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTGAAGGGGccctttt                                                                                                                                     >  1:1186074/1‑112 (MQ=255)
                                            gTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTGAAGGGGccctttt                                                                                                                                     <  2:1186074/116‑5 (MQ=255)
                                             tGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTCGGCCTTCGCCCGCCGCCCGGATGGGGCGGCCATGGt                                                                                                  >  2:669225/1‑150 (MQ=255)
                                                aTGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTaaa                                                                                               <  2:703785/150‑1 (MQ=255)
                                                 tgatATTGCGGCAAAGCCGCGGTGAGCAGCCCGCCCCGGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAAt                                                                                              >  2:1529337/4‑150 (MQ=255)
                                                 tGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATgggg                                                                                                            <  1:1598099/136‑1 (MQ=255)
                                                 tGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATgggg                                                                                                            >  2:1598099/1‑136 (MQ=255)
                                                 tGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAAt                                                                                              >  1:1735282/1‑150 (MQ=255)
                                                 tGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAAt                                                                                              <  2:1749277/150‑1 (MQ=255)
                                                    caccgcgGCATTGCCGCGATGACCAGCCCGCCCCTGTGGGGGGGCGGGGGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGACGTCTCCCGCCGCCCGGATGGGGCGTGCATGGGAAAGCGa                                                                                           >  1:1372322/4‑150 (MQ=255)
                                                     aTCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGaa                                                                                          <  2:1735282/150‑1 (MQ=255)
                                                     aTCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGaa                                                                                          <  2:885084/150‑1 (MQ=255)
                                                      tCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAAt                                                                                         >  1:497699/1‑150 (MQ=255)
                                                      tCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAAt                                                                                         <  1:1426290/150‑1 (MQ=255)
                                                       cgcgGCATTGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAAtt                                                                                        >  1:591431/1‑150 (MQ=255)
                                                        gcgGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTc                                                                                       >  1:582245/1‑150 (MQ=255)
                                                         cgGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCa                                                                                      >  2:33064/1‑150 (MQ=255)
                                                           gCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAgt                                                                                    <  1:967123/150‑1 (MQ=255)
                                                             aaTGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAgtgt                                                                                  >  2:121108/1‑150 (MQ=255)
                                                              aTGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTc                                                                                 <  1:1246783/150‑1 (MQ=255)
                                                                             gcccgcccCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGt                                                                                                         >  2:1778641/1‑111 (MQ=255)
                                                                             gcccgcccCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGt                                                                                                         <  1:1778641/111‑1 (MQ=255)
                                                                              cccgcccCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGCAGAGCAGAAGAGAGTGAATTCAGTGTCCTCGTTCAGCCAACgg                                                                 >  2:646455/1‑150 (MQ=255)
                                                                               ccgcccCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGc                                                                >  1:60726/1‑150 (MQ=255)
                                                                                   cccTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGaa                                                                                          >  1:143019/1‑120 (MQ=255)
                                                                                   cccTGGGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGaa                                                                                          <  2:143019/120‑1 (MQ=255)
                                                                                      tgtgGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGAttt                                                         >  1:937779/1‑150 (MQ=255)
                                                                                      tgtgGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGAttt                                                         <  2:1831650/150‑1 (MQ=255)
                                                                                        tgGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTc                                                       <  1:121108/150‑1 (MQ=255)
                                                                                         gggggggCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCg                                                      <  1:1628370/150‑1 (MQ=255)
                                                                                           aggggCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGAATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGtt                                                    <  1:101112/149‑1 (MQ=255)
                                                                                               gCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGt                                                <  1:496411/150‑1 (MQ=255)
                                                                                               gCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGt                                                >  2:496411/1‑150 (MQ=255)
                                                                                               gCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGt                                                <  1:1882783/150‑1 (MQ=255)
                                                                                                  ggtgggtgCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTa                                             >  2:1710413/1‑150 (MQ=255)
                                                                                                   gtgggtgCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTAt                                            <  1:33064/150‑1 (MQ=255)
                                                                                                       gtgCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGAt                                        <  2:1661805/150‑1 (MQ=255)
                                                                                                        tgCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATg                                       <  2:1809162/150‑1 (MQ=255)
                                                                                                        tgCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATg                                       >  1:1430332/1‑150 (MQ=255)
                                                                                                           aGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCa                                    >  1:1493800/1‑150 (MQ=255)
                                                                                                              tgatgaGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCg                                 <  1:223500/150‑1 (MQ=255)
                                                                                                                atgaGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCATGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGt                               >  2:609471/1‑150 (MQ=255)
                                                                                                                    ggTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGaa                                                                                          <  1:1843334/87‑1 (MQ=255)
                                                                                                                    ggTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGaa                                                                                          >  2:1843334/1‑87 (MQ=255)
                                                                                                                     gTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTc                                                                              >  2:691787/1‑98 (MQ=255)
                                                                                                                     gTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTc                                                                              <  1:691787/98‑1 (MQ=255)
                                                                                                                     gTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAg                          <  2:582245/150‑1 (MQ=255)
                                                                                                                     gTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAg                          <  2:497699/150‑1 (MQ=255)
                                                                                                                      tttGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGt                         <  2:1437436/150‑1 (MQ=255)
                                                                                                                        tGCCAAGGCATAAAAAAGGTCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCa                       <  2:1272836/150‑1 (MQ=255)
                                                                                                                         gCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGaa                                                                                          <  2:1854755/82‑1 (MQ=255)
                                                                                                                         gCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGaa                                                                                          >  1:1854755/1‑82 (MQ=255)
                                                                                                                            aaGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTcc                   >  2:1845167/1‑150 (MQ=255)
                                                                                                                                cATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTcc               <  2:1430332/150‑1 (MQ=255)
                                                                                                                                 aTAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAACTCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTTTATGGATGCCAGCGGTGCCAGTCATTCCGTccc              <  1:646455/150‑1 (MQ=255)
                                                                                                                                 aTAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGc                                                                >  2:552108/1‑100 (MQ=255)
                                                                                                                                 aTAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGc                                                                <  1:552108/100‑1 (MQ=255)
                                                                                                                                 aTAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCgg                                >  1:1318412/1‑132 (MQ=255)
                                                                                                                                 aTAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCgg                                <  2:1318412/132‑1 (MQ=255)
                                                                                                                                  tAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTAt                                            <  1:267208/119‑1 (MQ=255)
                                                                                                                                  tAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTAt                                            >  2:267208/1‑119 (MQ=255)
                                                                                                                                  tAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCa             <  2:591431/150‑1 (MQ=255)
                                                                                                                                  tAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCa             <  2:1493800/150‑1 (MQ=255)
                                                                                                                                     aaaaGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGtt          >  2:488425/1‑150 (MQ=255)
                                                                                                                                     aaaaGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGtt          >  2:679537/1‑150 (MQ=255)
                                                                                                                                        aggggCCCCTTGAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTAt        >  2:454401/3‑150 (MQ=255)
                                                                                                                                          gggCCCCTTGAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCa      >  2:234098/1‑150 (MQ=255)
                                                                                                                                             ccccTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACcg   >  1:579658/1‑150 (MQ=255)
                                                                                                                                              cccTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACcgc  >  2:1586958/1‑150 (MQ=255)
                                                                                                                                               ccTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACg                                                                  <  1:820610/85‑1 (MQ=255)
                                                                                                                                               ccTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACg                                                                  >  2:820610/1‑85 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
ACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACCGC  >  NC_002947/2388233‑2388523

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: