Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,907,375	+T	coding (61/210 nt)	cspC ←	transcription antiterminator and regulator of mRNA stability CspC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,907,375	1	.	T	97.4%	106.0 / ‑4.1	38	coding (61/210 nt)	cspC	transcription antiterminator and regulator of mRNA stability CspC
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base T (16/21);  total (17/21)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGT‑AA‑TGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGT‑CAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCG  >  NC_000913/1907142‑1907618                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      aTAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAACCt                                                                                                                                                                                                                                      >  1:356814/1‑251 (MQ=255)      gCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAACCTTTAGa                                                                                                                                                                                                                                 <  1:419656/251‑1 (MQ=255)          aGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAACCTTTAGACTCg                                                                                                                                                                                                                             <  2:201269/251‑1 (MQ=255)                        cATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGAACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTAACCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAAc                                                                                                                                                                                                               >  1:500289/1‑251 (MQ=255)                            cgacgaGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTg                                                                                                                                                                                                            >  2:172341/1‑250 (MQ=255)                            cgacgaGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTg                                                                                                                                                                                                            <  1:43844/250‑1 (MQ=255)                              acgaGGCTTAATAGCGAGAGTGAAGCGTATTCTCACTTCAGATCTGTGGATTCGATCAGATAGCTGTTACGTTAACATCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTg                                                                                                                                                                                                            <  2:475551/248‑1 (MQ=255)                                  ggCTTTATAGCGAGATTGAAGGCTATTCACACTTCAGATCAGTGGATTTGAACAGATAGCTGTGACGTAAACAGCTGCCGGACGATTCTGGCCGTCCTGAATTTCGAACGCAACGTTCTGACCTTCAGCCAGAGTTTTGAATCCATTACCCTGGATAGCGGAGAAGTGTACGAACTCATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAACCTTTAGACGCGTTGAACCACTTAACCTGACCTTTa                                                                                                                                                                                                     <  2:221077/251‑1 (MQ=255)                                    cTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAAt                                                                                                                                                                                                   >  1:100195/1‑251 (MQ=255)                                     tttATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATc                                                                                                                                                                                                  >  2:245519/1‑251 (MQ=255)                                          aGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTg                                                                                                                                                                                              <  2:337935/250‑1 (MQ=255)                                          aGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGc                                                                                                                                                                                             >  1:143848/1‑251 (MQ=255)                                                 ttGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAACCTTTAGACTCGTTGa                                                                                                                                                                                                                         >  1:209218/1‑216 (MQ=255)                                                 ttGAAGCGTATTCACACTTCAGAGCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAACCTTTAGACTCGTTGa                                                                                                                                                                                                                         <  2:209218/216‑1 (MQ=255)                                                           ttCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACATTTCTGGCCGTCCTGAATTTCGAACTAAACGTTCTGACCTTCAGCCAGCGTTTTGAAGCCATTACCATGGATAGCGGAGCAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCtt                                                                                                                                                                            >  2:237672/1‑251 (MQ=255)                                                             cacaCTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTGCTGGCCGTCCTGAATTTCGAACTCAACGGTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTa                                                                                                                                                                           <  1:278313/250‑1 (MQ=255)                                                                         cAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTttctt                                                                                                                                                               <  1:388940/250‑1 (MQ=255)                                                                                         ggaGCTGTTACGTTAACAGCTGCCGGACCTTTCTAGACGACCTGTATGTCGAACTCAACGATCTGACCTTCAGCCAGAGATTTGAAGCCATTACCCTGGATAGCGACGAAGAGTACGAACACATCTTTGCTGCCATCAGCCGGAGTTAA‑TGGAACACAAACCTTAAGACACGATGAACCACTTAACCTGCACTTTAATATTTGCCATTTGAAAAATTACTTAGATGGTTTTCTTCGCCCgacggaataac                                                                                                                                               >  2:378602/3‑240 (MQ=255)                                                                                         aTAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGACGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTTCCCTGGATAGCGGAGAAGTGTACGAACTCATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACa                                                                                                                                              <  2:157839/251‑1 (MQ=255)                                                                                           aGCTGTTACTTAAACGGCTGCCGTACCTTTCTGGCCGTCCTGAATTTCGAACCCAACGGTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTGCTTCGCCCGCAGGCATAACATa                                                                                                                                            <  2:51731/251‑1 (MQ=255)                                                                                                   cGTTAACAGCTGCCAGACCTTTATGGCCGTCCTGAATTTCGAACTAAACGTTCTGACCTTCAACCAGAGCTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGCACACATCTTTGCTGCCATCAGCAGGAGTTAA‑TGAAGCCAAAACCTTCAGACTCGTTGAACCACTTAACCTGACCTTTCATCTTTGCAATTTGATAAATTCCTTCGCTTGTTTTCTTCACCCGCAGTCATacaaacgataaaac                                                                                                                                    >  2:354629/1‑238 (MQ=255)                                                                                                                         cTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTTGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGcc                                                                                                              >  1:514079/1‑250 (MQ=255)                                                                                                                           ggCCGTCCTGAATTTCGAACCCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGGACGAACACATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCTTGAGGCAcc                                                                                                            <  2:308519/251‑1 (MQ=255)                                                                                                                           ggCCGTCCTGAATGTCGAACTCAACGTTCGGACCGTCAGCCAGATTTTTGAAGCCATTACCCTGGATAGCGGATAAGTGTGCGAACACATCTATGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAGCCTTTAGACTCGTTGAACCCCTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCCCCCGCACGCATAACATAGATAAAACTGACACAGTACTGCATGAGGc                                                                                                               <  1:490653/248‑1 (MQ=255)                                                                                                                               gTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGTCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAACTGTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCGTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCatta                                                                                                        <  2:356814/251‑4 (MQ=255)                                                                                                                                     aaTTGCGAACTCAACGTTCTGAACTTCAGCCAGAGTTTTTAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTGTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGc                                                                                                  <  2:81642/251‑1 (MQ=255)                                                                                                                                         tCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAACCGTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCATGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCg                                                                                               <  2:100195/250‑1 (MQ=255)                                                                                                                                          cGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCgg                                                                                              <  1:245519/250‑1 (MQ=255)                                                                                                                                                cAACGTTCTGACCTGCAGCCAGAGGTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTCCGAACACATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAACCTTTAGACTCGTTGAACCTCTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATATCATAGATAAAACTGACACATTACTTCAGGAGGCACCAATATAATGCTCGGCagag                                                                                         <  1:306517/249‑1 (MQ=255)                                                                                                                                                       cggACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCTGTa                                                                                 <  2:143848/248‑1 (MQ=255)                                                                                                                                                                                  aTTACCCTGGATAGCGAAGCAGAGTACGAACACATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAAACTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAACATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACAGACACATTACTGCATGAGGCACAAATATAAGGCTAGGCAGAGAAGCGGTATTCAACGT‑AAACGTGTTTACTCAGGAct                                                    >  2:529738/1‑251 (MQ=255)                                                                                                                                                                                      cccTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGT‑AAATGAAGCCAGAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGAGGTATTCATAGTACAA‑GTGTTTACTAAGGAattctt                                                >  2:445094/1‑251 (MQ=255)                                                                                                                                                                                                cGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTcaacgt‑caacgta                                                                 >  1:21841/1‑223 (MQ=255)                                                                                                                                                                                                cGGAGAAGTGTACGAACACATCTGTGCTGCCATCTGCCGGAGTTAA‑TGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGT‑CAACgtg                                                                 <  2:21841/224‑1 (MQ=255)                                                                                                                                                                                                  gagaAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGT‑CAACGTGTTTACTCAGGCCTTCTTTACTGAAAATGc                                    >  1:264401/1‑251 (MQ=255)                                                                                                                                                                                                             gAACACATCTTTGCTGCCATCAGCCGGAGTTAA‑TGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGT‑CAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATaaa                          >  1:474875/1‑250 (MQ=255)                                                                                                                                                                                                                   aTCTTTGATGACATCAGCAGGAGTTAA‑TGAAGCCAAACCCTTTAGACTCGTTAAACCACTTAACCTGAACTTTAATCATTGCCATTTGAAACATTCCTTCGATTGTTTTCTTCGCCCGAAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAAGATAAGGCTCTGCAGAGAAGCGATATTCAACGT‑CAACGTGTTTACTCAGGACGTCTTTACTATAAATGCcaaccataaccagtaat                   >  2:217099/1‑236 (MQ=255)                                                                                                                                                                                                                                    ccGGAGTTAA‑TGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGT‑CAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCg  <  1:109723/251‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      ATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGT‑AA‑TGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGT‑CAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCG  >  NC_000913/1907142‑1907618

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