Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	872,248	C→T	6.2%	T120M (ACG→ATG)	gsiD →	glutathione ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	872,248	0	C	T	6.2%	202.3 / 2.6	65	T120M (ACG→ATG)	gsiD	glutathione ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (27/33);  new base T (2/2);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

TCTGAATAACGGACCTTCTTTGCAGCACTGGTTTGGCGTCGATTCACTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGT  >  NC_000913/872108‑872389                                                                                                                                             |                                                                                                                                              tCTGAATAACGGACCTTCTTTGCAGCACTGGTTTGGCGTCGATTCACTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGtt                                                                                                                                            <  8:263978/144‑1 (MQ=255) tCTGAATAACGGACCTTCTTTGCAGCACTGGTTTGGCGTCGATTCACTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGtt                                                                                                                                            >  7:263978/1‑144 (MQ=255) tCTGAATAACGGACCTTCTTTGCAGCACTGGTTTGGCGTCGATTCACTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGGGTTTATCGGTGCGGCGATCGGGACGTTGCTg                                                                                                                                        >  8:186094/1‑148 (MQ=255)    gAATAACGGACCTTCTTTGCAGCACTGGTTTGGCGTCGATTCACTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGc                                                                                                                                     >  5:171046/1‑148 (MQ=255)     aaTAACGGACCTTCTTTGCAGCACTGGTTTGGCGTCGATTCACTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGc                                                                                                                                     <  1:194476/147‑1 (MQ=255)     aaTAACGGACCTTCTTTGCAGCACTGGTTTGGCGTCGATTCACTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCt                                                                                                                                    <  1:252987/148‑1 (MQ=255)     aaTAACGGACCTTCTTTGCAGCACTGGTTTGGCGTCGATTCACTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCt                                                                                                                                    <  4:263770/148‑1 (MQ=255)      aTAACGGACCTTCTTTGCAGCACTGGTTTGGCGTCGATTCACTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCtt                                                                                                                                   <  6:190793/148‑1 (MQ=255)          cGGACCTTCGTTGCAGCACTGGTTTGGAGTCGATTCACTGGGGATTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGAGCGTGTTTGCCGTTTTTATCCGTCCGGAGATCCGGAAGTTGCTTGGCTTGCTc                                                                                                                               >  7:327158/1‑148 (MQ=255)             aCCTTCTTTGCAGCACTGGTTTGGCGTCGATTCACTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCt                                                                                                                            >  1:404613/1‑148 (MQ=255)                     tGCAGCACTGGTTTGGCGTCGATTCACTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGAtatt                                                                                                                     <  2:404613/147‑1 (MQ=255)                        agcaCTGGTTTGGCGTCGATTCACTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGa                                                                                                                 <  8:447/148‑1 (MQ=255)                        agcaCTGGTTTGGCGTCGATTCACTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCGCGCTGGATATTATGa                                                                                                                 <  2:76455/148‑1 (MQ=255)                               gTTTGGCGTCGATTCACTGGGGCGGGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCtg                                                                                                           <  3:218165/147‑1 (MQ=255)                                  tGGCGTCGATTCACTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCtggtgg                                                                                                       >  3:302644/1‑148 (MQ=255)                                       tCGATTCACTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCg                                                                                                  <  7:186094/148‑1 (MQ=255)                                           ttgaCTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGGGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCGCGCTGGATATTATGAAGGCTGGTGGGATCGGCTg                                                                                              <  2:66835/145‑1 (MQ=255)                                           ttCACTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTg                                                                                              <  1:394572/148‑1 (MQ=255)                                              aCTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATc                                                                                           <  6:171046/148‑1 (MQ=255)                                                tGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCAt                                                                                         >  2:423518/1‑148 (MQ=255)                                                      gTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCAt                                                                                   >  5:7120/1‑148 (MQ=255)                                                        gACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGtt                                                                                                                                            >  5:299078/1‑89 (MQ=255)                                                        gACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGtt                                                                                                                                            <  6:299078/89‑1 (MQ=255)                                                        gACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGc                                                                                                             >  1:41833/1‑120 (MQ=255)                                                        gACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGc                                                                                                             <  2:41833/120‑1 (MQ=255)                                                        gACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTg                                                                                              >  4:187370/1‑135 (MQ=255)                                                        gACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTg                                                                                              <  3:187370/135‑1 (MQ=255)                                                        gACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCAttt                                                                                 >  3:316774/1‑148 (MQ=255)                                                          cATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGAt                                                                                            <  2:101753/135‑1 (MQ=255)                                                          cATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGAt                                                                                            >  1:101753/1‑135 (MQ=255)                                                           tttttCTGCCGTGTCCTGGTTGGTGCGCATATGTCTCTGGCGGCGGGCGGGTTTTCCGTGTTTATGGGTGCGGCGATCGTGACGTTGCGGGGCTTGCTCGCTGGATATTATCAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCg                                                                              <  1:129651/147‑1 (MQ=255)                                                             tttCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGATGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATc                                                                                           <  1:300841/133‑1 (MQ=255)                                                             tttCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGATGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATc                                                                                           >  2:300841/1‑133 (MQ=255)                                                                            ggttggtGCGCAAATCTCGCTGGCGGCGGGCGTTTTTGACGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTc                                                             >  2:371580/1‑148 (MQ=255)                                                                              ttggtGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCgg                                                                                                 <  7:211461/110‑1 (MQ=255)                                                                              ttggtGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCgg                                                                                                 >  8:211461/1‑110 (MQ=255)                                                                                ggtGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCg                                                          <  4:302644/147‑1 (MQ=255)                                                                                  tGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTg                                                                                              <  2:49925/109‑1 (MQ=255)                                                                                  tGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTg                                                                                              >  1:49925/1‑109 (MQ=255)                                                                                     gcAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGATGTTGCTGGGCTTGCTCGCTGGATTTTATGa                                                                                                                 <  2:322215/87‑1 (MQ=255)                                                                                     gcAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGATGTTGCTGGGCTTGCTCGCTGGATATTATGa                                                                                                                 >  1:322215/1‑87 (MQ=255)                                                                                     gcAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCg                                                          <  7:214333/142‑1 (MQ=255)                                                                                     gcAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCg                                                          >  8:214333/1‑142 (MQ=255)                                                                                      cAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCtg                                                                                                           >  2:349270/1‑92 (MQ=255)                                                                                      cAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCtg                                                                                                           <  1:349270/92‑1 (MQ=255)                                                                                       aaaTCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTAtttt                                                  <  1:423518/148‑1 (MQ=255)                                                                                         aTCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGAtatt                                                                                                                     <  3:55930/79‑1 (MQ=255)                                                                                         aTCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGAtatt                                                                                                                     >  4:55930/1‑79 (MQ=255)                                                                                                    ggcgggcgTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCt                                     <  6:7120/148‑1 (MQ=255)                                                                                                     gcgggcgTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCtg                                    <  4:316774/148‑1 (MQ=255)                                                                                                             gtTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCg                            >  8:4053/1‑148 (MQ=255)                                                                                                              tttGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTTGGGATCGCTGTTGTTGCgg                           >  6:242352/1‑148 (MQ=255)                                                                                                              tttGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTAGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCgg                           >  5:12159/1‑148 (MQ=255)                                                                                                                     tgtgTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCGGTTGTTGCGGGGTTggg                    <  1:371580/148‑1 (MQ=255)                                                                                                                       tttATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCg                                                                              >  1:143451/1‑88 (MQ=255)                                                                                                                       tttATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCg                                                                              <  2:143451/88‑1 (MQ=255)                                                                                                                        ttATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCAttt                                                                                 <  2:213153/84‑1 (MQ=255)                                                                                                                        ttATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCAttt                                                                                 >  1:213153/1‑84 (MQ=255)                                                                                                                        ttATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGgtg                         <  7:341703/140‑1 (MQ=255)                                                                                                                        ttATCGGTGCGGCGATCGGGACGTTCCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGTCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGgtg                         >  8:341703/1‑140 (MQ=255)                                                                                                                             ggTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGCGTTGAGAAGCGGCa            <  6:12159/148‑1 (MQ=255)                                                                                                                               tGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATAATGCGCATTTGCGATGTGCTGTTTGCCTTCCAGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGGTGGGAAGCGGCAtt          >  2:104314/1‑148 (MQ=255)                                                                                                                                       tCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTTGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGt  >  5:390747/1‑148 (MQ=255)                                                                                                                                           gACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTAtt                                                    >  6:289691/1‑94 (MQ=255)                                                                                                                                           gACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTAtt                                                    <  5:289691/94‑1 (MQ=255)                                                                                                                                             |                                                                                                                                              TCTGAATAACGGACCTTCTTTGCAGCACTGGTTTGGCGTCGATTCACTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGT  >  NC_000913/872108‑872389

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