Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,623,838	G→A	Q265Q (CAG→CAA)	ppk →	polyphosphate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,623,838	0	G	A	97.7%	130.4 / NA	43	Q265Q (CAG→CAA)	ppk	polyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (24/18);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATAAAGCCCAGTTCCGCAATGGTTTTGATGCCATTCGCGAACGCGATGTGTTGCTC  >  NC_000913/2623598‑2624072                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                           aCCGCCAGAGGCCCCGCGTCGACTCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACGGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGc                                                                                                                                                                                                                                  <  2:139610/251‑1 (MQ=255)       aGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGCGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATg                                                                                                                                                                                                                            >  2:64334/1‑251 (MQ=255)        gAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATGCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATgc                                                                                                                                                                                                                           <  2:288117/251‑1 (MQ=255)             gccgcGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTgg                                                                                                                                                                                                                      >  1:505276/1‑251 (MQ=255)                                           ggATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTAttt                                                                                                                                                                                        >  2:534150/1‑251 (MQ=255)                                                                   tgatgaTATTTGCAAAGGCTTTTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCGTGATGGAGTTGGTGTCTGCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGCCTACTCCCCGCTACGACTCCATCGTccac                                                                                                                                                                 <  2:310528/250‑3 (MQ=255)                                                                       gatATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCcggcgg                                                                                                                                                            >  1:499781/1‑251 (MQ=255)                                                                              tCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTAt                                                                                                                                                     >  1:329275/1‑251 (MQ=255)                                                                                   ggCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATa                                                                                                                                                 >  2:177772/1‑250 (MQ=255)                                                                                    gCTTCTTTGATTATGACGCGCTGAAGGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAAt                                                                                                                                               <  1:413452/251‑1 (MQ=255)                                                                                           tGATTATGACGCGCTGAATGCCTATTCACTGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATACTTTTaaa                                                                                                                                         >  2:525675/1‑250 (MQ=255)                                                                                               tATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGGCGTTATCATAATTTTAAAGACt                                                                                                                                     <  2:434250/250‑1 (MQ=255)                                                                                                  gACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGTTTTTGTGCATGAGATGGAAGACAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTTAAGTGGTACGCGAAAACCTGACTATTTCCCGCTACGACTCAATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTtta                                                                                                                                  <  1:313465/250‑1 (MQ=255)                                                                                                     gcgcTGAATGCCTATTCCATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTttatta                                                                                                                               >  1:119043/1‑250 (MQ=255)                                                                                                     gcgcTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTGAAGACTTTATTaa                                                                                                                              <  2:329275/251‑1 (MQ=255)                                                                                                          gAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTcc                                                                                                                         >  2:281683/1‑251 (MQ=255)                                                                                                           aaTGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTccc                                                                                                                        >  2:143598/1‑251 (MQ=255)                                                                                                            aTGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTcccc                                                                                                                       <  2:340479/251‑1 (MQ=255)                                                                                                                      aaTGAAGATGACCCGCGCTGCCGAATACGATTTCGTGCATGAGCTGGCAGCCAGCCTGATGGAGTTGATGTCTTACAGTCTCAAGCAGCGTTAAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGGTACGCGAAAAACTGACTATTTCCCGCTCCGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTCATTTCCCCAATGTCGGc                                                                                                              >  2:305410/1‑250 (MQ=255)                                                                                                                               gccccGCGATGCCTAATTCGATTTAGTGCATGAGATGGAAGCCAGCCGGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCTCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAAt                                                                                                     <  2:50045/248‑1 (MQ=255)                                                                                                                               gACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCAGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATc                                                                                                    >  1:317353/1‑251 (MQ=255)                                                                                                                                aCCCGGGATGCCGAATCCGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTACCGTATCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCCATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCCTAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCt                                                                                                   >  2:289979/1‑251 (MQ=255)                                                                                                                                       aTGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAAc                                                                                            >  1:26307/1‑251 (MQ=255)                                                                                                                                        tGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAAACTGGTGAAc                                                                                            >  1:31785/1‑250 (MQ=255)                                                                                                                                            gAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTCCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAAcc                                                                                       >  2:204574/1‑251 (MQ=255)                                                                                                                                             aaTACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAAcc                                                                                       >  1:25021/1‑250 (MQ=255)                                                                                                                                                      tAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGATGAGCCGGTGAGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACCAACCACTGCCGCGt                                                                             >  1:287761/1‑251 (MQ=255)                                                                                                                                                        gTGCATGAGATGGAAGCCAGCCTGATTGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAATCCACTGCCGCGttt                                                                           <  1:166677/251‑1 (MQ=255)                                                                                                                                                        gTGCATGAGAGGGAAGCCAGCCGGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTATCTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGGGTTACGCTAAAAACTTGCTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGttt                                                                           <  1:75977/251‑1 (MQ=255)                                                                                                                                                            aTGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGc                                                                       >  2:194541/1‑251 (MQ=255)                                                                                                                                                             agagaTGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCACAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGcc                                                                      >  2:529146/2‑251 (MQ=255)                                                                                                                                                                               gATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATaaa                                                     <  1:281683/250‑1 (MQ=255)                                                                                                                                                                                    aGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATAAAGccc                                                 >  1:253759/1‑249 (MQ=255)                                                                                                                                                                                               tCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATAAAGCCCAGTTCCGCAAt                                      <  2:25021/249‑1 (MQ=255)                                                                                                                                                                                               tCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATGTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCAGAGTTGGTTTGATAAAGCCCAGTTCCGCAATg                                     <  1:219476/250‑1 (MQ=255)                                                                                                                                                                                                      tcAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGCTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATAAAACCCAGTTCCGCAATGGTTTTGAt                             >  2:402255/1‑251 (MQ=255)                                                                                                                                                                                                              cGTTTATCTGCTGAGCCGGTGCGTTGTGTTTATCAACGCGATATGCCCTATGCGCTGGTGGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGGCGGCAAAGCCTGTCTTGTGAACAAACCTCTGCCGCGTTTACGCCATATTTGGTTTGATAAAGCCCAGTTCCGCAATGGTTTTGATGCCATTcg                     <  1:580982/251‑1 (MQ=255)                                                                                                                                                                                                               gTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATAAAGCCCAGTTCCGCAATGGTTTTGATGCCATTcgc                    <  1:194437/251‑1 (MQ=255)                                                                                                                                                                                                                   aaCTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATAAAGCCCAGTTCCGCAATGGTTTTGATGCCATTCGCGAAc                >  1:434940/1‑251 (MQ=255)                                                                                                                                                                                                                          gAGCCTGTGCGTTTTGTTTACCAGCGCGATATGCCCAATGCGCTGGGTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACCCCATCGTCCCCGGCGGGCGTTATCATAAGTTTAAAGACTGTATGAATTTCCCCAATGACGTCAAAGCCGATTTGGTTAATAAACTACGGCGGCGGGGACGGCGGATGTGGGGTGATAAATCCGGCTTCCGCAATGGTTTTGATGCCATCCGCCAACGCGAtg          <  2:252813/250‑1 (MQ=255)                                                                                                                                                                                                                                  gCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATAAAGCCCAGTTCCGCAATGGTGTTGATGCCATTCGCGAACGCGATGTGTTGctc  <  1:204508/250‑1 (MQ=255)                                                                                                                                                                                                                                       ttGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAg                                                                                                          <  2:460146/141‑1 (MQ=255)                                                                                                                                                                                                                                       ttGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAg                                                                                                          >  1:460146/1‑141 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                           ACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATAAAGCCCAGTTCCGCAATGGTTTTGATGCCATTCGCGAACGCGATGTGTTGCTC  >  NC_000913/2623598‑2624072

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