Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,026,541 T→C 5.1% intergenic (+16/+16) helD → / ← mgsA DNA helicase IV/methylglyoxal synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,026,5410TC 5.1% 360.7 / 6.8 118intergenic (+16/+16)helD/mgsADNA helicase IV/methylglyoxal synthase
Reads supporting (aligned to +/- strand):  ref base T (64/48);  new base C (2/4);  total (66/52)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.66e-01

Only 100 of 121 total aligned reads displayed.

AAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCG  >  NC_000913/1026392‑1026690
                                                                                                                                                     |                                                                                                                                                     
  gaaCGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGGTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGc                                                                                                                                                     <  2:1084570/150‑1 (MQ=255)
    aCGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCa                                                                                                                                                   >  2:1982022/1‑150 (MQ=255)
    aCGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTAGGGGCCCTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCa                                                                                                                                                   <  2:1429525/150‑1 (MQ=255)
      ggcggTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCAcc                                                                                                                                                 <  2:1558711/150‑1 (MQ=255)
         ggTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGCATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCCCCtgt                                                                                                                                              >  2:1355753/1‑150 (MQ=255)
         ggTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAATTACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCAAACGCCACCtgt                                                                                                                                              <  2:1395069/150‑1 (MQ=255)
          gTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCtgtg                                                                                                                                             <  2:522240/150‑1 (MQ=255)
           ttAATGTACGTGGCGCTGACCCTGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGc                                                                                                                                            <  2:1151249/150‑1 (MQ=255)
             aaTGTACGTGGCTCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAAAACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCAAACGCCACCtgtttct                                                                                                                                          <  1:447241/150‑5 (MQ=255)
              aTGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCAAACGCCACCtgtttctt                                                                                                                                         <  1:1467502/150‑6 (MQ=255)
                gTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCAAACGCCACCtgtttcttac                                                                                                                                       <  2:509094/150‑8 (MQ=255)
                 tatgTGGCGATGACCCGGTCACGCCAGCGGGTATTGGCACTGTTTAACAAAGAGTATCCCTCTCCCTTTGTGGATATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATAtt                                                                                                                                      <  2:1782493/147‑1 (MQ=255)
                 tACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATAtt                                                                                                                                      <  2:1480422/150‑1 (MQ=255)
                 tACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTATCCCTTTGTGGAATTACTGAAAAATCTTGATGTGCCGGTGGCGAGAAAACCGTAAGAAACATGTGGAGTTTGCCACCTGTGCAATAtt                                                                                                                                      >  1:872470/1‑150 (MQ=255)
                    gTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACt                                                                                                                                   <  2:1776394/150‑1 (MQ=255)
                     tGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACtt                                                                                                                                  <  1:659987/150‑1 (MQ=255)
                       gcgcTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCa                                                                                                                                <  1:493471/150‑1 (MQ=255)
                            gACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACgg                                                                                                                           >  2:313634/1‑150 (MQ=255)
                                 gggCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCcgc                                                                                                                      >  1:1472259/1‑150 (MQ=255)
                                   gCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCga                                                                                                                    <  2:1392813/150‑1 (MQ=255)
                                    cACGCCATCGGGTATGGGCACTGTTTAACAAAGAGTATGCCTCTCCCTTTGTGGAAATTCTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCgag                                                                                                                   <  1:157260/150‑1 (MQ=255)
                                    cACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAGAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCgag                                                                                                                   >  2:114397/1‑150 (MQ=255)
                                    cACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCgag                                                                                                                   <  2:1441268/150‑1 (MQ=255)
                                      cGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGAt                                                                                                                 >  1:1630486/1‑150 (MQ=255)
                                       gCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATa                                                                                                                >  1:640743/1‑150 (MQ=255)
                                        ccATCGGGTATGGACACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATaa                                                                                                               >  1:1452447/1‑150 (MQ=255)
                                         cATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAAc                                                                                                              >  1:1547496/1‑150 (MQ=255)
                                            cGGGTATGGGCACTGTTTAACAAAGAGACTCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGACCGCGAGATAACGCt                                                                                                           <  1:1853990/150‑1 (MQ=255)
                                                    ggCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCg                                                                                                   <  1:1147683/150‑1 (MQ=255)
                                                      cACTGTTTAACATAGACACTCCGTCTCCCTTTGTGGAAAAACTGAAAAATCTGGATGTGCCGGTGGCGCGAAAACCGTAAGAAACCGGGGGCGTTTGCCACCGGTGCAATATTACTTCAGCCGGTCCGCGAGATTACGCTGATAATgggg                                                                                                 <  2:104949/150‑1 (MQ=255)
                                                       acggTTTAACAAAGAGAATCCCTCTCCCTTTTTGGAAATACTGAATAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGGCCGCGAGATAACGCTGATAATCgggg                                                                                                <  1:136805/147‑1 (MQ=255)
                                                        cTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGGTAACGCTGATAATCGGGGa                                                                                               <  1:758152/150‑1 (MQ=255)
                                                        cTGTTTAACAAAGAGAATCCCTCTACCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGa                                                                                               <  1:371618/150‑1 (MQ=255)
                                                         tGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGAt                                                                                              <  1:600620/150‑1 (MQ=255)
                                                              aaCAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGaa                                                                                         <  1:1812953/150‑1 (MQ=255)
                                                              aaCAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGaa                                                                                         <  2:2032141/150‑1 (MQ=255)
                                                                cAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAata                                                                                       >  1:970931/1‑150 (MQ=255)
                                                                cAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGATGCTAATCGGGGATCAGAata                                                                                       >  2:1793444/1‑150 (MQ=255)
                                                                 aaaGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAatat                                                                                      >  2:501754/1‑150 (MQ=255)
                                                                 aaaGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAatat                                                                                      >  1:1816902/1‑150 (MQ=255)
                                                                  aaGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATc                                                                                     >  1:1881872/1‑150 (MQ=255)
                                                                      gaATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGAcc                                                                                 >  2:1370500/1‑150 (MQ=255)
                                                                         tCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTGAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACagcg                                                                              >  1:1051304/1‑150 (MQ=255)
                                                                         acccTCTCCCTTTGTGGAAATAATGAAAAATCTGGCTGTGCCGGTGGCGAGAAAACCGGAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACcgcg                                                                              <  2:402821/149‑1 (MQ=255)
                                                                          cccTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGAGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGcgt                                                                             >  1:150772/1‑150 (MQ=255)
                                                                            ctctCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGcgtcg                                                                           >  2:1345177/1‑150 (MQ=255)
                                                                               tcCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTg                                                                        >  1:1100629/1‑150 (MQ=255)
                                                                               tcCCTTTGTGCAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAACACCGTAAGAAACCGGTGGCGGTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGTTAACGCTGATAATCGGGGATCAGAATATCGACCGCGGCGTTg                                                                        <  1:1026603/150‑1 (MQ=255)
                                                                                cccTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATATCGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGa                                                                       >  1:1470654/1‑150 (MQ=255)
                                                                                cccTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGa                                                                       >  2:1835186/1‑150 (MQ=255)
                                                                                cccTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGa                                                                       >  2:683577/1‑150 (MQ=255)
                                                                                cccTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCGGATAATCGGGGATCAGAATATCGACCGCGTCGTTGa                                                                       >  2:1721780/1‑150 (MQ=255)
                                                                                 ccTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGaa                                                                      >  1:495859/1‑150 (MQ=255)
                                                                                       tgGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGTATCAGAATATCGACCGCGTCGTTGAAATgcgg                                                                >  2:1120073/1‑150 (MQ=255)
                                                                                       tgGAAATACTCAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGAGCAATATTACTACAGACGGTCCGCGAGATAACGCTGATAAACGGGGATCAGAATATCGACCGCGTCGTAGAAATccgg                                                                >  2:500092/1‑150 (MQ=255)
                                                                                        ggAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATgcggc                                                               >  2:682302/1‑150 (MQ=255)
                                                                                              aCTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTgg                                                         >  2:2019729/1‑150 (MQ=255)
                                                                                              aCTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTgg                                                         <  2:1855273/150‑1 (MQ=255)
                                                                                               cTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGGAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGa                                                        >  2:653332/1‑150 (MQ=255)
                                                                                                tGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGAt                                                       >  2:1437543/1‑150 (MQ=255)
                                                                                                tGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGAt                                                       <  2:694596/150‑1 (MQ=255)
                                                                                                tGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGAt                                                       <  2:1737743/150‑1 (MQ=255)
                                                                                                tGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGAt                                                       >  2:1889309/1‑150 (MQ=255)
                                                                                                 gAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGAtt                                                      >  1:1603518/1‑150 (MQ=255)
                                                                                                 gAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGAtt                                                      >  1:1365104/1‑150 (MQ=255)
                                                                                                 gAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGAtt                                                      <  2:1631011/150‑1 (MQ=255)
                                                                                                 gAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGAtt                                                      >  1:1980449/1‑150 (MQ=255)
                                                                                                 gAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAAAGCGGCGACTGGAtt                                                      >  2:48971/1‑150 (MQ=255)
                                                                                                   aaaaTCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTAt                                                    >  1:1084772/1‑150 (MQ=255)
                                                                                                    aaaTCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATg                                                   <  1:673615/150‑1 (MQ=255)
                                                                                                     aaTCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGa                                                  <  2:754580/150‑1 (MQ=255)
                                                                                                         tGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGCCGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGCCCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTc                                              >  2:333065/1‑150 (MQ=255)
                                                                                                             tgtgCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGcc                                          >  2:732403/1‑150 (MQ=255)
                                                                                                             tgtgCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGcc                                          <  2:1881871/150‑1 (MQ=255)
                                                                                                              gtgCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCg                                         <  1:210492/150‑1 (MQ=255)
                                                                                                               tgCCGGTGGCGAGAAAACCGTAAGAAACCGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGGGTCGTTGAAATGCGTCGACTGGATTATGAAGTCTGCAGt                                        >  1:57504/1‑150 (MQ=255)
                                                                                                                gCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGtt                                       >  2:944126/1‑150 (MQ=255)
                                                                                                                gCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGtt                                       >  1:89480/1‑150 (MQ=255)
                                                                                                                      ggCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCCCg                                 >  2:1521698/1‑150 (MQ=255)
                                                                                                                                 cGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGAc                      >  2:950905/1‑150 (MQ=255)
                                                                                                                                  gTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGAcc                     >  1:474703/1‑150 (MQ=255)
                                                                                                                                  gTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGAcc                     >  1:1572555/1‑150 (MQ=255)
                                                                                                                                   tAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGCCCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCg                    >  2:270551/1‑150 (MQ=255)
                                                                                                                                   tAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCg                    <  1:1277200/150‑1 (MQ=255)
                                                                                                                                      gAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGaa                 <  1:831355/150‑1 (MQ=255)
                                                                                                                                      gAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGaa                 >  1:1850508/1‑150 (MQ=255)
                                                                                                                                        caCAGGTGGCAAACGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATg               >  2:648919/2‑150 (MQ=255)
                                                                                                                                         aCAGGTGGCAAACGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGt              >  1:329067/1‑150 (MQ=255)
                                                                                                                                           aGGTGTCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTc            <  1:566747/150‑1 (MQ=255)
                                                                                                                                           aGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTc            <  2:1628300/150‑1 (MQ=255)
                                                                                                                                            gggggCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGATCGGAATGTTcc           <  2:1922398/147‑1 (MQ=255)
                                                                                                                                            ggTGGCGTTTGCCACCTGTGCAATATTACTTCAGCCGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGTCCGGAATGTTcc           >  1:1817330/1‑150 (MQ=255)
                                                                                                                                            ggTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTcc           <  2:1822708/150‑1 (MQ=255)
                                                                                                                                             ggggCGATTGCCACCTGTGCAATATTACTTCAGACGGGCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCa          <  1:384404/148‑1 (MQ=255)
                                                                                                                                               ggCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATa        >  2:2015853/1‑150 (MQ=255)
                                                                                                                                               ggCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTACCACGTTGGTGGCGACCGGAATGTTCCATa        <  1:313634/150‑1 (MQ=255)
                                                                                                                                                 cGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGAGTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGCCCGGAATGTTCCATAcc      >  2:150978/1‑150 (MQ=255)
                                                                                                                                                  gTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCg     <  1:1364418/150‑1 (MQ=255)
                                                                                                                                                   tttGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGTGACCGGAATGTTCCATACCGt    <  1:1785666/150‑1 (MQ=255)
                                                                                                                                                     tGCCACCAGTGCAATATTACTTCAGACGGTCCTCGAGATAACGCTGACAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCg  <  1:1280276/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCG  >  NC_000913/1026392‑1026690

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