Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,730,870	T→C	noncoding (288/1542 nt)	rrsG ←	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,730,870	0	T	C	97.7%	95.3 / ‑3.6	43	noncoding (288/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (22/20);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

CCTTCCTCCCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGT  >  NC_000913/2730704‑2731105                                                                                                                                                                       |                                                                                                                                                                                                                                            ccTTCCTCCCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCCGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCCAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAggac                                                                                                                                                         >  2:298988/1‑249 (MQ=14)    ccctccCCGCTGAAAGTACTTTACACCCCGAAGGCCTTCGTCATACCCGCGGCATGGCTGCATCAGGCTTGCGCCCATTGGGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTTCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCACGAAGGTcccc                                                                                                                                                      <  2:295767/250‑1 (MQ=11)     cctccCCTCTGAAAGTACTTTACCACCCGATGGCCTTCTTCATTCTCGCTTCATGGCTTCATCAGGCTTGCGCCCATTGTGCAATATTCCCCGCTGCTCCCGCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCGCTCAGACCAGCCAGGGATCGTCGCCTTGTTGAGCCGTTACCCCACCAACAAGCTACTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTacccc                                                                                                                                                     <  2:108128/251‑1 (MQ=11)     cctccCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTccccc                                                                                                                                                     >  1:26286/1‑251 (MQ=14)                aaaGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGCCCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGCTGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCtt                                                                                                                                          >  2:172874/1‑251 (MQ=14)                 aaGTACTTTACAAACCGAAAGCCTTCTTCATACACGCGGCATGGCTGCTTCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTACATCACGTAGGAGTCTGGGCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTCATCCCATCTGGGCACATACGATGGCAAGAGGCCCGAAGGTCTCCCCCTTTgggctcg                                                                                                                                         >  2:32610/1‑246 (MQ=11)                      cTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTTGGCACATCCGATGGCAAGAGGCCCGCAAGTCCCCCTCTTTGGTCTTGCgcc                                                                                                                                     >  2:433806/1‑248 (MQ=9)                                cGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGt                                                                                                                           >  1:228426/1‑250 (MQ=18)                                     gCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGc                                                                                                                     >  1:152580/1‑251 (MQ=18)                                       cttcttCATACACGCAGCATGACTGCATCAGGCTTGAGCCCATTGTGCAATATTCCTCACTGCTGCCTACCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGttttcgct                                                                                                                    >  2:407319/1‑243 (MQ=11)                                         tcttcATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCCCCAACAAGCTAATCCCATCTGGGCACATCCGCTGGCAAGAGGCCCGCAGCTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTAcc                                                                                                                 >  1:131413/1‑251 (MQ=255)                                             cATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGttt                                                                                                             <  1:98158/251‑1 (MQ=18)                                                acacGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTccc                                                                                                          >  2:472565/1‑250 (MQ=18)                                                 cacGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTcc                                                                                                           >  1:311517/1‑249 (MQ=18)                                                    gctgCATGGCCGCCTCCGGCGTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGGAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCTTCCTCTCAGCCCAGCCAGGGATCGACGCCTTGTTGCGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCCCATCCGATGGCAAGAGGCCCGTAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagtag                                                                                                      <  2:268751/248‑1 (MQ=11)                                                       gCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGtt                                                                                                    >  1:76150/1‑250 (MQ=18)                                                             cTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATccccc                                                                                             >  1:220393/1‑251 (MQ=18)                                                                   cAGGCTTGCGCCCATTGTGCAATATCCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCGCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCATTAGTTATCCCCCTCCATc                                                                                       <  2:415083/251‑1 (MQ=16)                                                                      gCTTGCGCCCATTGTGCAATATTCCCCCCTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAgg                                                                                    <  2:102814/251‑1 (MQ=18)                                                                            gcCCATTGTGCAATATTCCCCACTGCTGCCTNCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCtt                                                                              <  2:247283/251‑1 (MQ=21)                                                                             cccATTGTGCAATATTCCCCCCTGCTGCCTCNCGTAGGAGTCTGGACCGTGTCGCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTc                                                                             <  2:225322/251‑1 (MQ=19)                                                                                            tCCCCACTGCTGCCTCCCGTAGGACTCTGGACCGTGTCTCAGTTCCCGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTCGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCAc                                                              >  2:427789/1‑251 (MQ=255)                                                                                                ccctgatgcCTCCCGTAGGAGTCTGGACCGGGTCTCAGTTCCAGTGTGGCTGGTCAGCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACAACCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGTTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGt                                                          <  2:136746/245‑1 (MQ=11)                                                                                                 cctgctgCCTCCCGTAGGAGTCGGGACCGGGTCTCAGTTCCAGTGTGGCTGGGCATCCTCTCATGCCAGCCAGGGATCGTCGCCGTGGTGAGCCGTTCCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTGTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTGCCAGTAGTTATCCCCCCCCATCAGGCAGCTTCCCAGACATTACTCACCCGTc                                                         <  2:178498/250‑1 (MQ=11)                                                                                                    gctgcCTCCCGTCTGAGTCTTGACCGTTTCTCAGTTCCAGTTTGGCTGGTCTTCCTCTCAGACCAGCCAGGGATCGTCTCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCTCATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCg                                                       <  2:1547/250‑1 (MQ=11)                                                                                                       gcCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCAc                                                   >  1:46447/1‑251 (MQ=21)                                                                                                         cccccGTAGGAGGCTGGACCGTGTCTCATTTCCAGTGTGGCGGGTCACCCTCCCAGACCAGCCAGGGTTCGTCGCCTTGGTGAGCCGTTCCCCCCCCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTCCCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTc                                                 <  2:11810/249‑1 (MQ=11)                                                                                                               gaGGAGTCTGGACCGTGTCACAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGc                                           <  1:181315/250‑1 (MQ=21)                                                                                                                aGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCCACCCGCTAATCCCATCTGGGCACATCCGTTGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGt                                                                                                                                             >  2:138492/1‑152 (MQ=11)                                                                                                                aGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGt                                                                                                                                             <  1:138492/152‑1 (MQ=17)                                                                                                                         gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGa                                                                                                                                      >  2:165780/1‑150 (MQ=17)                                                                                                                         gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGa                                                                                                                                      <  1:165780/150‑1 (MQ=17)                                                                                                                             gtgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAgctg                             <  1:456076/251‑1 (MQ=37)                                                                                                                               gtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAgctgct                           >  2:142552/1‑251 (MQ=38)                                                                                                                               gtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACACCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAgctgct                           >  1:289886/1‑251 (MQ=255)                                                                                                                               gtCTCAGTCCCTGTGTGGCGGGTCACCCCCCAAAACCAGCCAGGGATGGTCGCCTTGGTGAGCCGTTACCCCACCAACAACCTACTCCCATCTGGGCCCACCCGATGGCAGGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCTACAGTATGCGGTATTAGCTCCCGTTTCCAGTAGTTATCCCCCCCCACCAGGCAGCGCCCCAGGCATGACTCACCCTTCCGCCACTCGTCAGAGAAGCAGCAAgctgct                           <  2:161671/251‑1 (MQ=2)                                                                                                                                tctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGAc                                                                                                                                     <  2:229054/144‑1 (MQ=17)                                                                                                                                tctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGAc                                                                                                                                     >  1:229054/1‑144 (MQ=17)                                                                                                                                                   ggTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACt         <  1:404670/249‑1 (MQ=39)                                                                                                                                                   ggTCATCCTCTCAGACCACCCATGGATCGTCGCCTTGGTAAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCACGAGGCCCGTAGGCCACCCTCTTTGGTCATGCGACATTATGCGGTATTAGCTACCGTTTCCCGTCGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGNACCAAGCTGCTTCCTGTTACCGCTCGACtcg       >  2:244626/1‑249 (MQ=14)                                                                                                                                                    gTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGTCACGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCCGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGCAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGc      >  1:77303/1‑251 (MQ=255)                                                                                                                                                      cATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCAt    >  2:257320/1‑251 (MQ=38)                                                                                                                                                        ccccctcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAACCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCAtgt  <  1:140162/247‑1 (MQ=255)                                                                                                                                                                       |                                                                                                                                                                                                                                            CCTTCCTCCCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGT  >  NC_000913/2730704‑2731105

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