Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,719,650	T→C	V58A (GTA→GCA)	yiaG →	HTH_CROC1 family putative transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,719,650	0	T	C	100.0%	302.2 / NA	79	V58A (GTA→GCA)	yiaG	HTH_CROC1 family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (52/27);  total (52/27)

AGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTT  >  NC_000913/3719503‑3719792                                                                                                                                                    |                                                                                                                                               aGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATc                                                                                                                                             >  1:959765/1‑151 (MQ=255)            agcCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATgg                                                                                                                                  <  2:385784/151‑1 (MQ=255)              ccTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTa                                                                                                                                <  1:26607/151‑1 (MQ=255)                   caCAGATTGATTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGCAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAgg                                                                                                                            >  2:305244/2‑150 (MQ=255)                   aaCAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAgg                                                                                                                            >  1:179083/1‑150 (MQ=255)                    aCAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGa                                                                                                                           >  2:1016298/1‑150 (MQ=255)                    aCAGACTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGa                                                                                                                           >  1:59993/1‑150 (MQ=255)                      aGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCTCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAACGGa                                                                                                                           >  2:1053689/1‑148 (MQ=255)                      aGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAAt                                                                                                                         >  2:474327/1‑150 (MQ=255)                      aGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATg                                                                                                                        >  1:378951/1‑151 (MQ=255)                        aTTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATggg                                                                                                                      >  2:265453/1‑151 (MQ=255)                        aTTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATggg                                                                                                                      >  2:201298/1‑151 (MQ=255)                        aTTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATggg                                                                                                                      >  2:885967/1‑151 (MQ=255)                           gTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGaa                                                                                                                    <  2:239943/150‑1 (MQ=255)                           gTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAAt                                                                                                                   >  2:956398/1‑151 (MQ=255)                            tttttAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATc                                                                                                                  >  1:615429/1‑151 (MQ=255)                               ttAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCa                                                                                                                <  1:79673/150‑1 (MQ=255)                                   aGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGAcg                                                                                                            >  2:206821/1‑150 (MQ=255)                                    gATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGAcgcg                                                                                                          >  1:82816/1‑151 (MQ=255)                                     aTGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGt                                                                                                         >  2:134018/1‑151 (MQ=255)                                      tGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTg                                                                                                        >  1:58829/1‑151 (MQ=255)                                       gAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTg                                                                                                        >  2:1030768/1‑150 (MQ=255)                                          aCGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGAGTGAAg                                                                                                     >  2:940229/1‑150 (MQ=255)                                           cGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGcc                                                                                                   >  2:619821/1‑151 (MQ=255)                                            gcagaaaaagacggggACGCACAGAACAACGGCCGGGACCGAAAAGGAGCAGGGACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGAGAGAAGCCt                                                                                                  >  2:512081/16‑151 (MQ=255)                                            gCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCt                                                                                                  >  1:75977/1‑151 (MQ=255)                                             cAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCtt                                                                                                 <  2:857392/151‑1 (MQ=255)                                             cAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCtt                                                                                                 <  2:266308/151‑1 (MQ=255)                                              aGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTc                                                                                                >  1:218954/1‑151 (MQ=255)                                              aGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTc                                                                                                <  2:630065/151‑1 (MQ=255)                                                aaaaTTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCa                                                                                               >  2:150716/1‑150 (MQ=255)                                                  caTTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGt                                                                                            >  1:238821/2‑151 (MQ=255)                                                      aCCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCg                                                                                        <  2:179077/151‑1 (MQ=255)                                                       cccTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCg                                                                                        >  2:519403/1‑150 (MQ=255)                                                           gACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACt                                                                                    <  2:861690/150‑1 (MQ=255)                                                            aCGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTa                                                                                   >  2:102026/1‑150 (MQ=255)                                                                acaGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAAtt                                                                              >  2:909669/1‑151 (MQ=255)                                                                 caGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCACGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTg                                                                             >  2:238233/1‑151 (MQ=255)                                                                  aGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTg                                                                             >  1:554610/1‑150 (MQ=255)                                                                  aGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTg                                                                             <  2:760129/150‑1 (MQ=255)                                                                          gTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTg                                                                    <  2:96400/151‑1 (MQ=255)                                                                            ccTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGAt                                                                  <  1:250863/151‑1 (MQ=255)                                                                            ccTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGAt                                                                  <  2:468791/151‑1 (MQ=255)                                                                             cTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGAtt                                                                 >  1:212501/1‑151 (MQ=255)                                                                              tGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTc                                                                >  1:906352/1‑151 (MQ=255)                                                                               gTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGAtt                                                                 <  2:58828/149‑1 (MQ=255)                                                                               gTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTc                                                                >  1:863663/1‑150 (MQ=255)                                                                               gTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCa                                                               <  1:157314/151‑1 (MQ=255)                                                                                      aaTTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAAc                                                        <  2:825104/151‑1 (MQ=255)                                                                                       aTTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAAcc                                                       >  1:194100/1‑151 (MQ=255)                                                                                             cAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCAt                                                 <  1:702859/151‑1 (MQ=255)                                                                                             cAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCAt                                                 >  1:1026506/1‑151 (MQ=255)                                                                                              aGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCAtt                                                >  1:158584/1‑151 (MQ=255)                                                                                                 tACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATtaa                                              <  2:136627/150‑1 (MQ=255)                                                                                                  aCGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATtaag                                             <  1:30434/150‑1 (MQ=255)                                                                                                      aaaGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGc                                        >  1:525205/1‑151 (MQ=255)                                                                                                       aaGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGc                                        >  1:915852/1‑150 (MQ=255)                                                                                                       aaGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCa                                       <  1:390427/151‑1 (MQ=255)                                                                                                        aGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGc                                        >  2:964963/1‑149 (MQ=255)                                                                                                          ggACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGtt                                    <  2:313827/151‑1 (MQ=255)                                                                                                           gACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTg                                   <  2:860537/151‑1 (MQ=255)                                                                                                            aCAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTg                                   >  2:50791/1‑150 (MQ=255)                                                                                                            aCAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGa                                  <  2:480932/151‑1 (MQ=255)                                                                                                              aGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATg                                >  2:357882/1‑151 (MQ=255)                                                                                                                gATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATgg                               >  1:728217/1‑150 (MQ=255)                                                                                                                gATTAAAAAGCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGa                              <  2:723839/151‑1 (MQ=255)                                                                                                                  ttAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAAt                            >  1:682923/1‑151 (MQ=255)                                                                                                                   tAAAAATCGATGATTTCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAAt                            <  1:210524/150‑1 (MQ=255)                                                                                                                                tttCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCAc              >  1:119736/1‑151 (MQ=255)                                                                                                                                 ttCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACt             >  1:129382/1‑151 (MQ=255)                                                                                                                                  tCGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACtt            >  1:317445/1‑151 (MQ=255)                                                                                                                                   cGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACttt           >  1:958904/1‑151 (MQ=255)                                                                                                                                   cGCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACttt           >  1:931985/1‑151 (MQ=255)                                                                                                                                    gCCCGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTa          <  2:891181/151‑1 (MQ=255)                                                                                                                                     cccGGGTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTAt         >  2:942874/1‑151 (MQ=255)                                                                                                                                          gTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTg    >  2:655914/1‑151 (MQ=255)                                                                                                                                          gTTTTGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTg    <  1:176386/151‑1 (MQ=255)                                                                                                                                            tttGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGtt  >  2:120349/1‑151 (MQ=255)                                                                                                                                            tttGGGCGCATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGtt  >  1:5006/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                               AGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTT  >  NC_000913/3719503‑3719792

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