Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	7,573	C→T	6.3%	P129P (CCG→CCA)	yaaJ ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	7,573	0	C	T	6.3%	300.9 / 6.5	95	P129P (CCG→CCA)	yaaJ	putative transporter
Reads supporting (aligned to +/- strand):  ref base C (45/44);  new base T (3/3);  total (48/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

CGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGC  >  NC_000913/7426‑7717                                                                                                                                                    |                                                                                                                                                 cGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCg                                                                                                                                                 <  7:4478/149‑1 (MQ=255) cGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCg                                                                                                                                                 <  3:77119/149‑1 (MQ=255)   ggggAAAATCAAAAGAAAAACTCATGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATTCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCCGCGCCTCCCCAGCCCGCGCGCCATATCCCATGCTGGt                                                                                                                                               >  8:153202/1‑149 (MQ=255)    gggAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTc                                                                                                                                              <  4:375330/149‑1 (MQ=255)    gggAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTc                                                                                                                                              <  6:185696/149‑1 (MQ=255)     ggAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTcc                                                                                                                                             >  6:374724/1‑149 (MQ=255)     ggAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTcc                                                                                                                                             >  5:351038/1‑149 (MQ=255)           taaaaaGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACACAACGCCCATCCAGCGCATCCCCACCCCCCGCGCCATAAACCATGCCGGTCCGCCAc                                                                                                                                        >  6:407910/3‑148 (MQ=255)           tCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCTGGTCCGCCAc                                                                                                                                        >  8:206410/1‑148 (MQ=255)           tCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTcc                                                                                                                                             >  5:101013/1‑143 (MQ=255)           tCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTcc                                                                                                                                             <  6:101013/143‑1 (MQ=255)             aaaaGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGa                                                                                                                                      >  6:128726/1‑148 (MQ=255)                 gAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTg                                                                                                                                 >  1:381633/1‑149 (MQ=255)                  aaaaaCTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGACCATCCCGCGCATCCCCAGCCCGAGCGCCATATACCATGCCGGTCCGCCACGAAACTGc                                                                                                                                >  5:239815/1‑149 (MQ=255)                    aaaCTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCCCGAAACTGcc                                                                                                                               >  5:3123/1‑148 (MQ=255)                       cTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGaaa                                                                                                                                    <  2:274368/140‑1 (MQ=255)                       cTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAAcagatcgg                                                                                                                            >  1:274368/1‑141 (MQ=255)                       cTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCAtt                                                                                                                           <  6:37108/149‑1 (MQ=255)                        tCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCAtt                                                                                                                           >  3:307267/1‑148 (MQ=255)                          aGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGAc                                                                                                                        <  1:139227/149‑1 (MQ=255)                             gcgcGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGt                                                                                                                      >  4:332311/1‑148 (MQ=255)                                  ggCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCCGCCCGCGCGCCATATACCATGCCGGTCCGCCACGACACTGCCCATTGCCGGCACGtt                                                                                                                >  8:417321/1‑149 (MQ=255)                                        ggCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCTGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTtat                                                                                                          <  7:206410/149‑1 (MQ=255)                                                cTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGTCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTtgt                                                                                                  <  2:381633/149‑1 (MQ=255)                                                cTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTtgt                                                                                                  <  6:335127/149‑1 (MQ=255)                                                 ttGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGtt                                                                                                    <  7:122901/146‑1 (MQ=255)                                                 ttGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGtt                                                                                                    >  8:122901/1‑146 (MQ=255)                                                 ttGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTtgtg                                                                                                 >  6:273398/1‑149 (MQ=255)                                                 ttGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTtgtg                                                                                                 >  2:360116/1‑149 (MQ=255)                                                 ttGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTGCTTTATAAAGTtgtg                                                                                                 >  8:191970/1‑149 (MQ=255)                                                    aaCTCCGCTGAAAATTCTGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCTGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGcc                                                                                               <  7:153202/148‑1 (MQ=255)                                                    aaCTACGCTGAAAATTATGCCATAGGCTATGAGCAAAAAGACGGCGCACAGAACGCCCATCCAGCGCATCCCCATCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGc                                                                                                >  6:371575/1‑147 (MQ=255)                                                      cTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCaga                                                                                            <  7:367171/149‑1 (MQ=255)                                                         cGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCAAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCagag                                                                                           <  2:61575/147‑1 (MQ=255)                                                            tGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCCGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACAt                                                                                      >  3:184650/1‑149 (MQ=255)                                                             gAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACAtt                                                                                     >  8:436721/1‑149 (MQ=255)                                                             gAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACAtt                                                                                     >  2:125752/1‑149 (MQ=255)                                                              aaaaTTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTc                                                                                    <  7:287885/149‑1 (MQ=255)                                                              aaaaTTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCAAGCGCATCACCAGCACGCGCGCCATATACCATGCCGGTCCGCAACGAAACTGCACATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTc                                                                                    >  5:17187/1‑149 (MQ=255)                                                               aaaTTTTGCCATAGGCGATGAGCAAAAAGACGGCGACCAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGAGGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCg                                                                                   <  6:127706/149‑1 (MQ=255)                                                               aaaTTATGCCATAGGCGATGAGCAAAACGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCg                                                                                   <  6:351038/149‑1 (MQ=255)                                                                        cATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGCACCTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCagagcg                                                                           >  7:243515/1‑143 (MQ=255)                                                                        cATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCgg                                                                                  <  8:243515/141‑1 (MQ=255)                                                                        cATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAg                                                                          >  7:341991/1‑149 (MQ=255)                                                                        cATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAg                                                                          <  5:273398/149‑1 (MQ=255)                                                                         aTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAgg                                                                         >  7:268262/1‑149 (MQ=255)                                                                             gCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCGTTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGc                                                                     >  5:171202/1‑149 (MQ=255)                                                                               gATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCa                                                                   <  5:374724/149‑1 (MQ=255)                                                                                   aGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGc                                                                <  1:125752/148‑1 (MQ=255)                                                                                   aGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTGCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGc                                                                <  6:244924/148‑1 (MQ=255)                                                                                      aaaaaCACGGCGAACACAACGCCCAGCCAGCGCATCCCCAGCCCGCGCGCCATAAACAATGCCGGGCCGCCACCAAACTGCCCATTGACGCCACGTGCTTTATAAAGTTGTGCCCCAGACCATTCGGCAAACGAGGTCGCCATGCCGAt                                                            <  6:3123/149‑1 (MQ=255)                                                                                        aaaGACGGCGAACAGAACGCCCCTCCAGCGCATCCCCAGCCCACGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGCCGTCACGTTCTGTATAAAGTTGTGCCAGAGAACATTCGGAAAACGAGGTCGCCATGCCGATa                                                           >  1:209167/1‑148 (MQ=255)                                                                                          aGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGTGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAAc                                                        >  4:279629/1‑149 (MQ=255)                                                                                             cGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAAcgcg                                                     <  3:362701/149‑1 (MQ=255)                                                                                             cGGCGAACAGAACGCCCATCCAGCGCATCCCAAGCCCCCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCCCGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAAcgc                                                      <  6:239815/148‑1 (MQ=255)                                                                                              ggCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCgg                                                    >  6:100140/1‑149 (MQ=255)                                                                                                 gAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCaa                                                 >  6:299296/1‑149 (MQ=255)                                                                                                       aaCGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGACGGTCCGCCACGAAACTGCCCATTGACGTCAAGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCca                                             >  1:249271/1‑147 (MQ=255)                                                                                                        aCGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACAt                                          <  1:360116/149‑1 (MQ=255)                                                                                                         cGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATc                                         >  8:70792/1‑149 (MQ=255)                                                                                                           cccATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATcc                                        >  2:58367/1‑148 (MQ=255)                                                                                                            ccATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACAATCGGCAAACGAGGTCGCCATGCCGATAAACTCGGCAACCCACATCCa                                       >  4:302766/1‑148 (MQ=255)                                                                                                              aTCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCaaaa                                    <  2:124240/149‑1 (MQ=255)                                                                                                              aTCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCaaaa                                    <  3:18784/149‑1 (MQ=255)                                                                                                              aTCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCaaaa                                    <  4:307267/149‑1 (MQ=255)                                                                                                               tCCAGCGCATCCCCAGCCCGCGCGCCATATACCCTGCCGGTCCGCCCCGAAACTGCCCATTGACGGCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCACCCCACATCCaa                                      >  3:37762/1‑146 (MQ=255)                                                                                                               tCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCTGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCaaag                                    >  3:306307/1‑147 (MQ=255)                                                                                                               tCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGtt                                                                                                    >  2:94720/1‑84 (MQ=255)                                                                                                               tCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGtt                                                                                                    <  1:94720/84‑1 (MQ=255)                                                                                                                 cAGCTCATCCCCAGCCCGCGCGCCATATACGATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGAc                                 <  5:407910/149‑1 (MQ=255)                                                                                                                 cAGCGCATCCCCCGCCCGCGCGCCATATACCATGCCGGGCCGCCACGAAACTGCCCATTCACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGAc                                 <  1:193362/149‑1 (MQ=255)                                                                                                                  aGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTCCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACg                                >  6:164469/1‑149 (MQ=255)                                                                                                                       aaccacAGCCCGCGCGCCATATACCATGCCGTGCCTCCACGAAACAGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCTGTGAAGATTCGGCAACCGAGGGCGCCCTGCCGATACACGCGGCAACCCACATCCACACGACGGCTcc                           <  2:100973/144‑1 (MQ=255)                                                                                                                       aTCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCGTTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACACGGCAACCCACATCCAAAAGACGGCTcc                           <  6:171202/149‑1 (MQ=255)                                                                                                                        tCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATACACGCGGCAACCCACATCCAAAAGACGGCTc                            >  8:176169/1‑147 (MQ=255)                                                                                                                         ccccCGCCCGCGCGCCCTAAACCCTCCCCCCCCGCCACGAAACTGCCCATTGACCTCCCGTCCTTTATAAAGGAGCGCCAGCGAACATTCGGCAAACGAGGGCGCCATGCCGATAAACCCGGCAACCCACATCCAAAAGACGGCTCCAg                         <  8:268262/149‑1 (MQ=255)                                                                                                                          cccAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAgg                        <  3:332311/149‑1 (MQ=255)                                                                                                                           ccAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCACCCGAGGTCGCCATGCCGATAAACGCGGCACCCCACCTCCAACAGCCGGCTCCAGGt                       >  6:163199/1‑149 (MQ=255)                                                                                                                           ccAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAg                         <  8:436846/147‑1 (MQ=255)                                                                                                                           ccAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAg                         >  7:436846/1‑147 (MQ=255)                                                                                                                                    gcgcCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACcggcg              <  8:341991/149‑1 (MQ=255)                                                                                                                                     cgcCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGATAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACcggcgg             <  5:371575/149‑1 (MQ=255)                                                                                                                                     cgcCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTGCACcggcgg             <  5:128726/149‑1 (MQ=255)                                                                                                                                          tataCCATGCTGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAGGACGGCTCCAGGTCCACCGGCGGTAATa        <  4:306307/149‑1 (MQ=255)                                                                                                                                          tataCCATGCCGGTCCGCCACGTAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATa        <  4:184650/149‑1 (MQ=255)                                                                                                                                           ataCCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATTAACGCGGCACCCCACCTCCAAAAGACGGCTCCAGGTCCACCGGCGGTTATAg       >  5:309439/1‑149 (MQ=255)                                                                                                                                            taCCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAg       <  3:279629/148‑1 (MQ=255)                                                                                                                                             aCCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGc      >  7:241828/1‑148 (MQ=255)                                                                                                                                              ccATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCa    <  7:436721/149‑1 (MQ=255)                                                                                                                                               cATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAg   >  8:329587/1‑149 (MQ=255)                                                                                                                                                aTGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAgc  >  3:191874/1‑149 (MQ=255)                                                                                                                                                 tGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTaa          >  1:206135/1‑140 (MQ=255)                                                                                                                                                 tGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTaa          <  2:206135/140‑1 (MQ=255)                                                                                                                                                  gCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCa    >  2:4858/1‑145 (MQ=255)                                                                                                                                                  gCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCa    <  1:4858/145‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 CGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGC  >  NC_000913/7426‑7717

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