Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,212,475 T→C 5.2% intergenic (‑171/+65) ymgK ← / ← ymgL protein YmgK/protein YmgL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,212,4750TC 5.2% 479.9 / 10.8 154intergenic (‑171/+65)ymgK/ymgLprotein YmgK/protein YmgL
Reads supporting (aligned to +/- strand):  ref base T (69/76);  new base C (4/4);  total (73/81)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

Only 100 of 160 total aligned reads displayed.

CGATAATAATGAATATAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATA  >  NC_000913/1212326‑1212622
                                                                                                                                                      |                                                                                                                                                   
cGATAATAATGAATATAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTNTNNCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACgtgt                                                                                                                                                     <  2:21896/150‑1 (MQ=255)
  ataataATGAATATAAGCATGACATGAGAATAAGGTTTCAATTTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGtgtg                                                                                                                                                    <  1:126401/150‑1 (MQ=255)
     ataatGAATATAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCAc                                                                                                                                                <  1:1218031/150‑1 (MQ=255)
      taatGAATATAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACg                                                                                                                                               <  1:653502/150‑1 (MQ=255)
        atGAATATAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGc                                                                                                                                             >  1:218959/1‑150 (MQ=255)
          taaTATAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCg                                                                                                                                           <  2:1091507/149‑1 (MQ=255)
            atatAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGCCGTGCCACACGTCCgtt                                                                                                                                         <  1:1307033/150‑1 (MQ=255)
             tataAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCgttg                                                                                                                                        <  1:1081549/150‑1 (MQ=255)
              ataAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGCCGTGCCACACGTCCgtttc                                                                                                                                       <  2:1756816/150‑3 (MQ=255)
               taAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGCCGTGCCACACGTCCgtttca                                                                                                                                      <  1:328627/150‑4 (MQ=255)
                aaGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGttt                                                                                                                                     <  2:1180168/150‑1 (MQ=255)
                aaGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGttt                                                                                                                                     <  2:124195/150‑1 (MQ=255)
                aaGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGttt                                                                                                                                     <  1:456018/150‑1 (MQ=255)
                 aGCTTGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGTAATGATTTTACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTTCGTTAAGTATTGAAACGGCCGTGGCACACGTCCgtttcaat                                                                                                                                    <  2:1414336/150‑6 (MQ=255)
                 aGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGtttt                                                                                                                                    <  1:1112820/150‑1 (MQ=255)
                  gCATGACATGACAATATGGTTTAAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTTTCTAATATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGttttt                                                                                                                                   <  2:1162332/150‑1 (MQ=255)
                   catgacatgaGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTta                                                                                                                                  >  1:156517/1‑150 (MQ=255)
                    atgacatgaGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTtat                                                                                                                                 <  2:1440283/150‑1 (MQ=255)
                     tgacatgaGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTtata                                                                                                                                <  2:882008/150‑1 (MQ=255)
                      cacatgaGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATaa                                                                                                                               >  1:270418/2‑150 (MQ=255)
                       acatgaGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATaaa                                                                                                                              <  2:1692560/150‑1 (MQ=255)
                        catgaGAATAAGGTTTCAATTTTTGATTTATATAGTAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAat                                                                                                                             <  2:1789862/150‑1 (MQ=255)
                        catgaGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAat                                                                                                                             <  2:1497724/150‑1 (MQ=255)
                        catgaGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAat                                                                                                                             >  2:1974874/1‑150 (MQ=255)
                         atgaGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAata                                                                                                                            <  2:1857087/150‑1 (MQ=255)
                          tgaGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAatat                                                                                                                           >  2:1845177/1‑150 (MQ=255)
                          tgaGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAatat                                                                                                                           <  2:977141/150‑1 (MQ=255)
                           tagaATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATg                                                                                                                          <  2:1746739/149‑1 (MQ=255)
                            agaATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGt                                                                                                                         >  2:1773760/1‑150 (MQ=255)
                               aTAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTaa                                                                                                                      <  2:1629852/150‑1 (MQ=255)
                                     tttCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTta                                                                                                                <  1:523839/150‑1 (MQ=255)
                                      ttCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTtat                                                                                                               >  2:1645016/1‑150 (MQ=255)
                                        cAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATaa                                                                                                             >  2:1190860/1‑150 (MQ=255)
                                        cAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATaa                                                                                                             >  1:1389239/1‑150 (MQ=255)
                                        cAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATaa                                                                                                             <  1:47289/150‑1 (MQ=255)
                                         aaTTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGCCCGTTGTTTTTATAAATATGTTAACCGTTATaaa                                                                                                            >  1:185972/1‑150 (MQ=255)
                                          aTTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATaaaa                                                                                                           <  2:919859/150‑1 (MQ=255)
                                            tttttGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAAt                                                                                                          <  2:1738819/150‑1 (MQ=255)
                                            tttttGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAAt                                                                                                          >  1:220751/1‑150 (MQ=255)
                                                tGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATTTGTTAACCGTTATAAAATAACg                                                                                                      <  2:1394474/150‑1 (MQ=255)
                                                 taGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGt                                                                                                     <  2:1439004/149‑1 (MQ=255)
                                                   gTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTAt                                                                                                   >  1:705912/1‑150 (MQ=255)
                                                   gTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTAt                                                                                                   >  1:292692/1‑150 (MQ=255)
                                                     tatataGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGAGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCa                                                                                                 >  1:1026762/1‑150 (MQ=255)
                                                       tataGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCaaa                                                                                               >  1:1050619/1‑150 (MQ=255)
                                                       tataGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACAAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCaaa                                                                                               >  1:696260/1‑150 (MQ=255)
                                                         taGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTTCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAg                                                                                             <  1:1897445/150‑1 (MQ=255)
                                                         taGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAg                                                                                             <  2:704067/150‑1 (MQ=255)
                                                          aGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGCCGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGt                                                                                            >  2:1583673/1‑150 (MQ=255)
                                                           ggAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTc                                                                                           <  2:1069245/150‑1 (MQ=255)
                                                           ggAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTc                                                                                           >  1:776990/1‑150 (MQ=255)
                                                             aaTGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCaa                                                                                         >  1:377411/1‑150 (MQ=255)
                                                             aaTGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCaa                                                                                         <  2:1389239/150‑1 (MQ=255)
                                                                 aTTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGa                                                                                     >  1:1457551/1‑150 (MQ=255)
                                                                     aaCCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATcac                                                                                 >  2:222843/1‑150 (MQ=255)
                                                                     aaCCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATcac                                                                                 >  1:852279/1‑150 (MQ=255)
                                                                     aaCCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATcac                                                                                 <  1:1976690/150‑1 (MQ=255)
                                                                     aaCCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATcac                                                                                 >  1:1929371/1‑150 (MQ=255)
                                                                          gTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTc                                                                            <  2:852279/150‑1 (MQ=255)
                                                                              cTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCaaaa                                                                                              >  1:150320/1‑128 (MQ=255)
                                                                               tGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAata                                                                       <  1:376380/150‑1 (MQ=255)
                                                                                ggCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAatat                                                                      >  1:817598/1‑150 (MQ=255)
                                                                                ggCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAatat                                                                      >  2:79517/1‑150 (MQ=255)
                                                                                  cTAAATTACATATAACCGGATTATGACTAAACCATAATACATTTGCTTTAAGTATTTAAACGGACGTGTGTCACGTCCGTTGTTTTTATAAATATTTTAACCTTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAatatga                                                                    <  1:909765/150‑3 (MQ=255)
                                                                                  cTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCa                                                                    <  2:610951/150‑1 (MQ=255)
                                                                                       aTACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCa                                                                    >  1:1486010/1‑145 (MQ=255)
                                                                                         aCATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGAt                                                             <  1:1639351/150‑1 (MQ=255)
                                                                                          cATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATa                                                            <  2:1050619/150‑1 (MQ=255)
                                                                                           atatAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAg                                                           <  1:938479/150‑1 (MQ=255)
                                                                                               aaCCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTATCCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTAtt                                                       <  2:270418/150‑1 (MQ=255)
                                                                                                 ccGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAg                                                     <  1:1760197/150‑1 (MQ=255)
                                                                                                       gatgaCTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGt                                               <  1:936205/150‑1 (MQ=255)
                                                                                                        atgaCTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTg                                              <  1:1773760/150‑1 (MQ=255)
                                                                                                            cTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGatta                                          >  2:474754/1‑150 (MQ=255)
                                                                                                              aaaCCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGattatt                                        <  2:1966675/150‑1 (MQ=255)
                                                                                                              aaaCCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGattatt                                        >  1:471530/1‑150 (MQ=255)
                                                                                                              aaaCCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGattatt                                        >  2:1762844/1‑150 (MQ=255)
                                                                                                                    aaaTACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATggg                                  <  2:1101281/150‑1 (MQ=255)
                                                                                                                    aaaTACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATggg                                  <  1:1439977/150‑1 (MQ=255)
                                                                                                                      aTACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTg                                <  2:1457551/150‑1 (MQ=255)
                                                                                                                      aTACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTg                                <  2:322664/150‑1 (MQ=255)
                                                                                                                       tACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGa                               >  2:192813/1‑150 (MQ=255)
                                                                                                                        aCATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGAc                              >  2:994787/1‑150 (MQ=255)
                                                                                                                        aCATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGAc                              >  2:891340/1‑150 (MQ=255)
                                                                                                                              gCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATaaa                        >  2:1375214/1‑150 (MQ=255)
                                                                                                                               ctttAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATaaaa                       <  1:1721686/148‑1 (MQ=255)
                                                                                                                                gTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATaaaaa                      >  1:16245/1‑150 (MQ=255)
                                                                                                                                 ttAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAg                     <  2:220751/150‑1 (MQ=255)
                                                                                                                                 ttAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAg                     <  2:156517/150‑1 (MQ=255)
                                                                                                                                         tGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTAtt             <  2:98029/150‑1 (MQ=255)
                                                                                                                                         aacaaCGGCCGTGCCACACGTCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTAtt             >  2:1354364/4‑150 (MQ=255)
                                                                                                                                         aacaaCGGCCGTGCCACACGTCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTAtt             >  2:162103/4‑150 (MQ=255)
                                                                                                                                         aacaaCGGCCGTGCCACACGTCCGTTGTTTTTAAAAATATGTTAACCGTTATAAAATAACGTAACAAAAGTCAAGTGATCAAATTTCAAATATCAAGTTGATAGTATTAGTATGGTGATTATTTATGGGTGACAATAAAAATACAGTAtt             >  2:549001/4‑150 (MQ=255)
                                                                                                                                          gAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTAtt                                                       >  2:705195/1‑107 (MQ=255)
                                                                                                                                          gAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTa            >  1:590479/1‑150 (MQ=255)
                                                                                                                                          gAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTa            >  2:1455574/1‑150 (MQ=255)
                                                                                                                                           aaaCGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTaa           >  1:436702/1‑150 (MQ=255)
                                                                                                                                                  cGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCa    <  1:192813/150‑1 (MQ=255)
                                                                                                                                                    tgtgGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATa  >  1:1305052/1‑150 (MQ=255)
                                                                                                                                                    tgtgGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATa  <  1:1190860/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
CGATAATAATGAATATAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATA  >  NC_000913/1212326‑1212622

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