Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,170,183	A→T	9.6%	intergenic (‑291/‑96)	katG ← / → PP_3669	catalase‑peroxidase/LysR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,170,183	0	A	T	9.6%	147.9 / 8.2	63	intergenic (‑291/‑96)	katG/PP_3669	catalase‑peroxidase/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (28/26);  new base T (2/4);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.44e-02

GCAGGTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTG  >  NC_002947/4170035‑4170323                                                                                                                                                     |                                                                                                                                              gCAGGTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCaa                                                                                                                                               <  1:86282/149‑1 (MQ=255)     gTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGCCAAGccc                                                                                                                                           >  3:236627/1‑149 (MQ=255)         tGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCACGCCCGCGa                                                                                                                                       >  1:148177/1‑149 (MQ=255)          gACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGAt                                                                                                                                      >  4:317950/1‑149 (MQ=255)           aCGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGAt                                                                                                                                      <  3:317950/148‑1 (MQ=255)             gcgcTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGAg                                                                                                                                   <  5:106443/149‑1 (MQ=255)             gcgcTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGAg                                                                                                                                   <  8:15661/149‑1 (MQ=255)             gcgcTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGAg                                                                                                                                   <  6:320797/149‑1 (MQ=255)               gcTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGAGGc                                                                                                                                 <  4:239054/149‑1 (MQ=255)                 tGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCg                                                                                                                               >  5:165156/1‑149 (MQ=255)                   cccTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑c                                                                                                                            <  6:187828/149‑1 (MQ=255)                             gAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATac                                                                                                                  <  6:291229/149‑1 (MQ=255)                              aaTTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATAGCGGTCAAGCCCGCGATGAGGCCGG‑cacaggcattc                                                                                                                  >  5:26127/1‑141 (MQ=255)                               ttgcgcgGGGCCACGATAAAGATTCGCCGCGGTTGGGCATTGGTCTATTAGAGTGGGTGTGGGCGCGTTATCGGGGGGCTCGGTCAGGCGCCGGCGTGTATCGGCCTCATCGCGGTCACGCCCGCGATGAGGCCGG‑CACATCCATacaa                                                                                                                <  2:151832/146‑1 (MQ=255)                                   gcgGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATacaaggca                                                                                                            >  5:64569/1‑149 (MQ=255)                                          cACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCACATAGCGGTGAACCCCCCGATGAGGGCCG‑CCCATCACAAcaaagcccaaagcc                                                                                                     >  1:183129/1‑149 (MQ=255)                                               gaAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGAt                                                                                                                                      >  1:361050/1‑112 (MQ=255)                                               gaAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGAt                                                                                                                                      <  2:361050/112‑1 (MQ=255)                                                     cTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGc                                                                                          <  5:200411/149‑1 (MQ=255)                                                          cGAGCTTGGGCAGAGGTCTAATTGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGCCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCaa                                                                                     <  4:236627/149‑1 (MQ=255)                                                          cGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGACACAAGGCCCCGCGCCATGCTTCaa                                                                                     <  2:41740/149‑1 (MQ=255)                                                          cGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATAGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCaa                                                                                     >  2:281154/1‑149 (MQ=255)                                                            aGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATa                                                                                   <  4:46942/149‑1 (MQ=255)                                                              cTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATaca                                                                                 >  8:235470/1‑149 (MQ=255)                                                               ttGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCATAAGGCCCCGCGCCATGCTTCAATaca                                                                                 >  3:342591/1‑148 (MQ=255)                                                                 gggCAGAGGTCTAATAGAGTGGCTGTGGGCGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCACGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATAcaca                                                                               <  4:143289/148‑1 (MQ=255)                                                                   gCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACCTCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGc                                                                             >  3:359799/1‑148 (MQ=255)                                                                   gCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGc                                                                             >  3:365501/1‑148 (MQ=255)                                                                         gTCTAATAGAGTGGCTGTTGTAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATAc                                                                        <  1:281154/148‑1 (MQ=255)                                                                            tAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCttt                                                                   >  2:74171/1‑149 (MQ=255)                                                                            tAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCttt                                                                   >  2:130098/1‑149 (MQ=255)                                                                              aTAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATAGCGGTAAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACCAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTga                                                                 >  1:334879/1‑148 (MQ=255)                                                                                 gagTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCa                                                              <  6:165156/149‑1 (MQ=255)                                                                                            ggAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGc                                                   >  6:325575/1‑149 (MQ=255)                                                                                                ccTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTc                                               <  6:64569/149‑1 (MQ=255)                                                                                                  tGATAGGCTGGGTCTGTCAGGCGCAGGCCTGTATCGGCCTCATAGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCAt                                             >  8:243779/1‑149 (MQ=255)                                                                                                         ctggctCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAATCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCg                                       <  6:200779/148‑1 (MQ=255)                                                                                                         ctggctCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGa                                                                                                                                       <  1:168066/53‑1 (MQ=255)                                                                                                         ctggctCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGa                                                                                                                                       >  2:168066/1‑53 (MQ=255)                                                                                                         ctggctCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGc                                      >  4:75116/1‑149 (MQ=255)                                                                                                           tgctCTGTCAGGCGTAGGCCTGTATCGGCCTCATGGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGGTTCAATAGACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCaa                                    <  4:143676/148‑1 (MQ=255)                                                                                                            gctCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAg                                                                              >  7:212571/1‑106 (MQ=255)                                                                                                            gctCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAg                                                                              <  8:212571/106‑1 (MQ=255)                                                                                                            gctCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCaa                                    <  3:75116/148‑1 (MQ=255)                                                                                                             ctctGTCAGGCGCAGGCCTGTATCGGCCTCAAAGCGGGCAAACCCGCGATGAGGCCGG‑CCCCAACAAACAAGGCACAAGGCCCCCCGCCCAGGGTAAAAAAACCCCCAAAACCGTTGCCCCGGACCCCCCGCCCGCAGGAACCGCAAcc                                  >  2:362948/1‑148 (MQ=255)                                                                                                                      gcgcAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGcgcc                         >  7:68383/1‑149 (MQ=255)                                                                                                                      gcgcAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGcgcc                         <  5:325575/149‑1 (MQ=255)                                                                                                                              cTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGGGGCCGG‑CCCATGCATACAAGGCCCAAGGCCCCGCGCCCTGCGTGAATACACAGGCAAACCCTTTGCCCAGGGCGCCGAGCCCGCCTGAACCGCCACGCCCGTCGCcgcggggg                  >  4:39798/1‑144 (MQ=255)                                                                                                                              cTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACCTCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGAc                 <  7:260328/149‑1 (MQ=255)                                                                                                                              cTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGa                  >  6:143827/1‑148 (MQ=255)                                                                                                                                   tCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAAc                                   <  3:202690/126‑1 (MQ=255)                                                                                                                                     ggCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCATAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAAc           <  4:342591/148‑1 (MQ=255)                                                                                                                                        cTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACctgc       <  7:243779/149‑1 (MQ=255)                                                                                                                                         tCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACctgct      >  2:35854/1‑149 (MQ=255)                                                                                                                                         tCATCGCGGGCTTGACCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACctgct      >  1:243261/1‑149 (MQ=255)                                                                                                                                           aTCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGa                  <  7:154720/135‑1 (MQ=255)                                                                                                                                           aTCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGa                  >  8:154720/1‑135 (MQ=255)                                                                                                                                           aTCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACctgctg     <  7:251013/148‑1 (MQ=255)                                                                                                                                           aTCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACctgct      <  8:68383/147‑1 (MQ=255)                                                                                                                                           aTCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTgg    <  5:143827/149‑1 (MQ=255)                                                                                                                                           aTCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCACGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTgg    >  4:35774/1‑149 (MQ=255)                                                                                                                                           aTCGCGGGCTTGACCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTgg    >  6:88868/1‑149 (MQ=255)                                                                                                                                             cgcgGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCATCATGAACCGCAACGACCTGCGACGCGTCGACATGAACctgctgctg  >  5:235757/1‑149 (MQ=255)                                                                                                                                              gcgGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTg  <  7:235470/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                              GCAGGTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTG  >  NC_002947/4170035‑4170323

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