Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,425,696	C→A	N59K (AAC→AAA)	priB →	primosomal replication protein N

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,425,696	0	C	A	95.7%	134.5 / ‑5.0	47	N59K (AAC→AAA)	priB	primosomal replication protein N
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (20/25);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCGCGATGATTTCGCAAACGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGACACGAAAACCAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGC  >  NC_000913/4425451‑4425946                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                            cGTCGCGTTGATTTCGCAAACGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTTTCAAATGCCCGGTATTGTTAGCGGAAACGAAAAACAGGc                                                                                                                                                                                                                                                        >  2:37201/1‑251 (MQ=255)  gTCGCGATGATTTCGCAAACGACACCGATGATGATGCTGACGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCCCCGTGTGCCGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGCGAATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGCCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGCAAAACAGGcc                                                                                                                                                                                                                                                       >  1:349959/1‑251 (MQ=255)    cgcgATGATTTCGCAAACGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCa                                                                                                                                                                                                                                                      <  1:73688/250‑1 (MQ=255)          gatgtCGCCAACGAAACCGCTGATGATGCTGAAGCTGGGGATGCTGAAGAGTAATTTCTGATGACCAACCGTCGGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCATTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTAcgca                                                                                                                                                                                                                                               <  2:130437/247‑4 (MQ=255)               cGCAAACGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGt                                                                                                                                                                                                                                           >  1:111358/1‑250 (MQ=255)                   aaCGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGACTTCCTCACTGCCCGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGGGTCAAATGCCCGTTATTGGTAGCGGAAACGAAAACCAAGCCATTACTCACAGGATaa                                                                                                                                                                                                                                       >  2:404360/1‑250 (MQ=255)                         aCCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCgg                                                                                                                                                                                                                                <  2:417292/251‑1 (MQ=255)                                  gatgCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCAt                                                                                                                                                                                                                       <  1:358360/251‑1 (MQ=255)                                                 gATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTTCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCTGCAGTCGCATAACCGTTCAGGGGtt                                                                                                                                                                                                        <  2:75312/251‑1 (MQ=255)                                                          gagTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATg                                                                                                                                                                                              <  1:279868/251‑1 (MQ=255)                                                                tttCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAAAGGTCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATGcccca                                                                                                                                                                                         >  2:195514/1‑247 (MQ=255)                                                                  tCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGGTATTGTTAGCGGAAACGAAAAACAGGCCATTAATAACAGTATAACGGTCGGCAGTCACATAACCGTTCAGGGGTT‑CATTTCATGcctaaagg                                                                                                                                                                                      >  1:218308/1‑245 (MQ=255)                                                                     gatgaCCAACCGCTTGGTGTGGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCTCCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATGCCACAAGGCaa                                                                                                                                                                                   <  2:258398/251‑1 (MQ=255)                                                                      atgaCCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGACGGCAGTCGCATAACCGTTCCGGGGTT‑CATTTCATGCCACAAGGcca                                                                                                                                                                                   >  2:2847/1‑248 (MQ=255)                                                                         cctcaccgccgggTGTTGTCCGGCACTGTGGGCAGGGCCCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGGCAGTTCGTGGTTGAGCATCGTTTTGTGCAGGAGGAAGCCGTATTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGAGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGCCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATGCCACAAGGCgaagaa                                                                                                                                                                               <  2:250492/240‑1 (MQ=255)                                                                               cGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAAACGTTCAGGGGTT‑CATTTCATGCCACACAGCAAAGAACGGACt                                                                                                                                                                         >  2:91063/1‑251 (MQ=255)                                                                                       gttgtCCGGCACCGTGTGCAGGGCTCCCCTGCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGACAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGGTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGACTGAGCaaa                                                                                                                                                                  <  2:333697/250‑1 (MQ=255)                                                                                               gCACCGTGTGCAGGGCTCCCCTTCGAACAGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTCTTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGACGCATAACCGTTCAGGGGTTACATT‑CCTGCCACAAGGCCACGAACGGACTGAgacaaatacatcgg                                                                                                                                                         >  2:207292/1‑238 (MQ=255)                                                                                               gCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGATTTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAAt                                                                                                                                                                <  1:207292/244‑1 (MQ=255)                                                                                                                  ccTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGGCAAATGCCCGTTATTGTTAGCGGAAACGCAAAACAGGCCCTTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGACTGTGCAAAATGGTTTTGCATTCCGAGCAGataaaat                                                                                                                                      >  2:47480/1‑246 (MQ=255)                                                                                                                        aaaGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAg                                                                                                                                <  2:111358/251‑1 (MQ=255)                                                                                                                         aaGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGAt                                                                                                                                  >  1:104090/1‑248 (MQ=255)                                                                                                                             tCAGTCCATCAGGAATTCCTCACTGCCAGTTAGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCt                                                                                                                           >  1:107128/1‑251 (MQ=255)                                                                                                                                    aTCAGGATTTCCACCCTGCCATTTCGTGCTTGCGCATCGTTCGGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGGGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACGCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACt                                                                                                                    <  2:206593/251‑1 (MQ=255)                                                                                                                                              ccTCCCTGCCAGTTCGTGCTGGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATgg                                                                                                          <  1:342090/251‑1 (MQ=255)                                                                                                                                              ccTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGTTCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATg                                                                                                           <  1:28978/250‑1 (MQ=255)                                                                                                                                                  aCTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATGCCACAAGGCaaa                                                                                                                                                                                  >  2:197970/1‑175 (MQ=255)                                                                                                                                                        aGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATGCCCCACGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCCCGtttttt                                                                                                >  2:63593/1‑251 (MQ=255)                                                                                                                                                        aGTTCGTGCTTGAGCATCGTTCTGGGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTAttt                                                                                                <  2:22060/251‑1 (MQ=255)                                                                                                                                                                       aTCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGGTAGATTCTGGAGACTAGCCATATGGCACGGTATTTCCGtcgtcgacattt                                                                                 >  2:371766/1‑245 (MQ=255)                                                                                                                                                                       aTCGTTCGGTGCAGGTGGAAGCCGGGTTTAACCGGCAGGCGTGGTGTCAAAGGCCCGTTATTGTGAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAGCCGTTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGtt                                                                                 <  2:147368/251‑1 (MQ=255)                                                                                                                                                                        tctttCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCgtcgtc                                                                                        <  2:366953/240‑1 (MQ=255)                                                                                                                                                                        tCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCgtcgtc                                                                                        >  1:366953/1‑243 (MQ=255)                                                                                                                                                                               gtgCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGtt                                                                         <  1:195514/251‑1 (MQ=255)                                                                                                                                                                               gggCAGGAGGAAGCCGGCTTCCACCGGCAGGCGTGGTGGCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCTTCGTCGCAAGTTCTGCCGt                                                                          <  2:38556/248‑1 (MQ=255)                                                                                                                                                                                            ccGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAgg                                                            <  1:91063/251‑1 (MQ=255)                                                                                                                                                                                             cGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGCCTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCCCCGCGGAAGGc                                                           >  2:396179/1‑251 (MQ=255)                                                                                                                                                                                                 tttCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTACCGGAAACGAAAAACAGGCCATTACTCACCGTATAACGGTCGGCAGTCGCATAACCGGTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGACTGCGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTc                                                       >  2:416805/1‑251 (MQ=255)                                                                                                                                                                                                                    gtCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGGTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCAGATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGa                                    <  1:37201/251‑1 (MQ=255)                                                                                                                                                                                                                      cAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGata                                  <  1:396179/251‑1 (MQ=255)                                                                                                                                                                                                                         atgaccGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCg                               <  1:348323/247‑1 (MQ=255)                                                                                                                                                                                                                           gccctctATTTTTAGCGGAAACGAAAAACAGGCCATTTCTCACAGTATAACGGCCTGCATTCGCATAACCGTTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGGCGCAATTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAGCGACTATAAAGATATCGCt                             <  1:365291/245‑1 (MQ=255)                                                                                                                                                                                                                               gTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCTGTCGTCGCAAGTTCTGCCGTTTACCCGCGGCATGCGTTCAAGAGATCGACTATAAAGATATctctccgc                         >  2:13356/1‑244 (MQ=255)                                                                                                                                                                                                                               gTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCACGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGCGATCGACTATAAAGCTATCGCTACGc                         >  2:202748/1‑251 (MQ=255)                                                                                                                                                                                                                                          cGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTAc              <  1:54756/251‑1 (MQ=255)                                                                                                                                                                                                                                                 gAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCg       >  1:184970/1‑251 (MQ=255)                                                                                                                                                                                                                                                    aaCCAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGaa     >  1:54238/1‑250 (MQ=255)                                                                                                                                                                                                                                                      ccAGGCCATTACTCACAGTATAACGGTCGGCAGTGGCATAACCGTTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAAAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGc  >  1:342193/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                            CGTCGCGATGATTTCGCAAACGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGACACGAAAACCAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTT‑CATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGC  >  NC_000913/4425451‑4425946

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: ‑