Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,242,308	G→A	10.2%	intergenic (+51/‑120)	PP_t30 → / → PP_1978	tRNA‑Glu/TetR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,242,308	0	G	A	10.2%	159.6 / 7.7	59	intergenic (+51/‑120)	PP_t30/PP_1978	tRNA‑Glu/TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (28/25);  new base A (3/3);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

AGTGGATCTTCGCGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAATATGATGAGAATAATATTTCGATCATCATATTTTGCTTGGTGCGTGCCAGATGAGCGACAAGAAGAGCAGAAC  >  NC_002947/2242160‑2242450                                                                                                                                                         |                                                                                                                                               aGTGGATCTTCGCGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATgg                                                                                                                                                <  6:333198/149‑1 (MQ=255)         ttCGCGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGGGAGCACAAAGGGGCCAGCCACGTTTAACCGCGTTTTTCGGGAATTTAAAAGCTCATTTGCTTATTTAAAggggca                                                                                                                                            >  1:117685/1‑146 (MQ=9)           cgcgAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCCCATTCAATAAGCGAATGAGCCtttaaa                                                                                                                                       <  5:104567/148‑4 (MQ=17)           cgcgAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCCCATTCAATAAGCGAATGAGCCttta                                                                                                                                         >  6:104567/1‑145 (MQ=17)            gcgAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGtt                                                                                                                                      <  8:80530/148‑1 (MQ=255)               aaGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGtttttt                                                                                                                                  <  6:346254/149‑1 (MQ=255)                       gTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTtctc                                                                                                                          <  4:79854/149‑1 (MQ=255)                       gTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTtct                                                                                                                           >  3:288662/1‑148 (MQ=255)                          cccTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCaa                                                                                                                        >  5:300740/1‑148 (MQ=255)                          cccTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCaa                                                                                                                        <  4:357421/148‑1 (MQ=255)                                     gTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTcccc                                                                                                             >  5:43052/1‑148 (MQ=255)                                     gTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTcccc                                                                                                             <  6:43052/148‑1 (MQ=255)                                             ggACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAAc                                                                                                     <  2:117685/148‑1 (MQ=255)                                                 acCGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCg                                                                                                >  8:200179/1‑149 (MQ=255)                                                  ccGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAA‑TTTCCCCCGACAAACGGCCGCa                                                                                              <  4:288662/149‑1 (MQ=255)                                                   cGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTcc                                                                                                                                >  2:131107/1‑115 (MQ=255)                                                   cGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACATTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCATTCTCAATAATTTTCCCCCGACAAACGGCCg                                                                                                >  2:374454/1‑147 (MQ=255)                                                    gCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCg                                                                                                >  5:295937/1‑146 (MQ=255)                                                     cccTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCACCTAGGGGCC‑GCCACGATAACCAGCGTATTCAGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTATTGTGTTTTCCAATATCAATAAATTTCCCCCGACAAAAGGCCg                                                                                                >  7:263585/1‑145 (MQ=25)                                                        tttCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGAttt                                                                                         >  6:148102/1‑149 (MQ=255)                                                            aCGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTttgct                                                                                     >  6:102443/1‑149 (MQ=255)                                                             cggcggTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTcc                                                                                                               >  8:84407/1‑122 (MQ=255)                                                             cggcggTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTcc                                                                                                               <  7:84407/122‑1 (MQ=255)                                                               gcggTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTttgcttgc                                                                                  <  2:50229/149‑1 (MQ=255)                                                                   tAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCaaa                                                                              <  2:281908/149‑1 (MQ=255)                                                                    aaCAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCaaa                                                                              >  6:92116/1‑148 (MQ=255)                                                                       aGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCaaaaaa                                                                           >  2:197185/1‑148 (MQ=255)                                                                       aGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCaaaaaa                                                                           >  3:289115/1‑148 (MQ=255)                                                                         gggTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAatat                                                                        >  6:378640/1‑149 (MQ=255)                                                                            ttCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAATatgat                                                                     >  5:71819/1‑149 (MQ=255)                                                                               gAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAATAtgatga                                                                   <  5:194352/148‑1 (MQ=255)                                                                                     tCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAatat                                                                        <  6:35783/140‑1 (MQ=255)                                                                                     tCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAatat                                                                        >  5:35783/1‑140 (MQ=255)                                                                                      ccccTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGGTTGCCAAAAAATATGATGAGaat                                                               >  8:132233/1‑148 (MQ=255)                                                                                        ccTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCTAAAAAATATGATGAGaataa                                                             >  1:15349/1‑148 (MQ=255)                                                                                        ccTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCTAAAAAATATGATGAGaataa                                                             <  2:15349/148‑1 (MQ=255)                                                                                            gggtAC‑GACACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAATATGATGAGAATAATAttt                                                        <  5:215936/149‑1 (MQ=255)                                                                                            ggggAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAATATGATGAGAATAATAtt                                                         >  7:172593/1‑148 (MQ=255)                                                                                               gAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGTCAAACGGCCGCAGATTTTGCTTGCCAAAAAATATGATGAGAATAATATTTCGa                                                     >  5:254930/1‑149 (MQ=255)                                                                                                    ccACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAATATGATGAGAATAATATTTCGatcat                                                 >  4:256529/1‑149 (MQ=255)                                                                                                          ttACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCg                                                                                                <  2:97646/96‑1 (MQ=255)                                                                                                          ttACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCg                                                                                                >  1:97646/1‑96 (MQ=255)                                                                                                               cgcgTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAATATGATGAGAATAATATTTCGATCATCATATTTTGCt                                      <  6:295937/149‑1 (MQ=255)                                                                                                               cgcgTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAATATGATGAGAATAATATTTCGATCATCATATTTTGCt                                      <  5:148102/149‑1 (MQ=255)                                                                                                                gcgTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAATATGATGAGAATAATATTTCGATCATCATATTTTGCtt                                     <  5:92116/149‑1 (MQ=255)                                                                                                                 cgTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAATATGATGAGAATAATATTTCGATCATCATATTTTGCTTg                                    <  6:300740/149‑1 (MQ=255)                                                                                                                 cgTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATATTGCTTGCCAAAAAATATGATGATCGAAATATTATTCTCATCATAttttttggc                                    <  6:196862/149‑6 (MQ=37)                                                                                                                     ttGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAATATGATGAGAATAATATTTCGATCATCATATTTTGCTTGgtg                                 <  5:378640/148‑1 (MQ=255)                                                                                                                        cGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAATATGATGAGAATAATATTTCGATCATCATATTTTGCTTGgtgcgtg                             >  8:229731/1‑149 (MQ=255)                                                                                                                        cGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAACTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAATATGATGAGAATAATATTTCGATCATCATATTTTGCTTGgtgcgt                              <  7:200179/148‑1 (MQ=255)                                                                                                                           gACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAATATGATGAGAATAATATTTCGATCATCATATTTTGCTTGGTGCGTGCCa                          <  3:256529/149‑1 (MQ=255)                                                                                                                               ttAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAACATGATGAGAATAATATTTCGATCATCATAttt                                          >  6:90765/1‑129 (MQ=255)                                                                                                                               ttAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAACATGATGAGAATAATATTTCGATCATCATAttt                                          <  5:90765/129‑1 (MQ=255)                                                                                                                                aaaaGGCCCATTCAATAAGCGAATGAGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAATATGATGAGAATAATATTTCGATCATCATATTTTGCTTGGTGCGTGCCa                          >  2:323925/2‑144 (MQ=255)                                                                                                                                aaaaGGCCCATTCAATAAGCGAATGAGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAATATGATGAGAATAATATTTCGATCATCATATTTTGCTTGGTGCGTGCCa                          <  1:323925/143‑1 (MQ=255)                                                                                                                                aaaaGGCCCATTCAATAAGCGAATGAGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAATATGATGAGAATAATATTTCGATCATCATATTTTGCTTGGTGCGTGCCAGAt                       >  3:96007/2‑147 (MQ=255)                                                                                                                                 aaaGGCCCATTCAATAAGCGAATGAGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAATATGATGAGAATAATATTTCGATCATCATATTTTGCTTGGTGCGTGCCAGATGAg                    <  4:96007/149‑1 (MQ=255)                                                                                                                                  aaGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAATATGATGATAATAATATTTCGATCATCATATTTTGCTTGGTGCGTGCCAGATGAGc                   <  4:289115/149‑1 (MQ=255)                                                                                                                                     gCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAATATGATGAGAATAATATTTCGATCATCATATTTTGCTTGGTGCGTGCCAGATGAGCGa                 >  4:135434/1‑148 (MQ=255)                                                                                                                                              cTTATTGAATGGGCCTTTTGTTTTTTCACTTCTCAATAATTTTCCCCAGACAAACGGCCGCAGATTTTGCTTTCCAAAAAATATGATGATAATAATTTTTCTATCATCATATTATGATTGGTTCGTGCAAGATGAGATACAAGAagag        >  5:250905/1‑148 (MQ=255)                                                                                                                                                     aaTGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAATATGATGAGAATATTATTTCGATCATCATATTTTGCTTGGTGCGTGCCAGATGAGCGACAAGAAGAGCAGAAc  >  2:169364/1‑147 (MQ=255)                                                                                                                                                         |                                                                                                                                               AGTGGATCTTCGCGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAATATGATGAGAATAATATTTCGATCATCATATTTTGCTTGGTGCGTGCCAGATGAGCGACAAGAAGAGCAGAAC  >  NC_002947/2242160‑2242450

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