Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	10,659	C→T	6.2%	S233N (AGC→AAC)	yaaW ←	UPF0174 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	10,659	0	C	T	6.2%	158.5 / 2.7	64	S233N (AGC→AAC)	yaaW	UPF0174 family protein
Reads supporting (aligned to +/- strand):  ref base C (25/35);  new base T (2/2);  total (27/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.95e-01

ATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACG  >  NC_000913/10513‑10805                                                                                                                                                   |                                                                                                                                                   atcGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGAGCTa                                                                                                                                                  >  2:921028/1‑149 (MQ=255)  tcGAATGATTGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCTGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTaa                                                                                                                                                 <  2:300476/149‑1 (MQ=255)  tcGAATGAATGTTGAAAAAATTTCCCTGAATAACTTTCGTGAGTTCTGGGCGCGGCATAATAATCAGACAGTGTGGCAGTGTCTACGATCTTTTGAGCGGAAAAAGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCTCTaa                                                                                                                                                 <  2:859393/149‑1 (MQ=255)  tcGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTTTCTACGATCTTTTGAGGGGAAAATGAAATTTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTaa                                                                                                                                                 <  1:9057/149‑1 (MQ=255)  tcGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTaa                                                                                                                                                 <  1:541277/149‑1 (MQ=255)    gaatgaatGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGTGGCAGTGTATACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCAGGTATCAGACATGAGTGGCGCTAAcc                                                                                                                                               >  2:964449/1‑149 (MQ=255)     aatgaatgGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCAT‑ATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTACCCAt                                                                                                                                             >  1:798738/1‑149 (MQ=255)         aatgGTGAAATAATTTCCCTGAATAACTGTAGTGGTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCg                                                                                                                                          <  1:979682/149‑1 (MQ=255)          atgGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCgg                                                                                                                                         <  2:457563/149‑1 (MQ=255)              tGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATATTCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCgca                                                                                                                                     >  2:471978/1‑149 (MQ=255)              tGAAATAATTTCCCTGAATAACTGTAGTGGTTACAGGGCGCGGCATTATAATCAGCCAGTGGGGCAGTGTCTACGATCTGTTGAGGGGAAAATGAATATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCgca                                                                                                                                     >  2:799033/1‑149 (MQ=255)                aaaTAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCGCGGTTTCCGGTATCAGACCGGAGTGGCGCTAACCATCCGGCgcggg                                                                                                                                   >  2:882064/1‑149 (MQ=255)                aaaTAATTTCCCTGAATAACTGTAGTGGTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGGAGTGTCTAGGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCgcagg                                                                                                                                   <  1:115647/149‑1 (MQ=255)                aaaTAATATCCCTGAATAACTGTAGTGTTTTGAGGGCGCGGCATATTAATCAGCCAGTGGGGCAGGGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCgcagg                                                                                                                                   <  1:527863/149‑1 (MQ=255)                 aataatTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCgcaggc                                                                                                                                  <  2:482577/149‑1 (MQ=255)                        tCCCTGAATAACTGTAGTGGTTTCAGGGCGCGGCCTAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGAt                                                                                                                           <  1:221100/149‑1 (MQ=255)                           cTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGAGGCGCTAACCATCCGGCGCAGGCAGGCGATTTg                                                                                                                        >  1:881206/1‑149 (MQ=255)                               aTAACTGTAGTGTTTTCAGGGCGCGGCAAAATAATCAGCCAGTGGGGCAGTGGCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGGGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGt                                                                                                                    <  1:894630/149‑1 (MQ=255)                                      ttgtgtTTTCAGGGCGATGCATAAATATCAGCCTGTGGGGCAGTGTCTACGATCGTTTGAGGGGAAAATGTAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCTTCCGGCGCAGGCAGGCGATTTGCAGTACGGCTg                                                                                                             <  1:567954/147‑1 (MQ=255)                                        ttggtttCAGGGCGCGGCATATTAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAATGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGa                                                                                                           <  1:809543/145‑1 (MQ=255)                                        gtgtTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTGTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGa                                                                                                           >  2:645216/1‑149 (MQ=255)                                        gtgcTTTCAGGGCGCGTCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCAGGTTTCCGTTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGa                                                                                                           <  1:290664/149‑1 (MQ=255)                                         ggtTTTCAGGGCGCGGCATAATCATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGaa                                                                                                          <  2:896245/148‑1 (MQ=255)                                             ttCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGt                                                                                                      <  1:287003/149‑1 (MQ=255)                                                aGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGCCCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCAc                                                                                                   >  2:875776/1‑149 (MQ=255)                                                aGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGAAGTACGGCTGGAATCGTCAc                                                                                                   >  2:832793/1‑149 (MQ=255)                                                    ggcgGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGGGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGGACGGCTGGAATCGTCACGCGa                                                                                               <  1:607008/148‑1 (MQ=255)                                                       ggCATAATAATCAGCCCGTGGGGCAGTGGCTACGAGCTTTTGAGGGGAACAGGAAACTGTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAg                                                                                            <  1:661839/149‑1 (MQ=255)                                                       ggCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGTTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAg                                                                                            >  1:506571/1‑149 (MQ=255)                                                       ggCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGTTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAg                                                                                            >  1:358143/1‑149 (MQ=255)                                                        gCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAgg                                                                                           <  1:842973/149‑1 (MQ=255)                                                        gCATAATAATCAGCCAGTGGGGCAGTGTCTACCATCTTTTGAGGCGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCAGAGAGGCGCTAACCATCCGGAGCAGGCAGACGATTTGCAGTACGGCAAGAATCGTCACGCCATAgg                                                                                           >  1:427798/1‑149 (MQ=255)                                                         cATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGgc                                                                                          >  1:383490/1‑149 (MQ=255)                                                          ataataATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGgcg                                                                                         >  2:409664/1‑149 (MQ=255)                                                              aaatCAGCCCGTGGGCCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTAGCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGc                                                                                     <  2:565146/148‑1 (MQ=255)                                                                    gCCAGTGGGGCAGGGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGa                                                                               <  2:279809/149‑1 (MQ=255)                                                                        ggggggCATGGGCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTCTCCGCCCCGCGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGc                                                                           <  2:103930/147‑1 (MQ=255)                                                                        gTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAAAGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGc                                                                           <  2:107602/149‑1 (MQ=255)                                                                             gCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTACCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCTATTTGAAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTac                                                                      >  1:640926/1‑148 (MQ=255)                                                                             gCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTACCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGCCAGGCGATTTGCAGTACGGCTGGAATCTTCCCGCGATATGCGCTGCCGCTGCCCGATTTaa                                                                      >  1:637181/1‑149 (MQ=255)                                                                                 tgtCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCGGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGACGCTGACCGCTTTAAcgcc                                                                  <  2:985591/149‑1 (MQ=255)                                                                                 tgtCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAAcccc                                                                  <  2:10120/149‑1 (MQ=255)                                                                                    cgaCGATCTTTTCAGGGGAAAATGAACATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCAtt                                                               <  1:559917/147‑1 (MQ=255)                                                                                     tACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTACCCCCAttt                                                              >  2:366994/1‑149 (MQ=255)                                                                                     tACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCAttt                                                              >  2:1006497/1‑149 (MQ=255)                                                                                     tACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGCCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTCGAATCGTCACGCGATAGGCGCTGCCGCTGACAGCTTTAACCCCAttt                                                              >  2:308383/1‑149 (MQ=255)                                                                                        gATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGATTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGGACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGt                                                           >  2:166565/1‑149 (MQ=255)                                                                                          tCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGc                                                         >  2:165674/1‑149 (MQ=255)                                                                                          tCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCGAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGc                                                         <  2:200071/149‑1 (MQ=255)                                                                                               tGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTAGCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGAAGGCGATTTGCAGAACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCAc                                                    <  1:966570/149‑1 (MQ=255)                                                                                                gAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCAcc                                                   <  2:798738/149‑1 (MQ=255)                                                                                                 aGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGTTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCt                                                  <  2:358143/149‑1 (MQ=255)                                                                                                 aGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGTTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCt                                                  <  2:506571/149‑1 (MQ=255)                                                                                                    ggAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGAGCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGCCCGCTTTAACACCATTTAGTGCCGCACCTACa                                               >  1:777047/1‑149 (MQ=255)                                                                                                      aaaatgataatTTTGCCCGGATTCCGGTATCAGTCCTGAGTGGCTCTAACATTCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCGATTTAGTGGCGCACCTACAgg                                             <  1:700061/149‑1 (MQ=255)                                                                                                       aaatgaaaatTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATGTGCAGGACGGCTGGAAACGTCCCGCGATAGGCGCTCCCCCTGACCGCTTTAACCCCCTATAGTGCCGCAACTACAggg                                            >  2:188082/1‑149 (MQ=255)                                                                                                                     cccGGTTTCCGGTCTCAGACGAGAGGGGCGCTACCCATCCGGCGCAGGCAGGCGAGTTGCAGTACGGCGGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCccgc                              <  1:464718/149‑1 (MQ=255)                                                                                                                     cccGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCccgc                              >  1:104553/1‑149 (MQ=255)                                                                                                                     cccGGTCTCCGGTATCAGACCTGAGTGGCGCTACCCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCccgc                              <  2:87990/149‑1 (MQ=255)                                                                                                                          ttgccGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCccgcgccgc                         <  1:895659/146‑1 (MQ=255)                                                                                                                               ggTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCgcagc                    >  1:241617/1‑149 (MQ=255)                                                                                                                                       aCCTTAGTGGCGCTAACCATCCGGCGCAGGCACGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTATTGACGCACCTCCAGGACCTCCCAGCCCCGCGCCGCGCAGCAAACCATg            >  2:795525/1‑149 (MQ=255)                                                                                                                                        ccTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGc           <  2:838350/149‑1 (MQ=255)                                                                                                                                              gggCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGt     <  1:471978/148‑1 (MQ=255)                                                                                                                                                 cgcTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACg  <  2:541002/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   ATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACG  >  NC_000913/10513‑10805

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