Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,644,942 C→T intergenic (‑40/+177) rem ← / ← hokD Qin prophage; uncharacterized protein/Qin prophage; small toxic polypeptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,644,9420CT100.0% 624.2 / NA 162intergenic (‑40/+177)rem/hokDQin prophage; uncharacterized protein/Qin prophage; small toxic polypeptide
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (86/76);  total (86/76)

Only 100 of 162 total aligned reads displayed.

TCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCATTTTCAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGACCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCGTTGAGGATGCCTGCCACATCAGAGGT  >  NC_000913/1644795‑1645089
                                                                                                                                                   |                                                                                                                                                   
tCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCATTTTCAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCAt                                                                                                                                                  >  1:430847/1‑151 (MQ=255)
   aaaaCACATAATGCGGTAATGACTAACAAGATTTCCATTTTCAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATcc                                                                                                                                                >  1:1451864/2‑150 (MQ=255)
    aaaCACATAATGCGGTAATGACTAACAAGATTTCCAGTTTCAACGCCCTCAGCATAGAGGCCATCATCGCAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGTATAATTAGCTCTCCCCTTTAAGGGATCAGCCct                                                                                                                                              <  3:353800/151‑1 (MQ=255)
     aaCACATAATGCGGTAATGACTAACAAGATTTCCATTTTCAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCctc                                                                                                                                             >  3:257127/1‑151 (MQ=255)
        acaTAATGCGGTAATGACTAACAAGATTTCCATTTTCAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCtt                                                                                                                                           >  3:21680/1‑150 (MQ=255)
          aTAATGCGGTAATGACTAACAAGATTTCCATTTTCAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATc                                                                                                                                        >  1:1590972/1‑151 (MQ=255)
             aTGCGGTAATGACTAACAAGATTTCCATTTTCAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATTATCCCTCTTATCCCt                                                                                                                                     >  2:27312/1‑151 (MQ=255)
              tGCGGTAATGACTAACAAGATTTCCATTTTCAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTg                                                                                                                                    >  4:470003/1‑151 (MQ=255)
               gCGGTAATGACTAACAAGATTTCCATTTTCAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTg                                                                                                                                    <  1:678317/150‑1 (MQ=255)
                cGGTAATGACTAACAGGATTTCCATTTTCAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTgcg                                                                                                                                  <  1:1679205/151‑1 (MQ=255)
                    aaTGACTAACAAGATTTCCATTTTCAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCt                                                                                                                              <  2:223533/151‑1 (MQ=255)
                      tGACTAACAAGATTTCCATTTTCAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTAc                                                                                                                            >  4:603541/1‑151 (MQ=255)
                        aCTAACAAGATTTCCATTTTCAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACtt                                                                                                                          <  1:1553178/151‑1 (MQ=255)
                           aaCAAGATTTCCATTTTCAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAg                                                                                                                       <  1:1413925/151‑1 (MQ=255)
                             cAAGATTTCCATTTTCAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTa                                                                                                                     <  2:188381/151‑1 (MQ=255)
                             cAAGATTTCCATTTTCAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTa                                                                                                                     <  4:1685550/151‑1 (MQ=255)
                             cAAGATTTCCATTTTCAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTa                                                                                                                     <  1:1030085/151‑1 (MQ=255)
                              aaGATTTCCATTTTCAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTAt                                                                                                                    >  3:135864/1‑151 (MQ=255)
                              aaGATTTCCATTTTCAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTAt                                                                                                                    >  2:355616/1‑151 (MQ=255)
                              aaGATTTCCATTTTCAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTAt                                                                                                                    >  4:241783/1‑151 (MQ=255)
                                   ttCCATTTTCAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTAttttt                                                                                                                <  1:1141319/150‑1 (MQ=255)
                                   ttCCATTTTCAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTAttttt                                                                                                                <  3:27518/150‑1 (MQ=255)
                                     ccATTTTCAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGAt                                                                                                             <  1:170066/151‑1 (MQ=255)
                                     ccATTTTCAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGAt                                                                                                             <  4:1648893/151‑1 (MQ=255)
                                       aTTTTCAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTc                                                                                                           <  4:481757/151‑1 (MQ=255)
                                        ttttCAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCt                                                                                                          >  1:1418664/1‑151 (MQ=255)
                                          ttCAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTAt                                                                                                        >  3:205130/1‑151 (MQ=255)
                                           tCAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTAtt                                                                                                       >  1:1075172/1‑151 (MQ=255)
                                           tCAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTAt                                                                                                        <  4:1831874/150‑1 (MQ=255)
                                            cAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTAtt                                                                                                       <  2:257134/150‑1 (MQ=255)
                                            cAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTAtt                                                                                                       <  3:359292/150‑1 (MQ=255)
                                            cAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTc                                                                                                      >  2:202955/1‑151 (MQ=255)
                                              aCGCACACAGCATAGAGGCCATCTTCACAAAAAGTTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCg                                                                                                    >  4:1544547/1‑151 (MQ=255)
                                               cGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCgg                                                                                                   >  1:134035/1‑151 (MQ=255)
                                                   cacaGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCAcc                                                                                               >  4:180583/1‑151 (MQ=255)
                                                   cacaGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTACTTTTGATTCTATTCCGGCAcc                                                                                               <  2:485200/151‑1 (MQ=255)
                                                    acaGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCg                                                                                              <  4:786993/151‑1 (MQ=255)
                                                        cATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCa                                                                                          <  2:331818/151‑1 (MQ=255)
                                                           agagGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGa                                                                                        <  4:475237/150‑1 (MQ=255)
                                                             agGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACt                                                                                     <  1:618206/151‑1 (MQ=255)
                                                               gCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTc                                                                                   <  1:1565001/151‑1 (MQ=255)
                                                                ccATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTc                                                                                   >  1:1686196/1‑150 (MQ=255)
                                                                ccATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTc                                                                                   >  4:1790165/1‑150 (MQ=255)
                                                                 cATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCaa                                                                                 <  4:650716/151‑1 (MQ=255)
                                                                 cATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCaa                                                                                 <  3:93919/151‑1 (MQ=255)
                                                                 cATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCa                                                                                  >  3:544712/1‑150 (MQ=255)
                                                                  aTCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCaa                                                                                 >  4:717025/1‑150 (MQ=255)
                                                                      tCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAAcgcg                                                                            >  4:1650162/1‑151 (MQ=255)
                                                                      tCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAAc                                                                               <  4:1613434/148‑1 (MQ=255)
                                                                       cacaAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTGTCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGt                                                                           >  4:663307/1‑151 (MQ=255)
                                                                       cacaAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGt                                                                           <  1:84186/151‑1 (MQ=255)
                                                                       cacaAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGt                                                                           >  4:1928635/1‑151 (MQ=255)
                                                                         caAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTg                                                                         >  4:1051109/1‑151 (MQ=255)
                                                                         acaaaaaTTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAAc                                                                               <  1:230167/143‑1 (MQ=255)
                                                                          aaaaaaTTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTg                                                                         >  3:74931/1‑150 (MQ=255)
                                                                          aaaaaaTTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGa                                                                        <  4:105402/151‑1 (MQ=255)
                                                                           aaaaaTTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTAGTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGaa                                                                       <  4:1813844/151‑1 (MQ=255)
                                                                                  ttACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATaaa                                                                 >  2:133313/1‑150 (MQ=255)
                                                                                  ttACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATaaa                                                                 >  1:152285/1‑150 (MQ=255)
                                                                                   tACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATaaaa                                                                <  4:674191/150‑1 (MQ=255)
                                                                                   tACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATaaaa                                                                <  1:1364506/150‑1 (MQ=255)
                                                                                    aCGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATaaaaac                                                              <  4:758133/151‑1 (MQ=255)
                                                                                        aGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATaaaaacaaaa                                                          <  1:1778055/151‑1 (MQ=255)
                                                                                          ttcttcGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATaaaaacaaaaac                                                        >  4:1511126/1‑151 (MQ=255)
                                                                                              tcGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACccg                                                     <  3:166618/150‑1 (MQ=255)
                                                                                                 aTGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACccgccg                                                  >  4:409017/1‑150 (MQ=255)
                                                                                                 aTGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACccgccg                                                  >  1:1020201/1‑150 (MQ=255)
                                                                                                  tGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGa                                                 >  1:75369/1‑150 (MQ=255)
                                                                                                   gTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGa                                                 >  4:531165/1‑149 (MQ=255)
                                                                                                   gTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGa                                                 >  2:352128/1‑149 (MQ=255)
                                                                                                         catcaGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCCCCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGtt                                         <  1:1268832/151‑1 (MQ=255)
                                                                                                         catcaGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGtt                                         >  1:420027/1‑151 (MQ=255)
                                                                                                         catcaGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGt                                          <  1:1740806/150‑1 (MQ=255)
                                                                                                         catcaGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGt                                          <  4:121901/150‑1 (MQ=255)
                                                                                                          atcaGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGACCTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTa                                        <  4:991359/151‑1 (MQ=255)
                                                                                                          atcaGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGtt                                         >  1:510524/1‑150 (MQ=255)
                                                                                                           tcaGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTa                                        >  4:1357191/1‑150 (MQ=255)
                                                                                                                     tCCGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGc                             >  1:450747/1‑151 (MQ=255)
                                                                                                                      ccGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGc                             <  1:1129776/150‑1 (MQ=255)
                                                                                                                      ccGGATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCg                            >  1:656499/1‑151 (MQ=255)
                                                                                                                        ggATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCGtt                          >  3:408588/1‑151 (MQ=255)
                                                                                                                        ggATAATTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGCCGCCGTCCAGACCTTCAACCGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCGt                           <  4:524952/150‑1 (MQ=255)
                                                                                                                            aaTTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCGTTGAgg                      >  1:1810454/1‑151 (MQ=255)
                                                                                                                            aaTTAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCGTTGAgg                      <  2:112391/151‑1 (MQ=255)
                                                                                                                            aaTTAGCTCCCCCCTTTAGGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCGTTGAgg                      <  1:107958/151‑1 (MQ=255)
                                                                                                                               tAGCTCTCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCGTTGAGGAtg                   >  1:352059/1‑151 (MQ=255)
                                                                                                                                  ctctCCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCGTTGAGGAtgcct                >  4:1889648/1‑151 (MQ=255)
                                                                                                                                     tcCCCTTTAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCGTTGAGGAtgcctgcc             >  4:561380/1‑151 (MQ=255)
                                                                                                                                      ccccTTTGAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCGTTGAGGATGCCTGCca            <  1:1887558/151‑1 (MQ=255)
                                                                                                                                           ttAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCGTTGAGGATGCCTGCCACATc        >  3:485396/1‑150 (MQ=255)
                                                                                                                                           ttAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCGTTGAGGATGCCTGCCACATc        >  4:1288337/1‑150 (MQ=255)
                                                                                                                                            tAAGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCGTTGAGGATGCCTGCCACATc        >  4:1882293/1‑149 (MQ=255)
                                                                                                                                             aaGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCTTTGAAAATAAAAACAAAACCCCGCCGAAGCGGGTTAAGTGCGGGTGCGTTGAGGATGCCTGCCACATCa       <  4:1409700/149‑1 (MQ=255)
                                                                                                                                             aaGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCGTTGAGGATGCCTGCCACATc        >  1:329462/1‑148 (MQ=255)
                                                                                                                                             aaGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCGTTGAGGATGCCTGCCACATCa       <  1:571277/149‑1 (MQ=255)
                                                                                                                                             aaGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCGTTGAGGATGCCTGCCACATCa       <  4:431000/149‑1 (MQ=255)
                                                                                                                                             aaGGGATCACCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATCCCGGCCCCGCCCAGAACTTAAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCGTTGAGGATGCCTGCCACATCa       <  1:247661/149‑1 (MQ=255)
                                                                                                                                              tgggATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCGTTGAGGATGCCTGCCACATCag      >  1:764079/2‑149 (MQ=255)
                                                                                                                                              aGGGATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCGTTGAGGATGCCTGCCACATCag      >  4:1434547/1‑149 (MQ=255)
                                                                                                                                                ggATCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCGTTGAGGATGCCTGCCACATCAGAGGt  >  1:112235/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCATTTTCAACGCACACAGCATAGAGGCCATCTTCACAAAAAATTTTACGCAGTTCTTCGATGTTCATCATCAGAATCCTTCCGGATAATTAGCTCTCCCCTTTAAGGGACCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCGTTGAGGATGCCTGCCACATCAGAGGT  >  NC_000913/1644795‑1645089

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: