Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,170,185	C→A	9.5%	intergenic (‑293/‑94)	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,185	0	C	A	9.5%	178.6 / 9.8	63	intergenic (‑293/‑94)	katG/PP_3669	catalase‑peroxidase/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (29/28);  new base A (2/4);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.04e-01

GTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTTCGAAG  >  NC_002947/4170039‑4170333                                                                                                                                                   |                                                                                                                                                      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                                                                                                                                             <  6:320797/149‑1 (MQ=255)         gcgcTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGAg                                                                                                                                             <  8:15661/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                                                                                                                                                <  2:361050/112‑1 (MQ=255)                                           gaAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGAt                                                                                                                                                >  1:361050/1‑112 (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:130098/1‑149 (MQ=255)                                                                        tAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCttt                                                                             >  2:74171/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                                                                                                                                                 >  2:168066/1‑53 (MQ=255)                                                                                                     ctggctCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGa                                                                                                                                                 <  1:168066/53‑1 (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)                                                                                                                                                  gCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGACCCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTTCg     >  2:14691/1‑147 (MQ=255)                                                                                                                                                   cccGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTTCGAAg  >  3:310853/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                      GTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTTCGAAG  >  NC_002947/4170039‑4170333

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