Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,845,792 C→T 6.6% R227H (CGT→CAT)  hycE ← formate hydrogenlyase subunit HycE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,845,7920CT 6.6% 425.7 / 10.5 136R227H (CGT→CAT) hycEformate hydrogenlyase subunit HycE
Reads supporting (aligned to +/- strand):  ref base C (70/56);  new base T (3/6);  total (73/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 = 1.00e+00

Only 100 of 139 total aligned reads displayed.

GCAGAATGGCGCGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCACGACGTTGTTTTTC  >  NC_000913/2845646‑2845941
                                                                                                                                                  |                                                                                                                                                     
gcagAATGGCGCGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGACAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACggg                                                                                                                                                    <  1:426048/150‑1 (MQ=255)
  agAATGGCGCGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGtt                                                                                                                                                  <  1:36882/150‑1 (MQ=255)
   gAATGGCGCGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGttt                                                                                                                                                 >  1:1511309/1‑150 (MQ=255)
      gtgcgcAGATGATCCGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGCTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCg                                                                                                                                              >  2:2004141/3‑150 (MQ=255)
       ggCGCGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCGTATGGGttataac                                                                                                                                             <  1:2246002/150‑6 (MQ=255)
       ggCGCGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCGTATGGGttataac                                                                                                                                             <  1:1140594/150‑6 (MQ=255)
       ggCGCGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGc                                                                                                                                             >  2:468923/1‑150 (MQ=255)
       ggCGCGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGc                                                                                                                                             >  2:2099105/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)
        gcgcGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGcc                                                                                                                                            <  2:1926736/150‑1 (MQ=255)
             gATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGttt                                                                                                                                       <  1:953304/150‑1 (MQ=255)
              atcatcTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGtttt                                                                                                                                      <  2:2265557/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)
                    tGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATg                                                                                                                                >  2:1110069/1‑150 (MQ=255)
                     gcgcACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTTGTGTAGGCGTTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGGTTTTCCATGc                                                                                                                               <  2:956062/150‑1 (MQ=255)
                           gTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGAt                                                                                                                         <  2:1347274/150‑1 (MQ=255)
                                ggCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACg                                                                                                                    >  1:786513/1‑150 (MQ=255)
                                ggCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACg                                                                                                                    >  1:1227479/1‑150 (MQ=255)
                                 gCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGt                                                                                                                   >  2:1879768/1‑150 (MQ=255)
                                      tGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAAc                                                                                                              <  1:2154659/150‑1 (MQ=255)
                                      tGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAAc                                                                                                              >  2:170726/1‑150 (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:1422960/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                                                                                                     >  2:1436172/1‑150 (MQ=255)
                                               aTCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAg                                                                                                     >  1:1377972/1‑150 (MQ=255)
                                                    gTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGc                                                                                                >  2:944127/1‑150 (MQ=255)
                                                       ttCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTc                                                                                             <  1:2209617/150‑1 (MQ=255)
                                                         ccACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGa                                                                                           >  1:1950879/1‑150 (MQ=255)
                                                          cACCGACGTGGTGTAGCCGGTGCTGTGGGCAAAGCCGCAGAACCCGCACAAACGGTCAGAGAGGAAGGTAACTTAGTTATAACAAATACGGGTTTCCGCCAGTTTTTCCATGCCGAGATGGACGTAGAACAGACGGTAGTCGGCGTCGAt                                                                                          >  2:1298989/1‑150 (MQ=255)
                                                            ccGACTTGGTGTAGGCGGTGCTGTCGTAAAACCCTCAGATCCCGCACACACGGTCAGAGAGGAACGTAACATCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTATAACAGACGGTAGTCGGCTTCGATaa                                                                                        >  2:2044614/1‑150 (MQ=255)
                                                             cGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAAt                                                                                       >  1:922681/1‑150 (MQ=255)
                                                              gACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGACGGACGTAGAACAGACGGTAGTCGGCGTCGATAATg                                                                                      <  1:1263110/150‑1 (MQ=255)
                                                                  tggtgTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGtttt                                                                                  >  2:628516/1‑150 (MQ=255)
                                                                   ggtgTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTc                                                                                 >  2:578927/1‑150 (MQ=255)
                                                                   ggtgTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTc                                                                                 >  2:570233/1‑150 (MQ=255)
                                                                    gtgtAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCg                                                                                <  1:1498109/150‑1 (MQ=255)
                                                                     tgtAGTCGGTGCTGTGGGCAAAGCCTCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGTTGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGc                                                                               <  1:841198/150‑1 (MQ=255)
                                                                     tgtAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGc                                                                               >  2:1751030/1‑150 (MQ=255)
                                                                     tgtAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGGCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGc                                                                               <  1:1875358/150‑1 (MQ=255)
                                                                      gtAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGcc                                                                              <  2:299703/150‑1 (MQ=255)
                                                                      gtAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGcc                                                                              >  1:133880/1‑150 (MQ=255)
                                                                           cGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTAACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATcga                                                                         >  1:2285265/1‑150 (MQ=255)
                                                                              tGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATcgacga                                                                      <  1:532886/150‑1 (MQ=255)
                                                                                cTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAAc                                                                    <  2:2168399/150‑1 (MQ=255)
                                                                                cTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAAc                                                                    <  1:565957/150‑1 (MQ=255)
                                                                                 tgtgGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACa                                                                   <  2:41750/150‑1 (MQ=255)
                                                                                   tgGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGa                                                                 >  1:1317303/1‑150 (MQ=255)
                                                                                    gggCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGAc                                                                <  1:1588582/150‑1 (MQ=255)
                                                                                    gggCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTGTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGAc                                                                >  1:2362620/1‑150 (MQ=255)
                                                                                     ggCAAAGCCGCAGATCCGGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACg                                                               >  2:568099/1‑150 (MQ=255)
                                                                                     ggCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACg                                                               >  1:1341973/1‑150 (MQ=255)
                                                                                      gCAAAGCCGCAGATCCCGCCCACACGGTCAGAGAGGAAGGTCACTACGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACgg                                                              >  2:856697/1‑150 (MQ=255)
                                                                                       cAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCTACAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGa                                                             >  2:1640400/1‑150 (MQ=255)
                                                                                       cAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGa                                                             >  1:2460856/1‑150 (MQ=255)
                                                                                       cAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGa                                                             <  2:2408995/150‑1 (MQ=255)
                                                                                         aaGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAg                                                           <  2:2225544/150‑1 (MQ=255)
                                                                                             cGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGc                                                       >  1:4749/1‑150 (MQ=255)
                                                                                               cAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGccc                                                     <  1:1110069/150‑1 (MQ=255)
                                                                                                 gATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCgg                                                   <  2:262450/150‑1 (MQ=255)
                                                                                                    cccGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCGCGGTtc                                                >  2:2103721/1‑150 (MQ=255)
                                                                                                     ccGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTtca                                               <  1:896733/150‑1 (MQ=255)
                                                                                                     ccGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCAGCAAGTTTTTCCATGCCTCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCTCCATCGACGAACAGACGGAAGTGGCCCGGTtca                                               >  2:1192433/1‑150 (MQ=255)
                                                                                                       gCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTtcatc                                             <  1:1879768/150‑1 (MQ=255)
                                                                                                        cacacaCGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTtcatca                                            <  2:2479042/150‑1 (MQ=255)
                                                                                                             acGGTCAGAGCGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGt                                       <  2:1409122/150‑1 (MQ=255)
                                                                                                              cGGTCAGAGAGGAAGGTCACTTAGTTATAAAAAATACGGGATTCCGCCAGTTTTTACATGCCGCCATGGACGTAGAAAAGAAGGTACTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTCGACCGGTTCATCAGAAGTg                                      >  2:465370/1‑150 (MQ=255)
                                                                                                                gTCAGAGAGGAAGGTCACTTCGTTATAACCCATAGGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGAc                                    <  1:2282815/150‑1 (MQ=255)
                                                                                                                    gagagGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGc                                >  2:2099282/1‑150 (MQ=255)
                                                                                                                        ggAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCg                            <  1:941266/150‑1 (MQ=255)
                                                                                                                        ggAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCg                            <  2:236020/150‑1 (MQ=255)
                                                                                                                           aGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGAc                         <  1:894282/150‑1 (MQ=255)
                                                                                                                           aGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGAc                         <  2:2479211/150‑1 (MQ=255)
                                                                                                                            ggTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGAcc                        <  1:1547500/150‑1 (MQ=255)
                                                                                                                            ggTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGAcc                        <  1:244686/150‑1 (MQ=255)
                                                                                                                            ggTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGAcc                        <  2:1520394/150‑1 (MQ=255)
                                                                                                                              tCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCaa                      >  2:2174951/1‑150 (MQ=255)
                                                                                                                              tCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCaa                      <  2:1613197/150‑1 (MQ=255)
                                                                                                                                 cTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCg                   <  1:2340139/150‑1 (MQ=255)
                                                                                                                                     gcggaaaCCCGTATGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCac               >  1:22946/6‑150 (MQ=255)
                                                                                                                                        ataACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCacgac            <  2:72257/150‑1 (MQ=255)
                                                                                                                                        ataACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCacgac            <  1:2099105/150‑1 (MQ=255)
                                                                                                                                         taACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCacgacg           >  2:1231662/1‑150 (MQ=255)
                                                                                                                                         taACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCacgacg           >  2:839762/1‑150 (MQ=255)
                                                                                                                                          aaCCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCACGACgt          >  1:1892224/1‑150 (MQ=255)
                                                                                                                                           aCCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCACGACatt         >  1:2480536/1‑150 (MQ=255)
                                                                                                                                            cccATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCACGACgttg        >  2:1459127/1‑150 (MQ=255)
                                                                                                                                                 aCGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCACGACGTTGttttt   <  1:2477923/150‑1 (MQ=255)
                                                                                                                                                  cGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCTGACCAATCGGCACGACGTTGTTTTTc  <  2:2460856/150‑1 (MQ=255)
                                                                                                                                                  cGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCACGACGTTGTTTTTc  >  1:69326/1‑150 (MQ=255)
                                                                                                                                                  cGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCACGACGTTGTTTTTc  <  1:356910/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCAGAATGGCGCGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCACGACGTTGTTTTTC  >  NC_000913/2845646‑2845941

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: