Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,242,528	G→T	9.4%	intergenic (‑62/+138)	bfr‑II ← / ← PP_1083	bacterioferritin/(2Fe‑2S)‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,242,528	0	G	T	9.4%	181.1 / 6.7	61	intergenic (‑62/+138)	bfr‑II/PP_1083	bacterioferritin/(2Fe‑2S)‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (34/21);  new base T (2/4);  total (36/25)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GTGCAGGAAGTACTGGTTGATCGCGACCAGCTCGTTTCCGAGGATCTTGTTGAGATGCTGGATGACGCTTACGTCGCCTTTCATGTCGGGTTCCTGCCCTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGAATTGCGGGCATAAAAAAACCGGACACCAGGTCCGGTTCTTGAAATTCTGCGTTCAGGCAGCGG  >  NC_002947/1242383‑1242676                                                                                                                                                  |                                                                                                                                                     gtgcAGGAAGTACTGGTTGATCGCGACCAGCTCGTTTCCGAGGATCTTGTTGAGATGCTGGATGACGCTTACGTCGCCTTTCATGTCGGGTTCCTGCCCTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGgtgtgt                                                                                                                                                   <  5:110091/149‑1 (MQ=255)          gTACTGGTTGATCGCGACCAGCTCGTTTCCGAGGATCTTGTTGAGATGCTGGATGACGCTTACGTCGCCTTTCATGTCGGGTTCCTGCCCTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTaa                                                                                                                                          <  4:194208/149‑1 (MQ=255)             cTGGTTGATCGCGACCAGCTCGTTTCCGAGGATCTTGTTGAGATGCTGGATGACGCTTACGTCGCCTTTCATGTCGGGTTCCTGCCCTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGt                                                                                                                                        >  8:129350/1‑148 (MQ=255)                gTTGATCGCGACCAGCTCGTTTCGGTGGATCTTTTTGAGATTCTGGATGACGCTTACGTCGCCTTTCATGTCGGTTTCCTGCCCTGTAGAATTGTGTCAATCCACTGAATTTTGAGCCCGATACTTAGGTTTGACACACCTAAGtatc                                                                                                                                     <  8:350691/148‑4 (MQ=255)                gTTGATCGCGACCAGCTCGTTTCCGAGGATCTTGTTGAGATGCTGGATGACGCTTACGTCGCCTTTCATGTCGGGTTCCTGCCCTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCtt                                                                                                                                    <  3:330923/149‑1 (MQ=255)                gTTGATCGCGACCAGCTCGTTTCCGAGGATCTTGTTGAGATGCTGGATGACGCTTACGTCGCCTTTCATGTCGGGTTCCTGCCCTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCtt                                                                                                                                    <  7:250334/149‑1 (MQ=255)                 ttGATCGCGACCAGCTCGTTTCCGAGGATCTTGTTGAGATGCTGGATGACGCTTACGTCGCCTTTCATGTCGGGTTCCTGCCCTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTTTGACACACCTAAGtatcgg                                                                                                                                   <  7:136559/149‑6 (MQ=255)                    aTCGCGACCAGCTCGTTTCCGAGGATCTTGTTGAGATGCTGGATGACGCTTACGTCGCCTTTCATGTCGGGTTCCTGCCCTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGaa                                                                                                                                 >  6:119826/1‑148 (MQ=255)                                tCGTTTCCGAGGATCTTGTTGAGATGCTGGATGACGCTTACGTCGCCTTTCATGTCGGGTTCCTGCCCTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGaaaa                                                                                                                    >  8:357621/1‑149 (MQ=255)                                      ccGAGGATCTTGTTGAGATGCTGGATGACGCTTACGTCGCCTTTCATGTCGGGTTCCTGCCCTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAACCTTaa                                                                                                                                          >  8:59713/1‑121 (MQ=255)                                      ccGAGGATCTTGTTGAGATGCTGGATGACGCTTACGTCGCCTTTCATGTCGGGTTCCTGCCCTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAAt                                                                                                               <  7:129350/148‑1 (MQ=255)                                       cGAGGATCTTGTTGAGATGCTGGATGACGCTTACGTCGCCTTTCATGTCGGGTTCCTGCCCTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAg                                                                                                             >  1:354709/1‑149 (MQ=255)                                        gAGGATCTTGTTGAGATGCTGGATGACGCTTACGTCGCCTTTCATGTCGGGTTCCTGCCCTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAg                                                                                                             >  4:155004/1‑148 (MQ=255)                                                 gttgAGATGCTGGATGACGCTTACGTCGCCTTTCATGTCGGGTTCCTGCCCTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAAt                                                                                                   >  1:164680/1‑149 (MQ=255)                                                   agagaTGCTGGATGACGCTTACGTCTCCTTTCATGTCGGGGTCCTGCCCTCTAGAAGTGTGTCAAACCACAGAAGAATGAGCCCGAAACATAGGTGAGTCAAACCTAAGTTCATGAATACTAAGTGAAAATTAAAAGGAATAAGAATgg                                                                                                 >  8:297696/2‑148 (MQ=255)                                                     agaTGCTGGATGACGCTTACGTCGCCTTTCATGTCGGGTTCCTGCCCTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGttt                                                                                               >  4:74524/1‑149 (MQ=255)                                                      gaTTCTTGATGAGGGTTACGTCGCATTTCATGTCGGGTTCATGCCCTGTAGAAGTGTTTCAATACACTGAAGTTTGAGCCCGATACTTAGGTTTTTCAAACCTAAGTTCTTTAATACTAAGTGAAAATTAATAGGAATAAGAATGttt                                                                                               <  7:297696/148‑1 (MQ=255)                                                           tGGATGACGCTTACGTCGCCTTTCATGTCGGGTTCCTGCCCTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATTCTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAAc                                                                                         >  5:26344/1‑149 (MQ=255)                                                                     ttACGTCGCCTTTCATGTCGGGTTCCTGCCCTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCaa                                                                               >  8:142559/1‑149 (MQ=255)                                                                        cGTCGCCTTTCATGTCGGGTTCCTGCCCTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGc                                                                            <  7:357621/149‑1 (MQ=255)                                                                            gCCTTTCATGTCGGGTTCCTGCCCTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTaa                                                                         <  6:26344/148‑1 (MQ=255)                                                                                ttCATGTCGGGTTCCTGCCCTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTaa                                                                                                                           >  4:283643/1‑94 (MQ=255)                                                                                ttCATGTCGGGTTCCTGCCCTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTaa                                                                                                                           <  3:283643/94‑1 (MQ=255)                                                                                   aTGTCGGGTTCCTGCCCTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACttgttg                                                                  <  3:155004/148‑1 (MQ=255)                                                                                    tGTCGGGTTCCTGCCCTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGaa                                                                >  1:152391/1‑149 (MQ=255)                                                                                      tCGGGTTCCTGCCCTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGAAtt                                                              >  2:286298/1‑149 (MQ=255)                                                                                        gggTTCCTGCCCTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGAATTGc                                                            >  8:93143/1‑149 (MQ=255)                                                                                          gTTCCTGCCCTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGAATTGCgg                                                          <  7:142559/149‑1 (MQ=255)                                                                                           ttCCTGCCCTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGAATTGCggg                                                         <  5:119826/149‑1 (MQ=255)                                                                                              cTGCCCTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACttgtt                                                                   >  2:245299/1‑136 (MQ=255)                                                                                              cTGCCCTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACttgtt                                                                   <  1:245299/136‑1 (MQ=255)                                                                                                 cccTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGc                                                                            >  7:13599/1‑124 (MQ=255)                                                                                                 cccTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGc                                                                            <  8:13599/124‑1 (MQ=255)                                                                                                  ccTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGAATTGCg                                                           <  3:243912/140‑1 (MQ=255)                                                                                                  ccTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGAATTGCg                                                           >  4:243912/1‑140 (MQ=255)                                                                                                  ccTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGAATTGCGGGCATaaaa                                                  <  7:93143/149‑1 (MQ=255)                                                                                                     gTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGAATTGCGGGCATaaaaa                                                 <  4:32921/147‑1 (MQ=255)                                                                                                     gTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGAATTGCGGGCATaaaaa                                                 >  3:32921/1‑147 (MQ=255)                                                                                                        gAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGtt                                                                                   <  2:164787/110‑1 (MQ=255)                                                                                                        gAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGtt                                                                                   >  1:164787/1‑110 (MQ=255)                                                                                                           gtgtgtCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGAATTGCGGGCATAAAAAAACCGGa                                          >  5:54998/1‑148 (MQ=255)                                                                                                              tgtCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGAATTGCGGGCATAAAAAAACCGGacctg                                      >  2:64060/1‑146 (MQ=255)                                                                                                              tgtCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGAATTGCGGGCATAAAAAAACCGGacac                                       >  1:358674/1‑148 (MQ=255)                                                                                                                  gaTCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGAATTGCGGGCATAAAAAAACCGGACACCAgga                                  >  2:227714/2‑148 (MQ=255)                                                                                                                  aaTCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGa                                                                                                     >  1:280702/1‑82 (MQ=255)                                                                                                                  aaTCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGa                                                                                                     <  2:280702/82‑1 (MQ=255)                                                                                                                   aTCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGAATTGCGGGCATAAAAAAACCGGACACCAGGTc                                 >  1:119610/1‑149 (MQ=255)                                                                                                                   aTCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGAATTGCGGGCATAAAAAAACCGGACACCAGGTc                                 >  6:349283/1‑149 (MQ=255)                                                                                                                   aTCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGAATTGCGGGCATAAAAAAACCGGACACCAGGTc                                 >  3:220797/1‑149 (MQ=255)                                                                                                                      cACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGAATTGCGGGCATAAAAAAACCGGACACCAGGTCCg                               <  3:74524/148‑1 (MQ=255)                                                                                                                       aCTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGAATTGCGGGCATAAAAAAACCGGACACCAGGTCCGGt                             >  7:215634/1‑149 (MQ=255)                                                                                                                       aCTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGAATTGCGGGCATAAAAAAACCGGACACCAGGTCCGGt                             >  3:57540/1‑149 (MQ=255)                                                                                                                         tGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAAt                                                                                                               >  4:137124/1‑65 (MQ=255)                                                                                                                         tGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAAt                                                                                                               <  3:137124/65‑1 (MQ=255)                                                                                                                         tGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGAATTGCGGGCATAAAAAAACCGGACACCAGGTCCGGTTc                           >  5:191033/1‑149 (MQ=255)                                                                                                                                   ttcaagaaCTTAGGTTTGACACACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGAATTGCGGGCATAAAAAAACCGGACACCAGGTCCGGTTCTTGAAATTCt                 >  3:343986/8‑149 (MQ=255)                                                                                                                                   gCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGAATTGCGGGCATAAAAAAACCGGACACCAGGTCCGGTTCTTGAAATTCt                 <  6:191033/149‑1 (MQ=255)                                                                                                                                       agaaCTTAGGTTTGACACACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGAATTGCGGGCATAAAAAAACCGGACACCAg                                    <  5:190796/123‑1 (MQ=255)                                                                                                                                       agaaCTTAGGTTTGACACACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGAATTGCGGGCATAAAAAAACCGGACACCAg                                    >  6:190796/4‑126 (MQ=255)                                                                                                                                         tACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGAATTGCGGGCATAAAAAAACCGGACACCAGGTCCGGTTCTTGAAATTCCGCGtt            >  6:143945/1‑148 (MQ=255)                                                                                                                                            ttAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGt                                                                                    >  7:314171/1‑73 (MQ=255)                                                                                                                                            ttAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGt                                                                                    <  8:314171/73‑1 (MQ=255)                                                                                                                                                gtgtgtCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGGAACGCTAACTTGTTGAATTGCGGGCATAAAAAAACCGGACACCAGGTCCGGTTCTTGAAATTCTGCGTTCAGGCAGc    <  2:152391/149‑1 (MQ=255)                                                                                                                                                gtgtgtCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGAATTGCGGGCATAAAAAAACCGGACACCAGGTCCGGTTCTTGAAATTCTGCGTTCAGGCAg     <  8:215634/148‑1 (MQ=255)                                                                                                                                                 tgtgtCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGAATTGCGGGCATAAAAAAACCGGACACCAGGTCCGGTTCTTGAAATTCTGCGTTCAGGCAGc    <  2:358674/148‑1 (MQ=255)                                                                                                                                                  gtgtCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGAATTGCGGGCATAAAAAAACCGGACACCAGGTCCGGTTCTTGAAATTCTGCGTTCAGGCAGCgg  >  3:162470/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                     GTGCAGGAAGTACTGGTTGATCGCGACCAGCTCGTTTCCGAGGATCTTGTTGAGATGCTGGATGACGCTTACGTCGCCTTTCATGTCGGGTTCCTGCCCTGTAGAAGTGTGTCAATCCACTGAAGTTTGAGCCCGATACTTAGGTGTGTCAAACCTAAGTTCTTGAATACTAAGTGAAAATTAATAGGAATAAGAATGTTTGTGAACCGCGTTGCAACGCTAACTTGTTGAATTGCGGGCATAAAAAAACCGGACACCAGGTCCGGTTCTTGAAATTCTGCGTTCAGGCAGCGG  >  NC_002947/1242383‑1242676

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