Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,627,084	G→A	P391P (CCG→CCA)	radA →	DNA recombination protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,627,084	0	G	A	97.9%	149.3 / NA	47	P391P (CCG→CCA)	radA	DNA recombination protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (19/27);  minor base C (1/0);  total (20/27)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGGGAGATATGTCGTCATTT  >  NC_000913/4626838‑4627326                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   gttgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTATTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAgg                                                                                                                                                                                                                                                >  1:152773/3‑251 (MQ=255)     ctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGAGGTCTGCAAATGGCCGATCAGGATGGGTTTGTGACCGTGTTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGAGTTACTGCTGGCGGTGGTTTCAGGCCTGCGGGACAGACCGCTGCCCCAGGGTCt                                                                                                                                                                                                                                            >  2:42788/1‑251 (MQ=255)           gagaTTCAGGCGCTGGTCGATCACTCGATGATGGCGAAGCCGCGCCGCGTGGCAGTGGGGCGGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTggtggt                                                                                                                                                                                                                                      <  2:200690/251‑1 (MQ=255)                         ggTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTc                                                                                                                                                                                                                         <  1:259074/250‑1 (MQ=255)                                                    gctgcgcgTGGCTGTGTGGCTGGTACAAACCCTTCTGGCAATCGTGCTGTCTGTCGTGCGCCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTATCCGAATCCAGCGCCGACTTAGCGTTACTGCTGGCGATGTTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGACCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGAAGATCCGCCCGGTGc                                                                                                                                                                                             <  2:150488/247‑1 (MQ=255)                                                       ccgcGTGGCGGTGGGGCTGGAACAAAACCGGCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCa                                                                                                                                                                                          <  2:64895/251‑1 (MQ=255)                                                           gTGGCAGTGGGGCTGGAACCAAACCGTCCGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGCTCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGCAACCAGCGCCGCCTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGGGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCCGCgg                                                                                                                                                                                      >  2:172839/1‑251 (MQ=255)                                                              gCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGACGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTc                                                                                                                                                                                    >  1:117287/1‑250 (MQ=255)                                                                     ggCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCCCGGTGGTCTGCAAATGGCCGCTCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGCCCGCTGCCACAGGATCTGGGGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCAGGGGCCCGGAGGTCAGGaccg                                                                                                                                                                             >  2:262005/1‑247 (MQ=255)                                                                            aCAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAtttc                                                                                                                                                                       >  1:166032/1‑246 (MQ=255)                                                                              aaaaCCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAtttcaga                                                                                                                                                                    >  2:221967/1‑244 (MQ=255)                                                                                 aCCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAgcg                                                                                                                                                                >  1:10762/1‑251 (MQ=255)                                                                                   cGTCTGGCAATCCTGCTGGCTGGCTTGCACCGTCACTGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAgcggc                                                                                                                                                              <  2:82455/251‑1 (MQ=255)                                                                                      cTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGa                                                                                                                                                            >  1:115028/1‑250 (MQ=255)                                                                                             tCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGt                                                                                                                                                    >  1:131679/1‑251 (MQ=255)                                                                                                 gctggctgTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTc                                                                                                                                                <  1:218039/251‑1 (MQ=255)                                                                                                           ttGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGAt                                                                                                                                      >  1:216232/1‑251 (MQ=255)                                                                                                             gCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGAtt                                                                                                                                     <  1:96571/250‑1 (MQ=255)                                                                                                                                cAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTGCCGGCGGCTAACGtacc                                                                                                                 <  1:221967/251‑4 (MQ=255)                                                                                                                                 aaaTGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCg                                                                                                                <  2:217413/251‑1 (MQ=255)                                                                                                                                            cAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGC‑ACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGAGCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGGTTGTTCCGGCGGCTAACGTTCCGAAAAAAGCGCCg                                                                                                    <  1:76894/251‑1 (MQ=255)                                                                                                                                             aGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCg                                                                                                    <  1:244410/251‑1 (MQ=255)                                                                                                                                                     tttGTAAAAGTGGTCGGCGGCGTGCAGGTAACCGCAAACAGAGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGACACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTACAGAAGCGGCGAAACACGGCTTTAGCCGGGAGATAGTTCCGGCGACTAACGTAACGAAAAAAgacacggaaggggt                                                                                            >  2:186850/1‑238 (MQ=255)                                                                                                                                                            ccGTGGTCGGCGGCGTGAAGGTAACCGAACCCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGAtt                                                                                     <  2:166032/250‑1 (MQ=255)                                                                                                                                                               tGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGAttttt                                                                                  >  1:79245/1‑251 (MQ=255)                                                                                                                                                                gggcggcggCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTGACTGCTGGCGATGGTTTCCAGCCTGCGGGACAGACCGCTGCCACAGGATCTGGTGGTGGTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGTCCGGGGGCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTGCGGCGTCTAACGGACCTAAAAAAGCGCCGGAAGGGATGCAGAttgtt                                                                                  <  2:243160/247‑4 (MQ=255)                                                                                                                                                                  gcggcggCTTGAATGTAACCGAATCCTGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTTACAGACCTCTGCCACAGGATCTGGTGGTGTGTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCATATTTTTGGc                                                                               <  1:96942/250‑1 (MQ=255)                                                                                                                                                                    ggcggcGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGt                                                                             >  1:57114/1‑251 (MQ=255)                                                                                                                                                                    ggcggcGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGt                                                                             >  2:165583/1‑251 (MQ=255)                                                                                                                                                                          gTGAGGGTAACCGAATCCAGCGCCGACTTAGCGTTACTGCTGGCGAGGGTTTCCAGGCTGCGTGACGGACCGCTGCCACAGGAGCTGGTGGTGTTTGGTGAAGTCGTGTTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAGTTCAGAAGCGGCGAAACACGGTTTTCGCGGGGGGATTGTTCCGGCGGCGAACGTACCGAAAAAATCGCCGGAAGGGATGCAGAGTTGTGGCGGTaaaaa                                                                       <  2:34466/251‑1 (MQ=255)                                                                                                                                                                               ggTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTAGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTAAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACAGAAAAAAGCGCCGGAAGGGATGCAGAATTTTGGAGTTAAAAAActct                                                                  >  2:72463/1‑251 (MQ=255)                                                                                                                                                                                  aaCCGAAACCCGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGACACAGGATCGGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGCTCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCg                                                               <  2:152773/251‑1 (MQ=255)                                                                                                                                                                                  aaCCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGTTCAGGAACGAATTTCAGAATCGGCGAAACACGGTTTGCGCCGGGGGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCg                                                               <  1:242046/251‑1 (MQ=255)                                                                                                                                                                                     cGAAACCAGCGCCGGCTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAATCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGAcg                                                            <  1:228019/251‑1 (MQ=255)                                                                                                                                                                                     cGAAACAAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAACGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGAcg                                                            >  2:235801/1‑251 (MQ=255)                                                                                                                                                                                              cgCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCg                                                   <  2:115028/250‑1 (MQ=255)                                                                                                                                                                                                 cGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCgtgt                                                <  1:234857/251‑1 (MQ=255)                                                                                                                                                                                                 cGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCgtgt                                                >  1:91523/1‑251 (MQ=255)                                                                                                                                                                                                 cGACTGAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCgtgt                                                <  1:172839/251‑1 (MQ=255)                                                                                                                                                                                                  gACTTAGCGTGACTGCTGGCGAGGGTTTCCAGCCTGCGTGACAGTCCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGtt                                               <  2:76199/251‑1 (MQ=255)                                                                                                                                                                                                      tAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTcg                                             <  1:128559/249‑1 (MQ=255)                                                                                                                                                                                                       aGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGCCCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGTCTGGCAGGGTAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTcgacg                                          <  1:120068/251‑1 (MQ=255)                                                                                                                                                                                                           ttACTGCTGGCGAGGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCgacgac                                        <  2:10762/249‑1 (MQ=255)                                                                                                                                                                                                                            ttCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGGGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGTGATGCATATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCTTGTTCGACGACATATAAttgggaatacgg                      <  2:183409/250‑12 (MQ=255)                                                                                                                                                                                                                                  ccTGCGTGACAGACCTGTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGgaga               <  1:103437/251‑1 (MQ=255)                                                                                                                                                                                                                                    tGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGat              >  2:29019/1‑250 (MQ=255)                                                                                                                                                                                                                                               ccGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCCGGGGAGATCCGCCCGGTGCCCCGCGGTAAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCAGGCGGCTAACGTACCGACAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTAGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTAAttt  >  2:120269/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   ACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGGGAGATATGTCGTCATTT  >  NC_000913/4626838‑4627326

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