Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,894,391	A→C	8.2%	L223R (CTG→CGG) ‡	yegH ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,894,391	0	A	C	8.2%	220.5 / 8.4	71	L223R (CTG→CGG) ‡	yegH	putative transporter
Reads supporting (aligned to +/- strand):  ref base A (31/34);  new base C (3/3);  total (34/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 = 9.06e-01

GGCAATCTCCTCGCCGACCTCGTCAGCCTCGACCCGGCGGCCGCCCAGCAGGCGCAGCACGGCATGGGCAGTACGTTCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGAAGCCCAGGCACAGCATGATCACCGTGGGGTGGGCGTTGACGAAGCGGGTCAG  >  NC_002947/3894243‑3894536                                                                                                                                                     |                                                                                                                                                  ggCAATCTCCTCGCCGACCTCGTCAGCCTCGACCCGGCGGCCGCCCAGCAGGCGCAGCACGGCATGGGCAGTACGTTCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCa                                                                                                                                                   <  3:180363/149‑1 (MQ=255)  gCAATCTCCTCGCCGACCTCGTCAGCCTCGACCCGGCGGCCGCCCAGCAGGCGCAGCACGGCATGGGCAGTACGTTCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAg                                                                                                                                                  <  7:115779/149‑1 (MQ=255)  gCAATCTCCTCGCCGACCTCGTCAGCCTCGACCCGGCGGCCGCCCAGCAGGCGCAGCACGGCATGGGCAGTACGTTCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAg                                                                                                                                                  >  5:153899/1‑149 (MQ=255)   cAATCTCCTCGCCGACCTCGTCAGCCTCGACCCGGCGGCCGCCCAGCAGGCGCAGCACGGCATGGGCAGTACGTTCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAg                                                                                                                                                  >  6:176773/1‑148 (MQ=255)         cctcGCCGACCTCGTCAGCCTCGACCCGGCGGCCGCCCAGCAGGCGCAGCACGGCATGGGCAGTACGTTCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCg                                                                                                                                             <  2:238116/147‑1 (MQ=255)         cctcGCCGACCTCGTCAGCCTCGACCCGGCGGCCGCCCAGCAGGCCCAGCACGGCATGGGCAGTACGTTCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCg                                                                                                                                             >  1:238116/1‑147 (MQ=255)          ctcGCCGACCTCGTCAGCCTCGACCCGGCGGCCGCCCAGCAGGCGCAGCACGGCATGGGCAGTACGTTCGCGCAACGGCTTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCgag                                                                                                                                           <  2:187015/148‑1 (MQ=255)          ctcGCCGACCTCGTCAGCCTCGACCCGGCGGCCGCCCAGCAGGCGCAGCACGGCATGGGCAGTACGTTCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCCTGATc                                                                                                                                              <  5:15071/145‑1 (MQ=255)          ctcGCCGACCTCGTCAGCCTCGACCCGGCGGCCGCCCAGCAGGCGCAGCACGGCATGGGCAGTACGTTCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATc                                                                                                                                              <  3:210652/145‑1 (MQ=255)             gCCGACCTCGTCAGCCTCGACCCGGCGGCCGCCCAGCAGGCGCAGCACGGCATGGGCAGTACGTTCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCCTGATCgagctgt                                                                                                                                       <  4:277314/149‑5 (MQ=255)             gCCGACCTCGTCAGCCTCGACCCGGCGGCCGCCCAGCAGGCGCAGCACGGCATGGGCAGTACGTTCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGaa                                                                                                                                         <  3:250661/147‑1 (MQ=255)              ccGACCTCGTCAGCCTCGACCCGGCGGCCGCCCAGCAGGCGCAGCACGGCATGGGCAGTACGTTCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGcc                                                                                                                                      <  5:35515/149‑1 (MQ=255)                  ccTCGTCAGCCTCGACCCGGCGGCCGCCCAGCAGGCGCAGCACGGCATGGGCAGTACGTTCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATg                                                                                                                                  <  3:328276/149‑1 (MQ=255)                  ccTCGTCAGCCTCGACCCGGCGGCCGCCCAGCAGGCGCAGCACGGCATGGGCAGTACGTTCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATg                                                                                                                                  <  5:176773/149‑1 (MQ=255)                   cTCGTCAGCCTCGACCCGGCGGCCGCCCAGCAGGCGCAGCACGGCATGGGCAGTACGTTCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCCCGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATgg                                                                                                                                 >  5:192986/1‑149 (MQ=255)                      gTCAGCCTCGACCCGGCGGCCGCCCAGCAGGCGCAGCACGGCATGGGCAGTACGTTCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGcc                                                                                                                               <  8:194327/148‑1 (MQ=255)                           ccTCGACCCGGCGGCCGCCCAGCAGGCGCAGCACGGCATGGGCAGTACGTTCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTAc                                                                                                                         <  6:153899/149‑1 (MQ=255)                            cTCGACCCGGCGGCCGCCCAGCAGGCGCAGCACGGCATGGGCAGTACGTTCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTa                                                                                                                          >  6:211803/1‑147 (MQ=255)                               gACCCGGCGGCCGCCCAGCAGGCGCAGCACGGCATGGGCAGTACGTTCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAg                                                                                                                       <  1:308668/147‑1 (MQ=255)                                  ccGGCGGCCGCCCAGCAGGCGCAGCACGGCATGGGCAGTACGTTCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGc                                                                                                                  >  7:85857/1‑149 (MQ=255)                                        gccgccCAGCAGGCGCAGCACGGCATGGGCAGTACGTTCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTgg                                                                                                            <  7:110258/149‑1 (MQ=255)                                         ccgccCAGCAGGCGCAGCACGGCATGGGCAGTACGTTCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTggg                                                                                                           <  2:63625/149‑1 (MQ=255)                                          cgccCAGCAGGCGCAGCACGGCATGGGCAGTACGTTCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTggg                                                                                                           >  2:111476/1‑148 (MQ=255)                                          cgccCAGCAGGCGCAGCACGGCATGGGCAGTACGTTCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTggg                                                                                                           >  4:296614/1‑148 (MQ=255)                                           gccCAGCAGGCGCAGCACGGCATGGGAAGTACGTTCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACCGGCAAGATGGTTGATCAGCTCGATAAGGATCGAGAAGCCGAGGGCCGCGTACAGGTAGCCTTTGGGGa                                                                                                         >  6:313889/1‑149 (MQ=255)                                                             ggCATGGGCAGTACGTTCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCtt                                                                                       >  4:254862/1‑149 (MQ=255)                                                                 tGGGCAGTACGTTCGCGCAACGGCCCGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTcggc                                                                                   >  1:105607/1‑149 (MQ=255)                                                                   ggCAGTACGTTCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTcggcgg                                                                                 <  3:296614/149‑1 (MQ=255)                                                                   ggCAGTACGTTCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTcggcgg                                                                                 <  3:324469/149‑1 (MQ=255)                                                                          cGTTCGCGCAACGGCCCGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCt                                                                          <  2:105607/149‑1 (MQ=255)                                                                           gTTCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCt                                                                          >  4:69083/1‑148 (MQ=255)                                                                            ttCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGa                                                                        <  8:85857/149‑1 (MQ=255)                                                                               gcgcAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAAc                                                                     >  4:300691/1‑149 (MQ=255)                                                                               gcgcAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAAc                                                                     >  4:305522/1‑149 (MQ=255)                                                                                  cAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGa                                                                  <  3:22553/149‑1 (MQ=255)                                                                                   aaCGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCACAGGCCTTCGGCGGTGAGGCTGAAACCGa                                                                  >  6:8322/1‑148 (MQ=255)                                                                                     cGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGatca                                                               >  2:308195/1‑149 (MQ=255)                                                                                          tgtgCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGaaa                                                                      >  1:339482/1‑137 (MQ=255)                                                                                          tgtgCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGaaa                                                                      <  2:339482/137‑1 (MQ=255)                                                                                            tgCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGatcatcatc                                                          <  8:196798/147‑1 (MQ=255)                                                                                            tgCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGAtcatcatca                                                         <  5:211803/148‑1 (MQ=255)                                                                                                 ttGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGaa                                                     <  1:308195/147‑1 (MQ=255)                                                                                                 ttGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGAAGc                                                   >  1:256606/1‑149 (MQ=255)                                                                                                   gCAGGCTGCGCTTGCGGCGGGCATGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGAAGcc                                                  >  2:322607/1‑148 (MQ=255)                                                                                                   gCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGAAGcc                                                  >  7:260693/1‑148 (MQ=255)                                                                                                    cAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGAAGccc                                                 <  8:260693/148‑1 (MQ=255)                                                                                                      ggCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGAAGCCCAgg                                              <  1:111476/149‑1 (MQ=255)                                                                                                                gcggcgGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGAAGCCCAGGCACAGCATg                                     <  3:69083/148‑1 (MQ=255)                                                                                                                    cggGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGAAGCCCAGGCACAGCATGATCAc                                <  3:300691/149‑1 (MQ=255)                                                                                                                    cggGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGAAGCCCAGGCACAGCATGATCAc                                <  3:305522/149‑1 (MQ=255)                                                                                                                       gCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGAAGCCCAgg                                              <  4:313178/132‑1 (MQ=255)                                                                                                                       gCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGAAGCCCAgg                                              >  3:313178/1‑132 (MQ=255)                                                                                                                        catgggCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGAAGCCCAGGCACAGCATGATCACCgtg                            <  1:322607/146‑1 (MQ=255)                                                                                                                        cACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGAAGCCCAGGCACAGCATGATCACCgtg                            >  3:108249/1‑149 (MQ=255)                                                                                                                            ggCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGAAGCCCAGGCACAGCATGATCACCgtggggt                        <  7:154507/149‑1 (MQ=255)                                                                                                                               cAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGAAGCCCAGGCACAGCATGATCACCgtggggtggg                     >  3:318310/1‑149 (MQ=255)                                                                                                                                aGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGATGCCCAGGCACAGCATGATCACCGTGGGGTGGGc                    >  8:111149/1‑149 (MQ=255)                                                                                                                                  cTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGAAGCCCAGGCACAGCATGATCACCGTGGGGTGGGCGt                  <  6:192986/149‑1 (MQ=255)                                                                                                                                     gTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGACACCGATTATAATCAGGAAGCCCAGGCACAGCAGGCTCACCGTGTGGTGGGCGtggc               >  4:52582/1‑146 (MQ=255)                                                                                                                                     gTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGACGCCCAGGTACAGCATGATCACCGTGGGGTGGGCGTTGa               >  6:4805/1‑149 (MQ=255)                                                                                                                                     gTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGAAGCCCAGGCACAGCATGATCACCGTGGGGTGGGCGTTGa               >  8:220042/1‑149 (MQ=255)                                                                                                                                       tcggcttcTCGATCCTGATCGAGCAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGAAGCCCAGGCACAGCATGATCACCGTGGGGTGGGCGTTGAc              >  6:57596/8‑148 (MQ=255)                                                                                                                                         aaCAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGGGAGGCTGACACCGATCATCCTCAGGAAGCCCAGGCAGAGCATGCTCACCCTGGGGTGGGCGTGGACGac           >  8:138198/1‑148 (MQ=255)                                                                                                                                           cttcTCGATCCTGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGAAGCCCAGGCACAGCATGATCACCGTGGGGTGGGCGTTGACGAAGc         >  5:257183/4‑149 (MQ=255)                                                                                                                                            ttcTCGATCCTGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGAAGCCCAGGCACAGCATGATCACCgtggggtg                       >  1:14796/3‑134 (MQ=255)                                                                                                                                            ttcTCGATCCTGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGAAGCCCAGGCACAGCATGATCACCgtggggtg                       <  2:14796/132‑1 (MQ=255)                                                                                                                                             gCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGAAGCCCAGGCACAGCATGATCACCGTGGGGTGGGCGTTGACGAAGCgg       <  4:318310/149‑1 (MQ=255)                                                                                                                                              cTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGAAGCCCAGGCACAGCATGATCACCGTGGGGTGGGCGTTGACGAAGCggg      >  1:315681/1‑149 (MQ=255)                                                                                                                                                 gATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGAAGCCCAGGCACAGCATGATCACCGTGGGGTGGGCGTTGACGAAGCGGGt     >  8:103769/1‑147 (MQ=255)                                                                                                                                                 gATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGAAGCCCA‑GCACAGCATGATCACCgtggg                          >  5:98534/1‑125 (MQ=255)                                                                                                                                                 gATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGAAGCCCA‑GCACAGCATGATCACCgtggg                          <  6:98534/125‑1 (MQ=255)                                                                                                                                                  aTCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGAAGCCCAGGCACAGCATGATCACCGTGGGGTGGGCGTTGACGAAGCGGGTCAg  <  7:220042/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GGCAATCTCCTCGCCGACCTCGTCAGCCTCGACCCGGCGGCCGCCCAGCAGGCGCAGCACGGCATGGGCAGTACGTTCGCGCAACGGCCTGTGCTGTTGCAGGCTGCGCTTGCGGCGGGCACGGGCCAGCTGGTTGAACAGCTCGATCAGGATCGAGAAGCCGATGGCCGCGTACAGGTAGCCTTTGGGGATATGGAAGCCCAGGCCTTCGGCGGTGAGGCTGAAACCGATCATCATCAGGAAGCCCAGGCACAGCATGATCACCGTGGGGTGGGCGTTGACGAAGCGGGTCAG  >  NC_002947/3894243‑3894536

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