Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,048,787	G→A	F182F (TTC→TTT)	gfcB ←	lipoprotein GfcB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,048,787	0	G	A	97.8%	169.0 / ‑6.0	46	F182F (TTC→TTT)	gfcB	lipoprotein GfcB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (21/24);  minor base C (0/1);  total (21/25)
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.96e-01

GTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCAGCAAAGGAAAAAGTAGAACTCAGGGTCGCAGAGCGAAAATGGTTGTCTTCACTCCAGAGAATATCGCGCGTCCAGCGGCTGCCATCTTTTATTGCCAGCGCGTTACGCAGAGGGTCCTGCCCGGAATTTGTCACTTCCAGC  >  NC_000913/1048541‑1049030                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    gtgtCACCGACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAAtg                                                                                                                                                                                                                                                 <  1:351648/251‑1 (MQ=255)    tCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCGAGTAGCCGAATCGACCCAAAATgtgt                                                                                                                                                                                                                                              >  2:447419/1‑251 (MQ=255)       ccAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGtttt                                                                                                                                                                                                                                           >  1:264638/1‑251 (MQ=255)             ggACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACt                                                                                                                                                                                                                                     <  2:225238/251‑1 (MQ=255)                 aacaTTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGttt                                                                                                                                                                                                                                 <  1:535269/251‑1 (MQ=255)                              gggCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTg                                                                                                                                                                                                                    >  2:124241/1‑251 (MQ=255)                                   aaCGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGt                                                                                                                                                                                                               >  1:387068/1‑251 (MQ=255)                                         aaGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACctc                                                                                                                                                                                                         <  2:223424/251‑1 (MQ=255)                                             gtCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTcc                                                                                                                                                                                                     >  2:308766/1‑251 (MQ=255)                                                       cGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACa                                                                                                                                                                                            <  1:144758/250‑1 (MQ=255)                                                        gCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGCATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACAtt                                                                                                                                                                                          >  1:299710/1‑251 (MQ=255)                                                             gCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTAcaca                                                                                                                                                                                     <  2:11366/251‑1 (MQ=255)                                                                 aTAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGcac                                                                                                                                                                                 <  2:215553/251‑1 (MQ=255)                                                                    aaTAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACAtt                                                                                                                                                                              <  1:103255/251‑1 (MQ=255)                                                                     aTAGTCACCGTTCCTTGCGCAAAGGCATTGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTg                                                                                                                                                                             >  2:419748/1‑251 (MQ=255)                                                                       aGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCg                                                                                                                                                                           <  1:53141/251‑1 (MQ=255)                                                                       aGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCg                                                                                                                                                                           >  1:92308/1‑251 (MQ=255)                                                                                tCCTTGCGCAAAGGCATGGGGTTTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAGGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCGTCCTTCCAGACATTACACAGCACATTGCGACTCGCaa                                                                                                                                                                   <  1:43276/250‑1 (MQ=255)                                                                                tCCTTGCGCAAAGGCATGGGGTGTAATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTCGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACAATGCGACCCGCact                                                                                                                                                                  >  2:115599/1‑249 (MQ=255)                                                                                        cAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCg                                                                                                                                                          <  2:454988/251‑1 (MQ=255)                                                                                                   ggTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTACGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCa                                                                                                                                               <  1:570982/251‑1 (MQ=255)                                                                                                        cATTACGTAAGGTACGCTGGCAATGAAATACGACTGTAATTTATTCATTGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACCTTTGCCGATTTTGACGAACCTGGCCAGTAGCCGAATCTACCCACAATGTGTTTTGCCTCTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACTTTACACAGCACATTGCGTCCCGCAATATTCAGCGTCTCTAATCCAGCaaa                                                                                                                                          <  1:273253/251‑1 (MQ=255)                                                                                                                          ggCAATGAAATACGACTGTAATTTATTCATGTTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGTCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCATCAAAGGAAAAAGTAGAACTc                                                                                                                          <  2:101085/249‑1 (MQ=255)                                                                                                                              aTGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCAGCAAAGGAAAAAGTAGAACTCAGGGTc                                                                                                                    <  2:227295/251‑1 (MQ=255)                                                                                                                                  aaTACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCAGCAAAGGAAAAAGTAGAACTCAGGGTCGCag                                                                                                                <  2:169947/251‑1 (MQ=255)                                                                                                                                    tACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCAGCAAAGGAAAAAGTAGAACTCAGGGTCGCag                                                                                                                >  1:500024/1‑249 (MQ=255)                                                                                                                                     aCGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCAGCAAAGGAAAAAGTAGAACTCAGGGTCGCagag                                                                                                              <  1:243839/250‑1 (MQ=255)                                                                                                                                      cGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATGTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTATACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCAGCAAAGGAAAAAGTAGAACTCAGGGTCGCAGAGCg                                                                                                            <  2:276667/251‑1 (MQ=255)                                                                                                                                          tGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCAGCAAAGGAAAAAGTAGAACTCAGGGTCGCAGAGCGaaa                                                                                                         >  1:575920/1‑250 (MQ=255)                                                                                                                                                   ttCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCAGCAAAGGAAAAAGTAGAACTCAGGGTCGCAGAGCGAAAATGGTTGTCt                                                                                               >  1:197923/1‑251 (MQ=255)                                                                                                                                                           gCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCAGCAAAGGAAAAAGTAGAACTCAGGGTCGCAGAGCGAAAATGGTTGTCTTCACTCCa                                                                                       <  1:173942/251‑1 (MQ=255)                                                                                                                                                                    aGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGATCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCAGCAAAGGAAAAAGTAGAACTCAGGGTCGCAGAGCGAAAATGGTTGTCTTCACTCCAGAGAATATc                                                                              <  1:523756/251‑1 (MQ=255)                                                                                                                                                                         cGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCAGCAAAGGAAAAAGTAGAACTCAGGGTCGCAGAGCGAAAATGGTTGTCTTCACTCCAGAGAATATcgcgcg                                                                         <  2:52782/251‑1 (MQ=255)                                                                                                                                                                          gTCATTTCTACGGGAATAACCCCTGCGCCTCACATTTACCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCAGCAAAGGAAAAAGTAGAACTCAGGGTCGCAGAGCGAAAATGGTTGTCTTCACTCCAGAGAATATCGCGCGt                                                                        >  1:271475/1‑251 (MQ=255)                                                                                                                                                                              tttCTACGGGAATAACCACTGAGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACTCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCAAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCAGCAAAGGAAAAAGAAGAACTAAGGGTCGCAGAGCGAAAATGGTTGTCTTCACTCCAGAGAATATcgcgcgacc                                                                      >  1:453482/1‑246 (MQ=255)                                                                                                                                                                                tCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCACTATTCAGCGTCTCTAATCCAGCAAAGGAAAAAGTAGAACTCAGGGTCGCAGAGCGAAAATGGTTGTCTTCACTCCAGAGAATATCGCGCGTCCAGCg                                                                  >  1:480707/1‑251 (MQ=255)                                                                                                                                                                                    cGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCAGCAAAGGAAAAAGTAGAACTCAGGGTCGCAGAGCGAAAATGGTTGTCTTCACTCCAGAGAATATCGCGCGTCCAGCGGCTg                                                              >  1:260181/1‑251 (MQ=255)                                                                                                                                                                                                     ccTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCAGCAAAGGAAAAAGTAGAACTCAGGGTCGCAGAGCGAAAATGGTTGTCTTCACTCCAGAGAATATCGCGCGTCCAGCGGCTGCCATCTTTTATTGcc                                               >  1:527553/1‑249 (MQ=255)                                                                                                                                                                                                          cATTTGCCGACTTTGACGAACCTGGCCCGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGCAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCAGCACAGGAAAAAGTAGAACTCAGGGTCGCAGAGCGAACATGGTTGTCTTAACTCCAGAGAATATCGCGCGTCCAGCGGCTGCCATCTTTTATTGCCAGCGCGt                                        >  2:584228/1‑251 (MQ=255)                                                                                                                                                                                                                           gAACCTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCAGCAAAGGAAAAAGTAGAACTCAGGGTCGCAGAGCGAAAATGGTTGTCTTCACTCCAGAGAATATCGCGCGTCCAGCGGCTGCCATCTTTTATTGCCAGCGCGTTACGCAGAGGGTCCTGc                       >  2:64611/1‑251 (MQ=255)                                                                                                                                                                                                                               cgggCCAGGAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCAGCAAAGGAAAAAGTAGAACTCAGGGTCGCAGAGCGAAAATGGTTGTCTTCACTCCAGAGAATATCGCGCGTCCAGCGGCTGCCATCTTTTATTGCCAGCGCGTTACGCAGAGGGTCCTGCCCgg                   <  2:439421/249‑1 (MQ=255)                                                                                                                                                                                                                               cTGGCCAGTAGCCGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCAGCAAAGGAAAAAGTAGAACTCAGGGTCGCAGAGCGAAAATGGTTGTCTTCACTCCAGAGAATATCGCGCGTCCAGCGGCTGCCATCTTTTATTGCCAGCGCGTTACGCAGAGGGTCCTGCCCgg                   >  1:93311/1‑251 (MQ=255)                                                                                                                                                                                                                                          ccGAATCGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCAGCAAAGGAAAAAGTAGAACTCAGGGTCGCAGAGCGAAAATGGTTGTCTTCACTCCAGAGAATATCGCGCGTCCAGCGGCTGCCATCTTTTATTGCCAGCGCGTTACGCAGAGGGTCCTGCCCGGAATTTGTCACt        <  1:358701/251‑1 (MQ=255)                                                                                                                                                                                                                                            gAATCGGCCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCAGCAAAGGAAAAAGTAGAACTCAGGGTCGCAGAGCGAAAATGGTTGTCTTCACTCCAGAGAATATCGCGCGTCCAGCGGCTGCCATCTTTTATTGCCAGCGCGTTACGCAGAGGGTCCTGCCCGGAATTTGTCACTTc      <  1:9898/251‑1 (MQ=255)                                                                                                                                                                                                                                                cGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCAGCAAAGGAAAAAGTAGAACTCAGGGTCGCAGAGCGAAAATGGTTGTCTTCACTCCAGAGAATATCGCGCGTCCAGCGGCTGCCATCTTTTATTGCCAGCGCGTTACGCAGAGGGTCCTGCCCGGAATTTGTCACTTCcagc  <  1:425773/251‑1 (MQ=255)                                                                                                                                                                                                                                                cGACCCAAAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCAGCAAAGGAAAAAGTAGAACTCAGGGTCGCAGAGCGAAAATGGTTGTCTTCACTCCAGAGAATATCGCGCGTCCAGCGGCTGCCATCTTTTATTGCCAGCGCGTTACGCAGAGGGTCCTGCCCGGAATTTGTCACTTCcagc  >  2:148456/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    GTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCAGCAAAGGAAAAAGTAGAACTCAGGGTCGCAGAGCGAAAATGGTTGTCTTCACTCCAGAGAATATCGCGCGTCCAGCGGCTGCCATCTTTTATTGCCAGCGCGTTACGCAGAGGGTCCTGCCCGGAATTTGTCACTTCCAGC  >  NC_000913/1048541‑1049030

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