Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	344,045:1	+T	9.2%	intergenic (+134/+55)	PP_5440 → / ← PP_t72	hypothetical protein/tRNA‑Phe

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	344,045	1	.	T	9.2%	245.1 / 18.8	65	intergenic (+134/+55)	PP_5440/PP_t72	hypothetical protein/tRNA‑Phe
Reads supporting (aligned to +/- strand):  ref base . (30/29);  new base T (2/4);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

GAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTAAA  >  NC_002947/343906‑344191                                                                                                                                              |                                                                                                                                                   gAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGC‑CTAGTGCGGGtttt                                                                                                                                           >  4:144534/1‑148 (MQ=255) gAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTc                                                                                                                                           >  5:329467/1‑149 (MQ=255)   gTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCtt                                                                                                                                         <  7:111107/149‑1 (MQ=255)      aGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGC‑CTAGTGCGGGTTTTttgtt                                                                                                                                      <  1:96606/149‑1 (MQ=39)      aGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGC‑CTAGTGCGGGTTTTttgtt                                                                                                                                      <  3:144534/149‑1 (MQ=39)        gTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTtgttgt                                                                                                                                    >  4:940/1‑149 (MQ=255)         tttGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCgg                                                                                                                                                >  8:135738/1‑136 (MQ=255)         tttGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCgg                                                                                                                                                <  7:135738/136‑1 (MQ=255)            ggCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGC‑CTAGTGCGGGTTTTTTGTTGTctg                                                                                                                                 <  7:1410/148‑3 (MQ=38)              cccGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCGTc                                                                                                                              >  8:247903/1‑149 (MQ=255)               ccGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCg                                                                                                                             >  4:297863/1‑149 (MQ=255)                     tCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGAATGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTa                                                                                                                         >  1:247350/1‑148 (MQ=255)                         cTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACttt                                                                                                                    <  3:940/148‑1 (MQ=255)                         cTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACtt                                                                                                                     <  3:57126/147‑1 (MQ=255)                          tCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCttgtt                                                                                                                                      <  1:169332/129‑1 (MQ=255)                          tCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCttgtt                                                                                                                                      >  2:169332/1‑129 (MQ=255)                          tCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGc                                                                                                                  <  6:6945/149‑1 (MQ=255)                                       tttCTCGATTGCATGGCTAGCCTTCAGGTTGCG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGc                                                                                                     >  6:255996/1‑149 (MQ=255)                                           tcGATTGCATGGCTAGCCTTCAGGTTGCG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCgg                                                                                                                            >  7:224247/1‑122 (MQ=255)                                           tcGATTGCATGGCTAGCCTTCAGGTTGCG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCgg                                                                                                                            <  8:224247/122‑1 (MQ=255)                                           tcGATTGCATGGCTAGCCTTCAGGTTGCG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTg                                                                                                  >  6:298395/1‑148 (MQ=255)                                             gATTGCATGGCTAGCCTTCAGGTTGCG‑TTGTGATGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGaa                                                                                                <  8:91807/148‑1 (MQ=255)                                             gATTGCATGGCTAGCCTTCAGGTTGCG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAg                                                                                               >  7:271006/1‑149 (MQ=255)                                             gATTGCATGGCTAGCCTTCAGGTTGCG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAg                                                                                               <  5:255996/149‑1 (MQ=255)                                              aTTGCATGGCTAGCCTTCAGGTTGCG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGc                                                                                                                                 >  6:121262/1‑114 (MQ=255)                                              aTTGCATGGCTAGCCTTCAGGTTGCG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGc                                                                                                                                 <  5:121262/114‑1 (MQ=255)                                                    tGGCTAGCCTTCAGGTTGCG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGc                                                                                        <  1:112615/149‑1 (MQ=255)                                                            cTTCAGGTTGCG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGAc                                                                                <  2:294469/149‑1 (MQ=255)                                                             ttCAGGTTGCG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTtgttgt                                                                                                                                    <  8:171567/96‑1 (MQ=255)                                                             ttCAGGTTGCG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTtgttgt                                                                                                                                    >  7:171567/1‑96 (MQ=255)                                                               cAGGTTGCG‑TTGTGAGGTCTAGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACgg                                                                              <  8:271006/148‑1 (MQ=255)                                                                      cG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCGTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGaa                                                                       <  7:247903/148‑1 (MQ=255)                                                                       g‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCa                                                                                     <  6:101294/133‑1 (MQ=255)                                                                       g‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCa                                                                                     >  5:101294/1‑133 (MQ=255)                                                                              aGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACttt                                                                                                                    >  4:142022/1‑96 (MQ=255)                                                                               ggTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGaca                                                              >  7:35365/1‑149 (MQ=255)                                                                                    cGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAg                                                          >  1:48530/1‑148 (MQ=255)                                                                                                 gccgccTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCctc                                            <  6:329467/149‑1 (MQ=255)                                                                                                     ccTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGct                                        <  5:298395/149‑1 (MQ=255)                                                                                                       tGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCTACACGAGGATTTTCAATCCTCTGctc                                       >  2:188974/1‑148 (MQ=255)                                                                                                       tGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGctc                                       >  4:133307/1‑148 (MQ=255)                                                                                                        gAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTa                                     >  1:54260/1‑149 (MQ=255)                                                                                                           gTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGAGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTAcc                                   >  4:288333/1‑148 (MQ=255)                                                                                                                aaaTCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGa                             <  2:54529/149‑1 (MQ=255)                                                                                                                   tCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCt                          <  2:54260/149‑1 (MQ=255)                                                                                                                        gACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCgg                                                                                                                            >  6:159336/1‑46 (MQ=255)                                                                                                                        gACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCgg                                                                                                                            <  5:159336/46‑1 (MQ=255)                                                                                                                         acaacaAGAAACCCGC‑CTAGTGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTgg                    >  4:194323/1‑149 (MQ=255)                                                                                                                         acaacaAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACgg                                                                              >  8:17435/1‑91 (MQ=255)                                                                                                                         acaacaAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACgg                                                                              <  7:17435/91‑1 (MQ=255)                                                                                                                         acaacaAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGa                             >  3:313646/1‑140 (MQ=255)                                                                                                                         acaacaAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGa                             <  4:313646/140‑1 (MQ=255)                                                                                                                         acaacaAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTAGCAACCACCAGAAGCTTGAAGCTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTATCAATCCTCTGCTCTACCGACTGAGCTATcagt                    >  1:48599/1‑146 (MQ=255)                                                                                                                           aacaaGAAACCCGC‑CTAGTGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGc                  <  3:194323/149‑1 (MQ=255)                                                                                                                           aacaaAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGc                  <  8:77024/149‑1 (MQ=255)                                                                                                                           aacaaAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGc                  >  6:39684/1‑149 (MQ=255)                                                                                                                             caaAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAAccgcc                                                                  >  8:87601/1‑99 (MQ=255)                                                                                                                             caaAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAAccgcc                                                                  <  7:87601/99‑1 (MQ=255)                                                                                                                             caaAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTATCTTCGGGTAACTTTGAAACAACCAGAAGCTTGAAGGTGGTGACCAGAGACGGAGTCGAAACGCCGAAACGATGATTTTCAATAATCTGATCTAAAGACTGAGCTATCTGGGcg                 >  5:150827/1‑148 (MQ=255)                                                                                                                              aaaaaaCCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGAtttt                                                    >  5:73093/1‑112 (MQ=255)                                                                                                                              aaaaaaCCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGAtttt                                                    <  6:73093/112‑1 (MQ=255)                                                                                                                                aaaaCCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGcgacga             <  2:247350/149‑1 (MQ=255)                                                                                                                                     ccGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGa        <  3:297863/149‑1 (MQ=255)                                                                                                                                       gCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACCTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGcgac               >  3:205853/1‑140 (MQ=255)                                                                                                                                       gCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACCTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGcgac               <  4:205853/140‑1 (MQ=255)                                                                                                                                             g‑gCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTaaa  <  3:133307/147‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                   GAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCG‑TTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTGAATTAAA  >  NC_002947/343906‑344191

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