Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,540,647 G→C 6.2% R648R (CGC→CGG narZ ← nitrate reductase Z subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,540,6470GC 6.2% 542.4 / 19.4 162R648R (CGC→CGGnarZnitrate reductase Z subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (73/79);  new base C (4/6);  total (77/85)
Fisher's exact test for biased strand distribution p-value = 7.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.39e-01

Only 100 of 167 total aligned reads displayed.

TGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCATTTGGCGCGGTGGTCGGTT  >  NC_000913/1540498‑1540796
                                                                                                                                                     |                                                                                                                                                     
 gttCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGgcgc                                                                                                                                                      >  1:1380637/1‑150 (MQ=255)
 gttCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGgcgc                                                                                                                                                      <  1:1496547/150‑1 (MQ=255)
    cGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGaa                                                                                                                                                   >  1:284871/1‑150 (MQ=255)
    cGCAGGCCATACGTAAATCGCCCGATTTCACCGCCTGCGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGaa                                                                                                                                                   >  2:558455/1‑150 (MQ=255)
        ggCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCAGTGGGGGATAATCCGGCCTTGTCCGCATCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATAAGTTCGGCGCCAACGta                                                                                                                                               >  1:1633040/1‑149 (MQ=255)
         gCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTg                                                                                                                                              >  2:141868/1‑150 (MQ=255)
          ccATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGa                                                                                                                                             <  2:663551/150‑1 (MQ=255)
          ccATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGCGCCGAACGtatg                                                                                                                                             <  1:2195615/150‑4 (MQ=255)
            aTACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATTCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGaaa                                                                                                                                           <  1:860482/150‑1 (MQ=255)
            aTACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGaaa                                                                                                                                           <  2:427027/150‑1 (MQ=255)
            aTACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGaaa                                                                                                                                           <  2:1293129/150‑1 (MQ=255)
            aTACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGCGCCGAACGtatggg                                                                                                                                           <  1:103018/150‑6 (MQ=255)
            aTACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGCGCCGAACGtatggg                                                                                                                                           <  2:95433/150‑6 (MQ=255)
              aCGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAaatc                                                                                                                                         >  2:1559730/1‑150 (MQ=255)
              aCGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAaatc                                                                                                                                         >  2:1155506/1‑150 (MQ=255)
                 tAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACTCCCCAGCTGCGTCGCATAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAaatcaat                                                                                                                                      <  2:1527270/150‑1 (MQ=255)
                 tAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAaatcaat                                                                                                                                      <  1:1172312/150‑1 (MQ=255)
                  aaaTCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAaatcaatc                                                                                                                                     >  1:2213205/1‑150 (MQ=255)
                  aaaTCGCCCGATTTCAGCGCCTGGGCGGTAAATTCCGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAaatcaatc                                                                                                                                     >  1:413084/1‑150 (MQ=255)
                   aaTCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAatcaatca                                                                                                                                    <  2:517250/150‑1 (MQ=255)
                   aaTCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAatcaatca                                                                                                                                    <  2:1163745/150‑1 (MQ=255)
                    aTCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAg                                                                                                                                   <  2:2213205/150‑1 (MQ=255)
                     tCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGa                                                                                                                                  <  2:1246933/150‑1 (MQ=255)
                       gCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATg                                                                                                                                >  1:54153/1‑150 (MQ=255)
                       gCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATg                                                                                                                                >  2:432819/1‑150 (MQ=255)
                       gCCCGATTTCAGCGCCTGGGCGATAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCCGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATg                                                                                                                                >  2:1095828/1‑150 (MQ=255)
                        cccGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGa                                                                                                                               >  2:157316/1‑150 (MQ=255)
                             tttCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCt                                                                                                                          >  1:1712673/1‑150 (MQ=255)
                                cAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGta                                                                                                                       >  2:2336678/1‑150 (MQ=255)
                                cAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGCCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGATGCGGCGCAGAAGGTCGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCAGta                                                                                                                       >  2:377345/1‑150 (MQ=255)
                                 aGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGtac                                                                                                                      <  2:201715/150‑1 (MQ=255)
                                 aGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGtac                                                                                                                      <  1:2058477/150‑1 (MQ=255)
                                  gcgcCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGtact                                                                                                                     <  2:216830/150‑1 (MQ=255)
                                  gcgcCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGtact                                                                                                                     >  2:650246/1‑150 (MQ=255)
                                     ccTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATACCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCAATACGTTAGGCGCGAAAGTTGAAATCAATAATATGACCGCTGtacttaa                                                                                                                  >  2:77056/1‑149 (MQ=255)
                                        gggCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTgg                                                                                                               <  1:2165880/150‑1 (MQ=255)
                                        gggCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTgg                                                                                                               <  1:1545760/150‑1 (MQ=255)
                                           cGGTAAATTCTGTGGGGTATAATCCGGCCTTGTCGGCTTCAGCTTTAATAACGCGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCCGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCAt                                                                                                            >  2:1938929/1‑150 (MQ=255)
                                            ggTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATc                                                                                                           >  1:911004/1‑150 (MQ=255)
                                                     tgtgGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCgg                                                                                                  <  1:6932/150‑1 (MQ=255)
                                                       tgGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTg                                                                                                >  2:1880580/1‑150 (MQ=255)
                                                           ggATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAg                                                                                            >  2:1101669/1‑150 (MQ=255)
                                                           ggATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAg                                                                                            >  2:1137924/1‑150 (MQ=255)
                                                               aaTCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAAc                                                                                        <  2:2081588/150‑1 (MQ=255)
                                                                  ccGGCCTTGTCGGCTTCAGCTTTAATCCCGAGTGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTcc                                                                                     >  2:1991935/1‑150 (MQ=255)
                                                                  ccGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTcc                                                                                     <  1:1800608/150‑1 (MQ=255)
                                                                            cGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAgaga                                                                           >  1:2137920/1‑150 (MQ=255)
                                                                            cGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAgaga                                                                           >  2:1704113/1‑150 (MQ=255)
                                                                                  cAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACttttt                                                                     <  1:1039783/150‑1 (MQ=255)
                                                                                  cAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACttttt                                                                     <  1:2485385/150‑1 (MQ=255)
                                                                                   aGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTc                                                                    <  1:641047/150‑1 (MQ=255)
                                                                                    gCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCa                                                                   <  1:780788/150‑1 (MQ=255)
                                                                                    gCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCa                                                                   <  1:1216646/150‑1 (MQ=255)
                                                                                    gCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCa                                                                   >  2:1980919/1‑150 (MQ=255)
                                                                                        tAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAgc                                                               <  2:1379434/150‑1 (MQ=255)
                                                                                         aaTCCCGAGCGGGTAACGCCCGAGCTGCGGCGCATAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTTTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAgcg                                                              <  1:1403329/150‑1 (MQ=255)
                                                                                            cccGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCa                                                           >  1:2282312/1‑150 (MQ=255)
                                                                                             ccGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCt                                                                                    <  2:2358734/124‑1 (MQ=255)
                                                                                             ccGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCt                                                                                    >  1:2358734/1‑124 (MQ=255)
                                                                                                   gggTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTg                                                    <  2:813138/150‑1 (MQ=255)
                                                                                                    ggTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTgg                                                   >  2:2353633/1‑150 (MQ=255)
                                                                                                    ggTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTgg                                                   <  1:2397820/150‑1 (MQ=255)
                                                                                                      ctgcgccgcAGCTGGGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGaa                                                 >  2:2181746/4‑150 (MQ=255)
                                                                                                       tACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAAt                                                >  1:364317/1‑150 (MQ=255)
                                                                                                       tACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAAt                                                <  1:418484/150‑1 (MQ=255)
                                                                                                       tACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAAt                                                <  1:81096/150‑1 (MQ=255)
                                                                                                        aCGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATg                                               >  2:711089/1‑150 (MQ=255)
                                                                                                         cGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGa                                              >  1:710695/1‑150 (MQ=255)
                                                                                                           ccccAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGAtt                                            >  1:8839/1‑150 (MQ=255)
                                                                                                           ccccAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGAtt                                            >  1:671726/1‑150 (MQ=255)
                                                                                                            cccATCTGCGGCTCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTg                                           <  2:861419/150‑1 (MQ=255)
                                                                                                            cccAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTg                                           <  1:141868/150‑1 (MQ=255)
                                                                                                               aGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAg                                        >  1:1273034/1‑150 (MQ=255)
                                                                                                                     ggcgCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAAc                                  >  1:2100131/1‑150 (MQ=255)
                                                                                                                      gcgcAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACg                                 >  1:290577/1‑150 (MQ=255)
                                                                                                                       cgcAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGa                                <  2:1619263/150‑1 (MQ=255)
                                                                                                                       cgcAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGa                                <  1:1491843/150‑1 (MQ=255)
                                                                                                                        gcAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAg                               <  2:303538/150‑1 (MQ=255)
                                                                                                                         cAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAgg                              <  2:1198892/150‑1 (MQ=255)
                                                                                                                         cAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAgg                              <  1:499602/150‑1 (MQ=255)
                                                                                                                         cAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAgg                              <  1:2343248/150‑1 (MQ=255)
                                                                                                                          aGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTGCTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGt                             <  2:936202/150‑1 (MQ=255)
                                                                                                                          aGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGt                             >  2:223126/1‑150 (MQ=255)
                                                                                                                           gAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTg                            >  2:386736/1‑150 (MQ=255)
                                                                                                                            aaGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGc                           <  1:991177/150‑1 (MQ=255)
                                                                                                                            aaGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGc                           <  2:290577/150‑1 (MQ=255)
                                                                                                                             aGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCTGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACTAGGTGCt                          <  1:1938929/150‑1 (MQ=255)
                                                                                                                             aGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCt                          >  1:1401303/1‑150 (MQ=255)
                                                                                                                              ggTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTg                         >  2:872046/1‑150 (MQ=255)
                                                                                                                              ggTAGCCAGCCCATAAGTTCGCCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTg                         >  2:393948/1‑150 (MQ=255)
                                                                                                                                  gccagccCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCa                     >  2:636730/1‑150 (MQ=255)
                                                                                                                                  gccagccCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCa                     <  1:432819/150‑1 (MQ=255)
                                                                                                                                     gatttcaaCGTTCGCGCCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCAttt                  >  1:1928325/8‑150 (MQ=255)
                                                                                                                                       tttcaaCGTTCGCGCCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCATTTgg                >  2:1841867/6‑150 (MQ=255)
                                                                                                                                       cccATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCATTTgg                >  1:2382819/1‑150 (MQ=255)
                                                                                                                                         tcaaCGTTCGCGCCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTAGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCATTTGgcg              >  2:1296717/4‑150 (MQ=255)
                                                                                                                                              gTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCATTTGGCGCggtg         >  1:2320013/1‑150 (MQ=255)
                                                                                                                                               ttCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGTAATGATTGTAGAAAAACGAGGTGCTGTTCATTTGGCGCggtgg        <  2:1484598/150‑1 (MQ=255)
                                                                                                                                                 cGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCATTTGGCGCGGTGGTc      <  2:1380637/150‑1 (MQ=255)
                                                                                                                                                     gcgAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGTAATGATTGTAGAAAAACGAGGTGCTGTTCATTTGGCGCGGTGGTCGttt  <  2:242015/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCATTTGGCGCGGTGGTCGGTT  >  NC_000913/1540498‑1540796

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: