Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,845,789 A→G 6.6% M228T (ATG→ACG)  hycE ← formate hydrogenlyase subunit HycE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,845,7890AG 6.6% 412.1 / 14.0 137M228T (ATG→ACG) hycEformate hydrogenlyase subunit HycE
Reads supporting (aligned to +/- strand):  ref base A (71/57);  new base G (3/6);  total (74/63)
Fisher's exact test for biased strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.64e-01

Only 100 of 141 total aligned reads displayed.

CCTCCAGCAGAATGGCGCGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCACGACGTTG  >  NC_000913/2845640‑2845935
                                                                                                                                                     |                                                                                                                                                  
ccTCCAGCAGAATGGCGCGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCa                                                                                                                                                    <  2:2221978/150‑1 (MQ=255)
ccTCCAGCAGAATGGCGCGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAAACCa                                                                                                                                                    >  1:2455014/1‑150 (MQ=255)
 cTCCAGCAGAATGGCGCGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCAt                                                                                                                                                   <  2:157121/150‑1 (MQ=255)
      gcagAATGGCGCGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGACAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACggg                                                                                                                                              <  1:426048/150‑1 (MQ=255)
        agAATGGCGCGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGtt                                                                                                                                            <  1:36882/150‑1 (MQ=255)
         gAATGGCGCGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGttt                                                                                                                                           >  1:1511309/1‑150 (MQ=255)
           aTGGCGCGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCGTATGGGttata                                                                                                                                         <  1:918173/150‑4 (MQ=255)
            gtgcgcAGATGATCCGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGCTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCg                                                                                                                                        >  2:2004141/3‑150 (MQ=255)
             ggCGCGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCGTATGGGttataac                                                                                                                                       <  1:2246002/150‑6 (MQ=255)
             ggCGCGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCGTATGGGttataac                                                                                                                                       <  1:1140594/150‑6 (MQ=255)
             ggCGCGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCGTATGGGttataac                                                                                                                                       <  2:2148218/150‑6 (MQ=255)
             ggCGCGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGc                                                                                                                                       >  2:468923/1‑150 (MQ=255)
             ggCGCGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGc                                                                                                                                       >  2:2099105/1‑150 (MQ=255)
             ggCGCGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGc                                                                                                                                       >  2:1093500/1‑150 (MQ=255)
             ggCGCGGATCATCAGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGGCAGAGAGGAAGGTCACTTCGTTATAACCCGTATGGGttataac                                                                                                                                       <  1:1753257/150‑6 (MQ=255)
              gcgcGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCGTATGGGttataacg                                                                                                                                      <  1:1052564/150‑7 (MQ=255)
              gcgcGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGcc                                                                                                                                      <  1:1681589/150‑1 (MQ=255)
                  ggATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGtt                                                                                                                                  >  1:1447397/1‑150 (MQ=255)
                   gATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGttt                                                                                                                                 <  1:953304/150‑1 (MQ=255)
                      catcTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTc                                                                                                                              <  1:678202/150‑1 (MQ=255)
                      catcTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTc                                                                                                                              >  2:124783/1‑150 (MQ=255)
                       atcTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTcc                                                                                                                             <  2:913302/150‑1 (MQ=255)
                           gcgcACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTTGTGTAGGCGTTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGGTTTTCCATGc                                                                                                                         <  2:956062/150‑1 (MQ=255)
                                cGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGa                                                                                                                    <  1:892437/150‑1 (MQ=255)
                                 gTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGAt                                                                                                                   <  2:1347274/150‑1 (MQ=255)
                                  ttCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATg                                                                                                                  <  1:2196494/150‑1 (MQ=255)
                                      ggCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACg                                                                                                              >  1:1227479/1‑150 (MQ=255)
                                      ggCACCTGAATACCCAGCGCGTTTTACACCGCCGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACCGGTTTCGGCAAGTTTTTCCATGCCGCGATGGACg                                                                                                              >  2:196569/1‑150 (MQ=255)
                                       gCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGt                                                                                                             >  2:1879768/1‑150 (MQ=255)
                                        cACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTa                                                                                                            >  1:72257/1‑150 (MQ=255)
                                            tGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAAc                                                                                                        >  2:170726/1‑150 (MQ=255)
                                            tGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAAc                                                                                                        <  1:2154659/150‑1 (MQ=255)
                                            tGAATACCCATCGCGTTTTCCACCGACGCGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAAc                                                                                                        >  1:1373096/1‑150 (MQ=255)
                                             gAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACa                                                                                                       <  2:124673/150‑1 (MQ=255)
                                                 aCCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACg                                                                                                   >  1:155048/1‑150 (MQ=255)
                                                  cccATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACgg                                                                                                  >  1:299703/1‑150 (MQ=255)
                                                  cccATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTACCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACgg                                                                                                  <  1:816180/150‑1 (MQ=255)
                                                   ccATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGt                                                                                                 >  2:1905164/1‑150 (MQ=255)
                                                   ccATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGt                                                                                                 >  2:2340139/1‑150 (MQ=255)
                                                     aTCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAg                                                                                               >  1:1377972/1‑150 (MQ=255)
                                                     aTCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAg                                                                                               >  2:1436172/1‑150 (MQ=255)
                                                        gcgTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCg                                                                                            <  2:1220448/150‑1 (MQ=255)
                                                         cgctttCCACCGACGGGGGGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATCGACGTAGAACAGACGGTAGTCgg                                                                                           <  2:129668/147‑1 (MQ=255)
                                                          gTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGc                                                                                          >  2:944127/1‑150 (MQ=255)
                                                               ccACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGa                                                                                     >  1:1950879/1‑150 (MQ=255)
                                                               ccACCGACGTGGTGTAGGCGGTGCTGTGGGAAAAGCCGAAAATCCCGCAAACACGGTCAGAGAGGAAGGTCAATTCGTTATAACCCATACGGGTTTCCGACAGTTTTTCCATGCCTCGATGGACGTAGAAAAGACGTTAGTCGGCGTCGa                                                                                     >  2:1687082/1‑150 (MQ=255)
                                                                  ccGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATaa                                                                                  <  2:541719/150‑1 (MQ=255)
                                                                    gACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGACGGACGTAGAACAGACGGTAGTCGGCGTCGATAATg                                                                                <  1:1263110/150‑1 (MQ=255)
                                                                        tggtgTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGtttt                                                                            >  2:628516/1‑150 (MQ=255)
                                                                         ggtgTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTc                                                                           <  1:1328980/150‑1 (MQ=255)
                                                                         ggtgTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTc                                                                           >  2:570233/1‑150 (MQ=255)
                                                                         ggtgTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTc                                                                           >  2:578927/1‑150 (MQ=255)
                                                                          gtgtAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCg                                                                          <  1:1498109/150‑1 (MQ=255)
                                                                           tgtAGTCGGTGCTGTGGGCAAAGCCTCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGTTGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGc                                                                         <  1:841198/150‑1 (MQ=255)
                                                                           tgtAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGGCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGc                                                                         <  1:1875358/150‑1 (MQ=255)
                                                                            gtAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGcc                                                                        >  1:133880/1‑150 (MQ=255)
                                                                            gtAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGcc                                                                        <  2:299703/150‑1 (MQ=255)
                                                                              aGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCAt                                                                      >  2:2481747/1‑150 (MQ=255)
                                                                               ggCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATc                                                                     >  1:2318502/1‑150 (MQ=255)
                                                                                 cGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTAACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATcga                                                                   >  1:2285265/1‑150 (MQ=255)
                                                                                    tGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATcgacga                                                                <  1:532886/150‑1 (MQ=255)
                                                                                     gCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGaa                                                               >  2:235534/1‑150 (MQ=255)
                                                                                      cTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAAc                                                              <  2:2168399/150‑1 (MQ=255)
                                                                                      cTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAAc                                                              <  2:786513/150‑1 (MQ=255)
                                                                                      cTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAAc                                                              <  1:565957/150‑1 (MQ=255)
                                                                                       tgtgGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACa                                                             >  2:1355698/1‑150 (MQ=255)
                                                                                       tgtgGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACa                                                             <  2:41750/150‑1 (MQ=255)
                                                                                          gggCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTGTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGAc                                                          >  1:2362620/1‑150 (MQ=255)
                                                                                           ggCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACg                                                         >  1:436053/1‑150 (MQ=255)
                                                                                           ggCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACg                                                         >  1:1341973/1‑150 (MQ=255)
                                                                                           ggCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACg                                                         >  1:1678669/1‑150 (MQ=255)
                                                                                            gCAAAGCCGCAGATCCCGCCCACACGGTCAGAGAGGAAGGTCACTACGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACgg                                                        >  2:856697/1‑150 (MQ=255)
                                                                                             cAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGa                                                       <  2:2408995/150‑1 (MQ=255)
                                                                                             cAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGa                                                       >  1:2460856/1‑150 (MQ=255)
                                                                                               aaGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAg                                                     <  2:2225544/150‑1 (MQ=255)
                                                                                                     cAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGccc                                               <  1:1110069/150‑1 (MQ=255)
                                                                                                           ccGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTtca                                         <  1:896733/150‑1 (MQ=255)
                                                                                                           ccGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCAGCAAGTTTTTCCATGCCTCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCTCCATCGACGAACAGACGGAAGTGGCCCGGTtca                                         >  2:1192433/1‑150 (MQ=255)
                                                                                                             gCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTtcatc                                       <  1:1879768/150‑1 (MQ=255)
                                                                                                              cacacaCGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACTGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTtcatca                                      <  2:2470627/150‑1 (MQ=255)
                                                                                                              cacacaCGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTtcatca                                      <  2:2479042/150‑1 (MQ=255)
                                                                                                                cacacGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGa                                    >  2:2059888/1‑150 (MQ=255)
                                                                                                                    cGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGACATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTg                                >  2:1831982/1‑150 (MQ=255)
                                                                                                                    cGGTCAGAGAGGAAGGTCACTTAGTTATAAAAAATACGGGATTCCGCCAGTTTTTACATGCCGCCATGGACGTAGAAAAGAAGGTACTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTCGACCGGTTCATCAGAAGTg                                >  2:465370/1‑150 (MQ=255)
                                                                                                                          gagagGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGc                          >  2:2099282/1‑150 (MQ=255)
                                                                                                                              ggAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCg                      <  1:941266/150‑1 (MQ=255)
                                                                                                                                 aGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGAc                   <  2:2479211/150‑1 (MQ=255)
                                                                                                                                  ggTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGAcc                  <  1:244686/150‑1 (MQ=255)
                                                                                                                                  ggTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGAcc                  <  2:1520394/150‑1 (MQ=255)
                                                                                                                                  ggTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGAcc                  <  1:1547500/150‑1 (MQ=255)
                                                                                                                                    tCACTTCGTTATAACCTATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCaa                >  2:2222992/1‑150 (MQ=255)
                                                                                                                                    tCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCaa                <  2:1613197/150‑1 (MQ=255)
                                                                                                                                        ttCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCgg            <  1:1905164/150‑1 (MQ=255)
                                                                                                                                           gcggaaaCCCGTATGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCac         >  1:22946/6‑150 (MQ=255)
                                                                                                                                            cggaaaCCCGTATGGGTTTCCGCCAGTTTTTCGATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCacg        >  1:2364131/5‑150 (MQ=255)
                                                                                                                                             tataACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCacga       >  1:2259604/1‑150 (MQ=255)
                                                                                                                                              ataACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCacgac      <  2:72257/150‑1 (MQ=255)
                                                                                                                                              ataACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCacgac      <  1:2099105/150‑1 (MQ=255)
                                                                                                                                                 aCCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCACGACatt   >  1:2480536/1‑150 (MQ=255)
                                                                                                                                                  cccATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCACGACgttg  >  2:1459127/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CCTCCAGCAGAATGGCGCGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCACGACGTTG  >  NC_000913/2845640‑2845935

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: