Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	433,649	T→G	28.5%	G65G (GGT→GGG)	ribD →	fused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	433,649	0	T	G	28.5%	89.4 / 8.8	73	G65G (GGT→GGG)	ribD	fused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase
Reads supporting (aligned to +/- strand):  ref base T (18/32);  new base G (21/2);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 8.83e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.46e-01

AATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCGCAGGTCG‑CTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCACC  >  NC_000913/433417‑433886                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     aaTTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGCGAcc                                                                                                                                                                                                                                                <  1:205063/239‑1 (MQ=255)                       ccTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAG‑CCAAAGGGGCGACACCCTATG‑TCA‑CACTCGaaa                                                                                                                                                                                                                           >  1:489214/1‑235 (MQ=255)                       ccTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAG‑CCAAAAGGGGGACCGCCTATG‑TCA‑CACTCGAAc                                                                                                                                                                                                                           >  1:118476/1‑236 (MQ=255)                        cTGGAGGCCTACGCCGTGCAGGACGAGTATTCCATGCCCCGGGCGCTAAAGCGGGCGCACCGAGGACTTTTTCCCACCCATCCCACCCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTCCCACCAACGGGCGGGTGAACCACATGCCGAAGTACACGCTTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTg                                                                                                                                                                                                                       <  2:40273/239‑1 (MQ=255)                           gaggaCTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGCGACCGCCCATG‑TCA‑CACTCGAAc                                                                                                                                                                                                                           >  1:353893/1‑232 (MQ=255)                                  aaGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGGAAAAG‑CCAAAAGTGCGACCGCCTATT‑TCA‑CACTCGAAACCTGGAGCcatcat                                                                                                                                                                                                             >  1:328640/1‑239 (MQ=255)                                  aaGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGGAAAAG‑CCAAAAGTGCGACCGCCTATG‑TCA‑CACTCGGAACCTGTAGCcat                                                                                                                                                                                                                >  1:1101638/1‑236 (MQ=255)                                        tGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAG‑CCAAAGGGGCGACCGCCCATG‑TCA‑CACTCGAACCCTGTAGCCATCCTGGTCG‑t                                                                                                                                                                                                      >  1:653211/1‑239 (MQ=255)                                              aCGAGTACTAAGGGGCGCGGGCCGTAACGCTGGCACAACGAGGACTTTTTACCCCGCATCCCACCCAGAAGGTCGGGTGCGTCATTGTAAAAGATGGCGAAATTGTCGGTGAAGGTTACCCCCAACGTGCGGGTGAACCCCATGCCGAAGTACCCGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑Ta                                                                                                                                                                                                     <  2:653211/234‑1 (MQ=255)                                               cgtgttttACAGGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACTTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACCCGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCa                                                                                                                                                                                                <  1:1040566/232‑1 (MQ=255)                                               cGAGTATTACAGGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTAATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTCCACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCa                                                                                                                                                                                                <  1:619418/238‑1 (MQ=255)                                                  gTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAG‑CCAAAAGGGCGACCCCCCATT‑TCA‑CACCCGAACCCTGTAGCCATCATGGTCG‑TACGcccc                                                                                                                                                                                               >  1:917718/1‑234 (MQ=255)                                                   tATTCCATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTCCCACCGACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGCGa                                                                                                                                                                                                                                                  <  1:580850/187‑1 (MQ=255)                                                   tATTACAGGGCCCGGGCACAAAAGCTGGCGCAAGGAGAACTTTTTACCACGCATCCAACCCCGAATGTCGGGTGCGTCATTGTAAAAGATGGCAAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACGTGCCGAAGTACCCGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCcat                                                                                                                                                                                                                <  2:144096/219‑1 (MQ=255)                                                     ttACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCACCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTc                                                                                                                                                                                                                               <  1:1185338/202‑1 (MQ=255)                                                       aCATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGGAAAAA‑CCAAAGGGGGGACCGCCCATG‑GCC‑CACTCCAACCCTGTAGGCATCACCGTCG‑TACGCCACCggtcttt                                                                                                                                                                                       >  1:1044103/1‑233 (MQ=255)                                                        cATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGACCCACATGCCGAAGTACACGCGTTGCGTAGGGC‑GGGTGAAAAAG‑CCAAAGGTGCGa                                                                                                                                                                                                                                                  <  1:275545/182‑1 (MQ=255)                                                             cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCGACGTGCGGGTGAACCACATGCCGAAGTACAAGCGTTGCGTATGGCGGGGGGAAAAAG‑CAAAAGGGGGGACccccctgt                                                                                                                                                                                                                                         >  2:580850/1‑183 (MQ=255)                                                             cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGGAAAAG‑CCAAAAGGGGGAACGCCCCTGTTCA‑CAC‑CCAACCCTGTGTGCATCATGGTCGTTACGCCCCCGGCCtgtagtg                                                                                                                                                                                   >  2:414568/1‑234 (MQ=255)                                                             cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGGGAAAA‑CCAAAAGGGGGGCCGCCCATT‑TCA‑CAACCGGAACCCGTAGACCTCATGTGCG‑GTCGCGACCCTGCcgtggagc                                                                                                                                                                                  >  2:1116983/1‑234 (MQ=255)                                                             cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGGAAAAA‑GCAAAAGTGGGACCGCCCATG‑TTACCACCCAAACCC‑GTAGTCCACATCGTCG‑TCAGACACCCGGCtgtggctc                                                                                                                                                                                  >  2:144577/1‑234 (MQ=255)                                                             cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑Ta                                                                                                                                                                                                     >  1:144096/1‑219 (MQ=255)                                                             cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAA‑CCAAAGGGGCGACCGCCCATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑GACGCCACCCGGCtgtggagc                                                                                                                                                                                  >  1:717164/1‑234 (MQ=255)                                                                  gcgcTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGGGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCCCCCTGCGGTGAGGCACCCa                                                                                                                                                                             >  1:1243701/1‑238 (MQ=255)                                                                  gcgcTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAG‑CCAAAGGGGCGACCGCCCATG‑TCA‑CACTCCAACCCTGTAGCCATCATGGTCG‑TACGCCACCCGGCTGTGACGCAcccc                                                                                                                                                                             >  1:977158/1‑235 (MQ=255)                                                                              gcgcAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GCCGt                                                                                                                                                               >  1:964277/1‑239 (MQ=255)                                                                                    cGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAG‑CCAAAGGGGCGACCGCCCATG‑TCA‑CACTCGAACCCCGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCCCT‑GGCGTAgcgc                                                                                                                                                          >  1:87833/1‑238 (MQ=255)                                                                                                  ctggactCCCAACCGGAATGTGGGGTGCGTCATTGACAAAGATGGCAAAGTTGCGGGAGAGGGTTACCCCCACCGGGCGGGTGAACCACATGCCAAAGTCCCCGCGTGGCGTAGGGC‑GGGTGAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGa                                                                                                                                          <  2:1243701/233‑1 (MQ=255)                                                                                                     gCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAG‑CCAAAGGGGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTTACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGGGGGTGCCTCGATG‑c                                                                                                                                      >  1:100745/1‑239 (MQ=255)                                                                                                      cATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGAGG‑Ca                                                                                                                                     >  1:58672/1‑239 (MQ=255)                                                                                                       aTCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAA‑CCCAAGGGGCGACCCCCCATG‑TCA‑CCCTCCAAACCCGTAGCCATCCAGGTCG‑CACGGCGCTCGGCTCTGGAGCCCTCATCTCCCCTCGTCGTACCGCG‑CGGGTTTGCC‑CGCTG‑Caa                                                                                                                                    >  2:1205209/1‑239 (MQ=255)                                                                                                       aTCCCAACCCGAATGGCGGGTGCGTCATTGTCAAAGAGGGCGAAATTGTCGGGGAAGGTTCCCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑Ca                                                                                                                                     <  1:1205209/238‑1 (MQ=255)                                                                                                                 gAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCCCGATG‑CAAGATCCTaaa                                                                                                                           >  1:782203/1‑237 (MQ=255)                                                                                                                  aaTGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAA‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCCTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGGTGCCTCGATG‑CAAGATCCTAAcc                                                                                                                          >  1:254303/1‑238 (MQ=255)                                                                                                                   aggTCGGGTGCGTCATGGAAAAAGAGGGCAAAATTGGCGGGAAAGGTTACCACCAGCGGGGGGGTGACCCACATGCCAAAGTACACGCGTGGCGTAGGGC‑GGGTAAAAAAG‑CAAAAGGGGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCt                                                                                                                              <  2:1101638/231‑1 (MQ=255)                                                                                                                            gctcaatggccAAAGAGGGCGAATTTGTCGGTGAAGGTTCCCACACACGTGCGGGTGAACCACATGCCGAAGTCCACGCGTTGCGTATGGC‑GGGTAAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCGCAGGTCG‑c                                                                                                               <  1:906836/228‑1 (MQ=255)                                                                                                                             cGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCcatca                                                                                                                                                                                                              >  2:949241/1‑147 (MQ=255)                                                                                                                             cGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCcatca                                                                                                                                                                                                              <  1:949241/147‑1 (MQ=255)                                                                                                                                aagtgaaaaggATGGGAAAATTGCCGGTGAGGGTTCCCACCAAGGGGCGGGTGAACCAAATGCCGAAGTACACGCGTTGCGTAGGGG‑GGGTAAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCGCAGGTCG‑CTGGGc                                                                                                          <  2:489214/229‑1 (MQ=255)                                                                                                                                       aaaGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGGGCGGGTGAACCACATGCCGAAGTACACGCGTGGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCGCAGGTCG‑CTGGGCGTGGACt                                                                                                   <  1:1116983/239‑1 (MQ=255)                                                                                                                                            tgtcaaaaGTGGCGGTGAAGGTTACCACCAACGTGCGGGTGACCCACAGGCGGAAGGACACGCGTGGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCGCAGGTCG‑CTGGGCGTGGACTTTa                                                                                                <  1:272645/232‑1 (MQ=255)                                                                                                                                            tGGCGAATTTGTCGGTAAAGTTTACCACCAACGTGCGGGTAAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCGCAGGTCG‑CTGGGCGTGGACTTTa                                                                                                <  1:66872/237‑1 (MQ=255)                                                                                                                                             ggCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGGGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCGCAg                                                                                                                     >  1:31359/1‑216 (MQ=255)                                                                                                                                               cGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGGAAAAG‑CCCAAAGGGCGACCCCCCATG‑TCA‑CACTCGAACCCTGTAGCCATCATTGGCG‑GACCCCCCCGTGGTGGGGCGCACCCCTCGCCGCT‑GGCGCAGCGCGCCGTGGTTG‑CTCGATGCCAAGTTCCTAA‑CCGCACGTCG‑CCGGGGGTgggcgcttgtg                                                                                             >  2:489609/1‑228 (MQ=255)                                                                                                                                               cGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑Ta                                                                                                                                                                                                     >  2:487655/1‑137 (MQ=255)                                                                                                                                               cGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑Ta                                                                                                                                                                                                     <  1:487655/137‑1 (MQ=255)                                                                                                                                               cGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAG‑CCAAAAGGGGGACCGCCCATG‑TCA‑AACCCGAAACCTGTAGCCATCAAGGTCG‑TACGCCACCGGGCTGGTACGCACTCAACGCCGCC‑GGCGTAGGGCG‑CGCGGGTGCCTCGATG‑CAAAATCCTAACCCGCACGGCG‑CTGGGGGGGGGCGTTAACGt                                                                                            >  2:968789/1‑238 (MQ=255)                                                                                                                                                    ttGCCGGTGAAGGTAACCAGCAACGGGCGGTTGAACCAATTGCCGAAGTCCACGTGTGGCTTAGGGC‑GGTGGAAAAAG‑CCAAGGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCGCAGGTCG‑CTGGGCGTGGACTTTACCGTCTGCaa                                                                                      <  2:254303/239‑1 (MQ=255)                                                                                                                                                      ggaggTGAATGTCACCACAAACGGGCGGATAAACCACTTGCAAAAGTCCACGTGTGGCGTATGGC‑GGGTAAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCGCAGGTCG‑CTGGGCGTGGACTTTACCGTCTGCAACa                                                                                    <  2:100745/236‑1 (MQ=255)                                                                                                                                                                 ttACCACCAGCGTGCGGGTAAACCAAGTGCCGAAGTCCCCGTGTTGCGTAGGGG‑GGGTAAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCGCAGGTCG‑CTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCAtt                                                                          <  2:917718/238‑1 (MQ=255)                                                                                                                                                                 ttACCACCAACGTGGGGTTGAACCACATGCCGAAGTCCACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CAAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCGCAGGTCG‑CTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCAtt                                                                          <  1:968789/238‑1 (MQ=255)                                                                                                                                                                 ttACCACCAACGGGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCGCAGGTCG‑CTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCAt                                                                           <  1:338398/237‑1 (MQ=255)                                                                                                                                                                 tacccaccaACGGGCGGGACACCCCCGTGCCGAGGTACACGCGTCGCGTAGGGC‑GGG‑GAAAAAGCCCAAGGGGGCAACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCGCAGGTCG‑CTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCAt                                                                           <  2:620600/234‑1 (MQ=255)                                                                                                                                                                   accaccAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGCGACCGCCTCTG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCGCAGGTCG‑CTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTTAc                                                                       >  1:59446/1‑239 (MQ=255)                                                                                                                                                                     caccaACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAA‑CCAAAAGGGCGACCGCCTAAT‑TCA‑CACTCGCACCCCGTAGCCATCATGGTCG‑TACGCCACCCGGCTGTTACCCACTCCCCTCCCCC‑TGCGCAGCCGC‑CGTGGCTGGCTCGATG‑CCAAATCCTAACCCCCCAGTCGCCGGGGCGGGGGTTTAACAGTC‑GCAACACGCTGGCACGGACg                                                                      >  2:4475/1‑238 (MQ=255)                                                                                                                                                                     caccaACGTGCGGGTAAACCACATGCCGAAGTCCAGGGGTTGCGTATGGG‑GGGTGAAAAAG‑CCAAGGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCGCAGGTCG‑CTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGt                                                                     <  1:746567/239‑1 (MQ=255)                                                                                                                                                                     cacaatggggCGGGTAACCAACGTGCCGAATTACAGGCGTTGCGTAGGGC‑GGGTGAAAACG‑CAAAGGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCGCAGGTCG‑CTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGt                                                                     <  2:118476/230‑1 (MQ=255)                                                                                                                                                                          aCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGGGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCGCCGGTCG‑CTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGc                                                                 >  1:684632/1‑238 (MQ=255)                                                                                                                                                                          aCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAG‑CCAAAGGTGCGACCGCCCATG‑TCA‑CACTCGAACCCCGTAGCCATCATGGGCG‑TACGCCCCCCTGCTGTGACGCACTCATCGCCGCT‑GGCGTAACGCG‑CGGGGTTGGCTCGCAG‑CAAGATCCTAAACCGCCGGTCG‑CTGGGGGTGGGCTTTACCGTCTGCAACAGGCTGGCGttttatgcagc                                                                 >  2:942440/1‑229 (MQ=255)                                                                                                                                                                             tGCGGGTAACCCACTGGCCGAAGTCCCCGCTTTGCGTATGGC‑GGGTAAAAAAC‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCGCAGGTCG‑CTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGcc                                                                <  1:1201640/236‑1 (MQ=255)                                                                                                                                                                               cGGGTAACCCACAGGCCGAAGTACACGCGTTGCGTAGGGC‑GGGTAAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCGCAGGTCG‑CTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa                                                               <  1:924206/235‑1 (MQ=255)                                                                                                                                                                                            tGCCGAAGTACCCGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCGCAGGTCG‑CTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa                                                               <  1:448166/222‑1 (MQ=255)                                                                                                                                                                                            tGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGGGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCGCCGGTCG‑CTGGGCGGGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCCCGGCCtgatgacgtgc                                               >  1:1018652/1‑233 (MQ=255)                                                                                                                                                                                                      cgcggttgccGTAGGGC‑GGGAAAAAACC‑CAAAAGGTCCTACCGCCTATG‑TCA‑CTCTCGAACCCTGTAGCCATCACGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCGCAGGTCG‑CTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGc                                     <  2:1044103/229‑1 (MQ=255)                                                                                                                                                                                                      cacGCGTTGCGTATGGC‑GGGGGAAAAAG‑CCAAAGGGGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGGGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCCCAGGGCG‑CTGGGCGTGGGCTTTACCGTCTGCAACAGGCTGGGATTGAAGTCAGGCCCGGcccggtg                                                     >  2:448166/1‑217 (MQ=255)                                                                                                                                                                                                               cGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCGCAGGTCG‑CTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGc                                                           <  1:1006143/207‑1 (MQ=255)                                                                                                                                                                                                                          gtaaaaaaG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑Ta                                                                                                                                                                                                     <  2:512854/60‑1 (MQ=255)                                                                                                                                                                                                                          gTGAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑Ta                                                                                                                                                                                                     >  1:512854/1‑63 (MQ=255)                                                                                                                                                                                                                          gTGAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCGCAGGTCG‑CTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGGGAAGCCGAGCAATTGAATa                            >  1:158104/1‑228 (MQ=255)                                                                                                                                                                                                                          gTGAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGGCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGAACCTAACCCCCAGGGCG‑CTGGGGGGGGGCCTTTCCGCCTGCAAAAGGCTGGCATTGACGTCCGCCACGGCCcgatgagta                                                 >  2:1006143/1‑204 (MQ=255)                                                                                                                                                                                                                            gAAAAAG‑CCAAGGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCGCAGGTCG‑CTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCa                <  1:414568/238‑1 (MQ=255)                                                                                                                                                                                                                                 aG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCGCAGGTCG‑CTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGAAGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGa          >  1:663025/1‑239 (MQ=255)                                                                                                                                                                                                                                    ccAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAGCCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCGCAGGTCG‑CTGGGCGGGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAAAACAGGCTTTATCCAGCTGAt         >  2:974304/1‑238 (MQ=255)                                                                                                                                                                                                                                        agttGCCACCGCCTAGG‑TCC‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCGCAGGTCG‑CTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCa    <  2:125835/236‑1 (MQ=255)                                                                                                                                                                                                                                           tGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCGCAGGTCG‑CTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCAcc  <  1:1007703/238‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     AATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGCGACCGCCTATG‑TCA‑CACTCGAACCCTGTAGCCATCATGGTCG‑TACGCCACCGTGCTGTGACGCACTCATCGCCGCT‑GGCGTAGCGCG‑CGTGGTTGCCTCGATG‑CAAGATCCTAACCCGCAGGTCG‑CTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCACC  >  NC_000913/433417‑433886

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