Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,639,053	T→G	T292T (ACT→ACG)	creD →	putative inner membrane protein CreD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,639,053	0	T	G	93.0%	113.4 / ‑0.5	43	T292T (ACT→ACG)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base G (19/21);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.15e-01

GTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGT  >  NC_000913/4638808‑4639289                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             gtggtgCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGAc                                                                                                                                                                                                                                           <  1:420908/249‑1 (MQ=255)          gccgggcgTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTa                                                                                                                                                                                                                                 <  1:342657/248‑1 (MQ=255)                    acgtgttctcttccggtctcAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGGCGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCAtt                                                                                                                                                                                                                      <  1:443572/232‑1 (MQ=255)                          aaaTGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGAAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATAACGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGAGACGACGCCAGCCGATCAATACCCATTAACGGACCGGGCGAATAAGTACGCCATTctgctg                                                                                                                                                                                                                >  1:622077/1‑251 (MQ=255)                                   tagccAGAAACTGGCCTCATCCCATTTTTTTAGGTGATTTTCTACCATCCAAACGGGAAGGTATCTAGTCAGGGTGTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAAAAATCGCGGTGGGCGTTTGGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCTTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACt                                                                                                                                                                                                        <  2:649176/247‑1 (MQ=255)                                       cAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTCCCACCCAACCGGGAGGTTAGCGGGTCAGGTTTTCAGGCGCATTGCCAAAGCAGCTGGTTTGATAATAATCTCGGTGAGCGTTTTGCTTCAGGCAAGGATACCGGCTGGGAAAACTTCCCCGCCTTTAGCGCCCCAGGGCCAACGCCACCCGATCAATACCAAATAAGGGAACGGGCGAATAAGTACGACAtataacaggaagagcgccggg                                                                                                                                                                                                   >  2:443572/1‑231 (MQ=255)                                        aGCATCTGGCCGCAGCCCAGTTTTTGAGGTGATTTGTTGCGATCCAAACGGGAAGTGAGCGTGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGGTCGTTTTGCTTCTGGCAATGATACCGGCTGGGAAAACGTCCCCGGGTTTAGCGGCGCAGTGTGGACGTCAGCCGATCAACACCAATTAACGGTCCGGGCGCCTAAGGACGCCATTCTGCTGATTGCACTGACttt                                                                                                                                                                                                  <  1:415290/251‑1 (MQ=255)                                                       cccAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTg                                                                                                                                                                                    <  1:477872/250‑1 (MQ=255)                                                          aGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCCGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTTCTGATTGCCCTGACTTTTTTGGCGTTCTTTGtttt                                                                                                                                                                                >  2:220162/1‑251 (MQ=255)                                                                      gATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCa                                                                                                                                                                     <  1:219591/250‑1 (MQ=255)                                                                               ccAGCCAAACGGGACGTTAGCGAGTCAGGTTTTCCGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGGCGCCGCCAGCCGATCAATCCCCATTCACGGCCAGGGCGCCTAAGAACGCCATTCTGCTGATTGCACTTAATTTTATGGCGTTTTTTGttttttacaagctcacagcgacaac                                                                                                                                                           >  2:614152/1‑231 (MQ=255)                                                                                 aGCCAGACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGtt                                                                                                                                                         >  1:86822/1‑251 (MQ=255)                                                                                 aGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGtt                                                                                                                                                         >  1:203731/1‑251 (MQ=255)                                                                                  gCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGttt                                                                                                                                                        <  1:351637/251‑1 (MQ=255)                                                                                         gggAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGCTCAATACCAATTAACGGACCCGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGACAAGCTCACCGCGCAACGTTTACAccc                                                                                                                                                  >  2:471254/1‑250 (MQ=255)                                                                                                  gCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCCatat                                                                                                                                        >  1:400170/1‑251 (MQ=255)                                                                                                         aGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCtgg                                                                                                                                 <  1:220162/251‑1 (MQ=255)                                                                                                           gTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCCATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGCCtt                                                                                                                                                                                                   >  2:108677/1‑183 (MQ=255)                                                                                                           gTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACtt                                                                                                                                                                                                   <  1:108677/183‑1 (MQ=255)                                                                                                                      caatGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTGCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCAtt                                                                                                                    <  1:212603/248‑1 (MQ=255)                                                                                                                       aatGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCCGTGACGACGCCAGCCGATCAATACCAATTAAAGGACCGGGCGGCTAAGTACTCCATTCTGCTGCTTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACCACCACCGCGCAACGTTTTAAACCAATGCAATTTTTTATGGTGGGGCttttcttg                                                                                                                   >  2:93008/3‑246 (MQ=255)                                                                                                                        atGGCAAAGAAGCTTGTTTGCTAATAATCGCGGTGCGCGTTTTGCTCTAGGCAATGATACCGGTTGGTAAAACTGCCCCGCTTTTAGCGTCGCAGTGACGCCGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGTTGGGGCTTTCATTgg                                                                                                                  <  2:510358/250‑1 (MQ=255)                                                                                                                          ggCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGGACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTg                                                                                                                <  1:640611/251‑1 (MQ=255)                                                                                                                              aaGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGCAGCCAGCCGATCAATACCAATTAACGGAACGGGCGACTAAATACGCCATTCTGCTGATTGCACTGCCTTTTTTGGCGTTTTTTGTTTTTGAAACGCTCAACGCGCAACGTTTACAACCAATGCAACTTTTTCTGGGGGGGCTTTTCTTGGtgttgt                                                                                                            >  2:474330/1‑251 (MQ=255)                                                                                                                                   gcTGGTTTGCTAATAATCTCGGTGAGCGTTTTGATTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACCCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCCTTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGGTTTTGAAACCCTCACCGCGCAACGTTTACACCCAATGCCATATTTGCTTGTGGGGCTTTCATTGGTGCTGTTTTAt                                                                                                       >  2:41578/1‑251 (MQ=255)                                                                                                                                      ggTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTAtt                                                                                                      <  1:99315/249‑1 (MQ=255)                                                                                                                                             taataaTCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCACCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCAATTCTGCTGAATGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCCTTGGTGATGTTTTATTAGATCTTGc                                                                                             >  1:374295/1‑251 (MQ=255)                                                                                                                                              aataatCTCGGTGATGGTTTTGCTTCAGGCAATGATACCGGCTGGTAAAACTTCCCCGCTTTTAGCGTCGCTGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATAGTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCt                                                                                            <  1:327000/251‑1 (MQ=255)                                                                                                                                                      cGGTTAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAAc                                                                                                                                                            >  1:115832/1‑179 (MQ=255)                                                                                                                                                      cGGTTAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTACCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAAc                                                                                                                                                            <  2:115832/179‑1 (MQ=255)                                                                                                                                                      cGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCttt                                                                                    <  1:471254/251‑1 (MQ=255)                                                                                                                                                        gTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCt                                                                                  >  1:76486/1‑251 (MQ=255)                                                                                                                                                                    cTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCCACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGt                                                                      >  1:641437/1‑251 (MQ=255)                                                                                                                                                                                aTACCGTCTGGGAAACCTTCCCCGCTTTTGGCTGCGCAGTGACGACTCCAGCCGATCAATACCAATTAACGGACCGGGCTACTAAGTACGCCATTCTGCTGATTGCACTGACTGTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTGACACCCAATGCAATATTTGCTGGTGGGGCTTTCAGTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCa                                                          <  2:76486/251‑1 (MQ=255)                                                                                                                                                                                    cGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGa                                                      <  1:614152/251‑1 (MQ=255)                                                                                                                                                                                         gggAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATc                                                 <  1:430129/251‑1 (MQ=255)                                                                                                                                                                                          ggAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCg                                                >  1:638004/1‑251 (MQ=255)                                                                                                                                                                                                                gTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAAt                          <  1:93008/251‑1 (MQ=255)                                                                                                                                                                                                                        cacgagcccAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGta                  <  2:638004/244‑1 (MQ=255)                                                                                                                                                                                                                                cAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTg                                                                                                                                     >  1:140233/1‑128 (MQ=255)                                                                                                                                                                                                                                cAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTg                                                                                                                                     <  2:140233/128‑1 (MQ=255)                                                                                                                                                                                                                                cAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTACACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTACTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCGTTCTGATTGGAGCGTTAATGAACGGAATTTATTTg          >  1:379025/1‑251 (MQ=255)                                                                                                                                                                                                                                        cAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGt  >  1:12310/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             GTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGT  >  NC_000913/4638808‑4639289

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