Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,170,180	T→G	12.4%	intergenic (‑288/‑99)	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,180	0	T	G	12.4%	128.2 / 12.3	64	intergenic (‑288/‑99)	katG/PP_3669	catalase‑peroxidase/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (27/27);  new base G (3/5);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.62e-01

CACGCAGGTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTG  >  NC_002947/4170032‑4170323                                                                                                                                                     |                                                                                                                                                 cACGCAGGTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGt                                                                                                                                                  <  5:166012/149‑1 (MQ=255)  aCGCAGGTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGt                                                                                                                                                  >  5:187828/1‑148 (MQ=255)    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                                                                                                                                   <  8:15661/149‑1 (MQ=255)                gcgcTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGAg                                                                                                                                   <  6:320797/149‑1 (MQ=255)                gcgcTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGAg                                                                                                                                   <  5:106443/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                                                                                                                                       >  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                         <  5:325575/149‑1 (MQ=255)                                                                                                                         gcgcAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGcgcc                         >  7:68383/1‑149 (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)                                                                                                                                                     |                                                                                                                                                 CACGCAGGTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGG‑CACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTG  >  NC_002947/4170032‑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: ‑