Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,903	A→G	I216T (ATC→ACC)	rph ←	truncated RNase PH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,903	0	A	G	98.6%	240.6 / NA	71	I216T (ATC→ACC)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (41/29);  total (41/30)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCC‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCAC  >  NC_000913/3815664‑3816147                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     gaagaaATAGGTGCTTTTTCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACTCCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATAGGGTGTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGct                                                                                                                                                                                                                                             <  2:446974/250‑1 (MQ=255)         aGGGGCTTTTGCGCCCGGNTTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGt                                                                                                                                                                                                                                                 <  2:238927/238‑1 (MQ=255)         aGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGt                                                                                                                                                                                                                                                 >  1:238927/1‑238 (MQ=255)         aGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCCTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATgg                                                                                                                                                                                                                                     <  1:183833/250‑1 (MQ=255)            ggCTTTTGCGCCCGGATTTCAGCGTAAACACNCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCAGCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGa                                                                                                                                                                                                                                 <  2:249792/251‑1 (MQ=255)              cTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAACGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAAc                                                                                                                                                                                                                               >  1:472566/1‑251 (MQ=255)                tttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACg                                                                                                                                                                                                                              <  1:428267/250‑1 (MQ=255)                tttGCGACCGGATTTCCGCGTAAACTCGCAAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGGGAACgg                                                                                                                                                                                                                             >  1:429917/1‑251 (MQ=255)                  ttcgcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCAC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGc                                                                                                                                                                                                                            <  2:35436/248‑1 (MQ=255)                  tGCGCCCGTATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCTCAAATTCAATAAGCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAGAGGCGGCTCATCAGGCGCATTAAAAAGCAGTTTGCCAGCGCCGCCTTGTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGGCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGc                                                                                                                                                                                                                            <  2:143015/250‑1 (MQ=255)                                   gTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTcccc                                                                                                                                                                                                          >  1:91726/1‑251 (MQ=255)                                   gTAAACCCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTcccc                                                                                                                                                                                                          <  2:91726/251‑1 (MQ=255)                                         ccgcgtaacgttaccacCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACt                                                                                                                                                                                                    <  1:253298/237‑1 (MQ=255)                                           gCCAAACTTTAACACCTGCTTGCTAAGCGCCAATTCAATAACCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAACAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACACGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTc                                                                                                                                                                                                  >  2:480647/1‑251 (MQ=255)                                             cAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGACTTAGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGTGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCa                                                                                                                                                                                                 >  2:402158/1‑250 (MQ=255)                                                    tAACACCTGCTTGCTAAGCGCAAATTCAATTAACGGGCGCTGATATGGTTTCATGCCTTCGCTCCTCAACGTACTTTTCGACAGACAAAAAAAAGGCGACACATCAGTCGCCTTAAAAATCAGTTTTCCAGCGCCGCCTTCTGCGTCGCGCCAATGTATTCGAGTCC‑CCC‑TC‑GGGCCAGCGCCAACAAGGGGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTACCCTGCACTGCAATGATgc                                                                                                                                                                                          <  2:335057/250‑1 (MQ=255)                                                     aaCACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATgcgc                                                                                                                                                                                        <  1:372187/251‑1 (MQ=255)                                                                gcaacgcgaaaattcaagatccttgcgcTGAAATGTTTTCATGCCTTCGCTCCGCATCTTACTTTTCTACAGACAAAAAAAAGGCGACGCAGCAGTCTCCGTAAAAATCAGTTTGCCATCGCCGCCATCTGCTTCGCTACAATGGAGTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTATCTCTTCATGGGTGAACGGCTCGCCTTCTGTCGTCCCCTGCACTTCTATGATGCGCCCGTCTTCGGt                                                                                                                                                                             <  2:282294/227‑1 (MQ=255)                                                                                  aaaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGCCGCTAAAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATAACTACGTTCATGTc                                                                                                                                                           >  2:178271/1‑251 (MQ=255)                                                                                  aaaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCCCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTc                                                                                                                                                           <  1:276112/251‑1 (MQ=255)                                                                                       ggCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGtct                                                                                                                                                      >  1:336019/1‑251 (MQ=255)                                                                                          gcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGGCTTCGGTCATCACTACGTTCATGTCGGTCTCTg                                                                                                                                                   >  1:480168/1‑251 (MQ=255)                                                                                            tGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCg                                                                                                                                                 >  1:267505/1‑251 (MQ=255)                                                                                                   gTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACCCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGt                                                                                                                                          >  1:322144/1‑251 (MQ=255)                                                                                                      tCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCCTCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCCCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGCGTCtt                                                                                                                                       >  1:260646/1‑251 (MQ=255)                                                                                                            cTTCGCTACTCATCTTACTTTTCTACAGACAACAAACAGGCGACTCATCAGTCGCCTTAACAATCAGTTTGCCAGAGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGACAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCCATGCTGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCCGAGGCTTCAACgg                                                                                                                                 >  2:206813/1‑250 (MQ=255)                                                                                                               cGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTAtt                                                                                                                              >  1:414893/1‑251 (MQ=255)                                                                                                                      cATCTTACTTTTCTACAGACAAAAAAAAGNCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCCGCGCCGCCGTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACTTATTCCAGa                                                                                                                         <  2:249094/249‑1 (MQ=255)                                                                                                                               tttCTACAGACAAAAAAAAGGCGCCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCaaatt                                                                                                                <  1:117613/249‑3 (MQ=255)                                                                                                                                ttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCaaatttct                                                                                                             >  2:322001/1‑246 (MQ=255)                                                                                                                                  cTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCaaattgcttc                                                                                                           >  2:511074/1‑244 (MQ=255)                                                                                                                                    aCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCaaattgcttcac                                                                                                         <  1:389571/251‑10 (MQ=255)                                                                                                                                    aCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCaaattgcttc                                                                                                           <  1:471834/249‑8 (MQ=255)                                                                                                                                       gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGCACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGt                                                                                                      >  2:182970/1‑251 (MQ=255)                                                                                                                                       gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGCCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTATTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGt                                                                                                      >  1:155674/1‑251 (MQ=255)                                                                                                                                          aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGGCTTCAACGTATTCCAGATCGCAAATTGCTTCCCCGttt                                                                                                    >  2:422116/1‑249 (MQ=255)                                                                                                                                           aaacaaaGGCGCCGCATCAGCCGCCTTAAAAATCAGTTGGCCAGCGCCGTCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCGC‑GGGCCAGCTCCAACAAGGTGAGTAGCTCTTCATGGGTTCACGGCTCGCCGTCTGCCATCCCCTGCACTGCAATTATGCGCCCGGCTTCGGTCATCACTCCGTGCATGTCGGTCTCTGCGGCTTATTCTTCAACGGATTCCAGATCGCAAAGTGCTTCACCGtttac                                                                                                  <  2:210915/247‑4 (MQ=255)                                                                                                                                               aaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTCCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTCACCC‑TC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGCAATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAAt                                                                                               >  1:183388/1‑250 (MQ=255)                                                                                                                                                aaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAAtac                                                                                             >  1:262752/1‑249 (MQ=255)                                                                                                                                                  ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCg                                                                                           <  1:440499/251‑1 (MQ=255)                                                                                                                                                  ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCg                                                                                           <  1:344055/251‑1 (MQ=255)                                                                                                                                                    cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGAc                                                                                         >  1:197973/1‑251 (MQ=255)                                                                                                                                                     gaatcataaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCATC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCCACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGAc                                                                                         >  1:293667/4‑250 (MQ=255)                                                                                                                                                              gTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCg                                                                               >  1:510088/1‑251 (MQ=255)                                                                                                                                                                cGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGACCGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGGCTCTGCGGCTGAGTCTTCAACGTATTCCCGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGc                                                                             >  2:384908/1‑251 (MQ=255)                                                                                                                                                                cGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCg                                                                                                                                                 >  2:229719/1‑183 (MQ=255)                                                                                                                                                                cGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCg                                                                                                                                                 <  1:229719/183‑1 (MQ=255)                                                                                                                                                                    ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGAATCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGGCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCACATTGCTTCACCGTTTACAATACCGACAGAACCTGCGGCTaac                                                                         >  2:303507/1‑249 (MQ=255)                                                                                                                                                                      aaaaaTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGCTTACAATACCGACAGACACTGCGGCTAccct                                                                       >  2:142792/1‑249 (MQ=255)                                                                                                                                                                        aaaTCAGTTTGCAAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCACAAAGGCGAGTAGCTCTTGATGGGTGAACGGCTCGCCTTCTTCCGTCCCCTGCCCTTTAATTGTGCCCCCATCGATGGTGTTCACTACCTTCATGTCGGTCTATTCTGCTGACTCTTCAATGTATCGCAGGTCCCAAATTGCTTCAGCGTTNACAATCCCGCCAGAAACTGCGGccaccaaac                                                                     >  2:255360/1‑248 (MQ=255)                                                                                                                                                                                 gtcccgcgcCGCATTCTGCGTCGCTACAATGTATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCAAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCGGCTCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCTAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCAGCCCCTTCATCg                                                            <  2:488228/246‑1 (MQ=255)                                                                                                                                                                                  gCCAGCGCCGCCTTCTTCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGCCATCACTACGTTCATGTCGGCCTCTGCGGCTGAGGCTTCACCGTATTCCAGATCGCAACTTGCTTCCCCGTTTACACTACCTACAGAAACTGCGGCTACCACCCCTTTCCTCgg                                                           >  2:66177/1‑251 (MQ=255)                                                                                                                                                                                     aGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCCCTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGAAAGAAACTGCCGCCACCATCCCTTTCATCggggt                                                        >  2:285705/1‑248 (MQ=255)                                                                                                                                                                                     aGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCgggtt                                                        <  1:322001/251‑4 (MQ=255)                                                                                                                                                                                         ccgccTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGCTGCGCCCGTCTTCGGTCATCACTCCGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAACCTGCGGCTACCCTCCCTTTTATCCggttggtt                                                    >  2:37546/1‑251 (MQ=255)                                                                                                                                                                                         ccgccTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCCTCGggttggtt                                                    >  1:308149/1‑251 (MQ=255)                                                                                                                                                                                                cTGCGTCGCTACAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCt                                             >  2:378524/1‑251 (MQ=255)                                                                                                                                                                                                          acaagagctgcgatgcc‑ccc‑TC‑GGGCCAGCGCCAACATGGGTAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCTCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTTTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCTTCCCTTTCAACGGGTTGGTTTTCAGCTTGCCGGTTTc                                   <  2:366404/235‑1 (MQ=255)                                                                                                                                                                                                          aCAAGGGATTCGA‑TTC‑CCC‑CCCGGGCCAGCGCCAACCAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCGCTTCAATGAGGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCTGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATATCGACAGAAACTGCGGCTACAATCCCTTTCATCGGGTTGGTGTTCAGCTTGCCGtttt                                    <  2:193/250‑1 (MQ=255)                                                                                                                                                                                                           cAATGGATTCGA‑TTC‑CCCCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAACCTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTTCCGTTTTcc                                  >  2:19336/1‑251 (MQ=255)                                                                                                                                                                                                                   tcga‑ttc‑CCCCTC‑GGGCCAGCGCCAACAAGGTGCGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGACACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCCCCAGCtt                          >  2:64700/1‑251 (MQ=255)                                                                                                                                                                                                                                 c‑GGGCCAGCGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCGTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGATCATCTCTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGCCAGAAACTGCGTCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATc              <  1:10630/251‑1 (MQ=255)                                                                                                                                                                                                                                    ggCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCt                                             >  2:53690/1‑218 (MQ=255)                                                                                                                                                                                                                                    ggCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCt                                             <  1:53690/218‑1 (MQ=255)                                                                                                                                                                                                                                    ggCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCtg            >  1:45310/1‑250 (MQ=255)                                                                                                                                                                                                                                     gCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCCATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCtgc           >  1:314931/1‑249 (MQ=255)                                                                                                                                                                                                                                     gCCAGCGCCAAAAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCGGACTTCCCCTGCACTTCAATCATGCGCCCGTCTTCGGGCATCACTACGTTCATGTCGGTCTCGGCGGCTGAGTCTTCAACGGATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTGCCGCCAGCTTCTGTAGCGCATCtgc           <  2:440436/251‑3 (MQ=255)                                                                                                                                                                                                                                        agcgCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCtgcc          >  1:126448/4‑246 (MQ=255)                                                                                                                                                                                                                                        agcgCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCAGCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTGCTGTAGCGCATCTGCCAg        <  1:127587/248‑1 (MQ=255)                                                                                                                                                                                                                                          cgCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAgcg      >  1:209484/2‑251 (MQ=255)                                                                                                                                                                                                                                              aacaaGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCAc  <  1:511074/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     GAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTC‑CCC‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCAC  >  NC_000913/3815664‑3816147

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