Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,279,031	A→G	I185V (ATT→GTT)	ghxP →	guanine/hypoxanthine transporter GhxP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,279,031	0	A	G	97.8%	156.6 / ‑3.7	46	I185V (ATT→GTT)	ghxP	guanine/hypoxanthine transporter GhxP
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (25/19);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

GCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGT  >  NC_000913/4278784‑4279272                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  gCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTa                                                                                                                                                                                                                                                <  1:213419/251‑1 (MQ=255)    aCATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGt                                                                                                                                                                                                                                             <  2:8440/251‑1 (MQ=255)      catatAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTc                                                                                                                                                                                                                                            >  2:88817/5‑250 (MQ=255)        atTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCAc                                                                                                                                                                                                                                          <  2:326602/250‑1 (MQ=255)           aGCGTACCTGTCGCGCGGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTgg                                                                                                                                                                                                                                       <  2:76434/250‑1 (MQ=255)                           tGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacg                                                                                                                                                                                                                       >  1:93607‑M1/1‑240 (MQ=255)                            gggTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCTCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgac                                                                                                                                                                                                                     <  1:484704‑M1/251‑13 (MQ=255)                              gttccgtgTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCCCGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaa                                                                                                                                                                                                                   <  1:441843‑M1/248‑15 (MQ=255)                                               ggTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaa                                                                                                                                                                                                          <  2:344026‑M1/243‑24 (MQ=255)                                               ggTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaa                                                                                                                                                                                                          >  1:344026‑M1/1‑220 (MQ=255)                                                             tgtAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtg                                                                                                                                                                                    >  2:233060‑M1/2‑206 (MQ=255)                                                                aaTTTCTGCAACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAATCCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaa                                                                                                                                                                                  <  1:120790‑M1/250‑48 (MQ=255)                                                                 aTTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatc                                                                                                                                                                                >  1:270650‑M1/1‑202 (MQ=255)                                                                 aTTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatc                                                                                                                                                                                <  2:111996‑M1/251‑50 (MQ=255)                                                                    tCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgt                                                                                                                                                                              <  2:93607‑M1/250‑52 (MQ=255)                                                                          aCGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccctacccaaaaggtagtaaatcgtttgcttt                                                                                                                                                                       >  2:107561‑M1/1‑193 (MQ=255)                                                                             ggTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgctttta                                                                                                                                                                     >  1:200070‑M1/1‑190 (MQ=255)                                                                                    gtacctggaTTTTGCGCAACTTGCCTCACGGCGTGGCGCACGGGACGGGGCTTGGTATCGGTCTGTTCCTGCTGCTCCTTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCataactaatgacaacaggaaaaaaaaccatccccaaacgggagcgcatcgttggcttttaattggggg                                                                                                                                                             >  2:235146‑M1/5‑183 (MQ=255)                                                                                    gTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactgggg                                                                                                                                                              >  1:98521‑M1/1‑183 (MQ=255)                                                                                      aGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAGCCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggat                                                                                                                                                            <  2:51153‑M1/250‑70 (MQ=255)                                                                                      aGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggatt                                                                                                                                                           <  1:127588‑M1/251‑71 (MQ=255)                                                                                      aGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggatt                                                                                                                                                           >  1:268296‑M1/1‑181 (MQ=255)                                                                                       gCTGGATTTTGCGCAACTTGCCTCACGGTGCGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGGTCATTGCCGATAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCTGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaacaaaccatac                                                                                                                                                                                                 >  1:514659‑M1/1‑180 (MQ=255)                                                                                        cTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGCGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgt                                                                                                                                                         >  1:533741‑M1/1‑179 (MQ=255)                                                                                          ggATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgtt                                                                                                                                                        <  1:133021‑M1/250‑74 (MQ=255)                                                                                                         tGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGTGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcg                                                                                                                                          <  2:70837‑M1/249‑88 (MQ=255)                                                                                                              cACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTGCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgtttt                                                                                                                                   <  2:489495‑M1/251‑95 (MQ=255)                                                                                                               aCGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcctttaactggggattgttcagggccggaatggccacagct                                                                                                                                  >  1:107678‑M1/1‑156 (MQ=255)                                                                                                                           acggcacggGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGAATTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagact                                                                                                                      >  1:557970‑M1/1‑144 (MQ=255)                                                                                                                                        tGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggatagttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggaga                                                                                                          >  2:75944‑M1/1‑131 (MQ=255)                                                                                                                                              cGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggagaatctgt                                                                                                    <  2:233430‑M1/250‑126 (MQ=255)                                                                                                                                                                                      cTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggagaatctgtggcaggaaccgcctctggtatagggggaggcgaagatagc                                                            >  1:110666‑M1/1‑85 (MQ=255)                                                                                                                                                                                                 aaaCCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggagaatctgtggcaggaaacgcccctggtataggggggggcgaaggtagcattatttccgtt                                                >  2:393147‑M1/1‑74 (MQ=255)                                                                                                                                                                                                     ccGCTTGATGTTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGGGTTTTTGTCCCTGGTAGGTCatccctaacgacaacaggaaaaaaaaccatccccaaaaggtagtgactcgttggtttttaactggggattgttcgggtggctttgcggtgttttgaatgttaagcctgtgggagggagaatctggggcaggaaccgcctctggtatagggggaggcgaagatagcattatttcctctccct                                            <  1:144350‑M1/251‑182 (MQ=255)                                                                                                                                                                                                       gCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGACATTATGGgcg                                          >  2:528623/1‑251 (MQ=255)                                                                                                                                                                                                       gCTTGATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggagaatctgtggcaggaaccgcctctggtatagggggaggcgaagatagcattatttcctctccctct                                          >  1:219671‑M1/1‑68 (MQ=255)                                                                                                                                                                                                           gATGGTCTGCCCGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggagaatctgtggcaggaaccgcctctggtatagggggaggggaagatagcattattttctctccctcttctt                                      >  2:457452‑M1/1‑64 (MQ=255)                                                                                                                                                                                                                     ccGTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgataacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggagaatctgtggcaggaaccgcctctggtatagggggaggcgaagatagcattatttcctctccctcttcttcgctgtacc                             >  1:526503‑M1/1‑54 (MQ=255)                                                                                                                                                                                                                       gTTGCGCTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggagaatctgtggcaggaaccgcctctggtatagggggaggcgaagatagcattatttcctctccctcttcttcgctgtacc                             >  1:34099‑M1/1‑52 (MQ=255)                                                                                                                                                                                                                            gcgggttgacttcgcgacctccCCGGTGGTTATGTCACTGGTAGTTCatacctaacgaaaacagaaaaaaaaaccatccccaaaaggtagtgaatcgtttgtttttaacgggggattgttcaggtggcgtgggggtttttggaattttaagactgtgggagggaaaatctgtgtcaggaaccgcctctggtatagggggaggcgaagatagcattatttcctctccctcttcttcgctgtaccagataacc                     <  2:200070‑M1/241‑205 (MQ=255)                                                                                                                                                                                                                            gcTGGGTGACTTCGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggagaatctgtggcaggaaccgcctctggtatagggggaggcgaagatagcgttatttcctctccctcttcttcgctgtaccagataacc                     >  2:126119‑M1/1‑47 (MQ=255)                                                                                                                                                                                                                                   gactccgccaccttcccggTGGTTATGTCACTGGTAGGTCatcccaaacgaaaacagaaaaaaaaaccatccccaaaaggtagtgattcgtttgtttttaactggggtttgttcaggtggcgttggggttttttgaatgttaagacggtgggagggagaatctgtggcaggaaccgcctctggtatagggggaggcgaagatagcattatttcctctccctcttcttcgctgtaccagataaccggcacc               <  2:110666‑M1/241‑211 (MQ=255)                                                                                                                                                                                                                                     cttcGCGACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggagaatctgtggcaggaaccgcctctggtatagggggaggcgaagatagcattatttcctctccctcttcttcgctgtaccagataaccggcaccatt            >  1:174080‑M1/2‑38 (MQ=255)                                                                                                                                                                                                                                         gcgACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggagaatctgtggcaggaaccgcctctggtatagggggaggcgaagatagcattatttcctctccctcttcttcgctgtaccagataaccggcaccattaatt        >  2:398676‑M1/1‑34 (MQ=255)                                                                                                                                                                                                                                          cgACCTTCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggagaatctgtggcaggaaccgcctctggtatagggggaggcgaagatagcattatttcctctccctcttcttcgctgtaccagataaccggcaccattaattt       >  1:154578‑M1/1‑33 (MQ=255)                                                                                                                                                                                                                                                tCCCGGTGGTTATGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatcccaaaaaggtagtaaatcgtttgcttttaactggggatttttcaggtggcgttgcggtgttttaaatgttaagactgtgggagggagaatctgtggcgggaaccgcctctggtatagggggaggcgaagatagcattatttcctctccctcttcttcgctgtaccagataaccggcgccattaatttatagc  <  1:88817‑M1/250‑224 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  GCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGT  >  NC_000913/4278784‑4279272

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