Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,242,290	T→C	9.1%	intergenic (+33/‑138)	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,290	0	T	C	9.1%	182.4 / 5.2	65	intergenic (+33/‑138)	PP_t30/PP_1978	tRNA‑Glu/TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (29/30);  new base C (3/3);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGGGTTCACGCAAGTGGATCTTCGCGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAATATGATGAGAATAATATTTCGATCATCATATTTTGCTTGGTGCGTGCCAGATGAGCGA  >  NC_002947/2242146‑2242435                                                                                                                                                     |                                                                                                                                                  ccAGGGTTCACGCAAGTGGATCTTCGCGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCAt                                                                                                                                                <  1:252325/148‑1 (MQ=255)  cAGGGTTCACGCAAGTGGATCTTCGCGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTc                                                                                                                                              >  7:80530/1‑149 (MQ=255)          aCGCAAGGGGATCTTCGCGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTAtt                                                                                                                                       <  5:254386/148‑1 (MQ=255)             cAAGTGGATCTTCGCGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAAt                                                                                                                                   <  3:173114/149‑1 (MQ=255)               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‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAAt                                                                                                                                   <  1:187648/125‑1 (MQ=255)                                     gTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAAt                                                                                                                                   >  2:187648/1‑125 (MQ=255)                                     gTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTtctc                                                                                                           <  4:79854/149‑1 (MQ=255)                                     gTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTtct                                                                                                            >  3:288662/1‑148 (MQ=255)                                        cccTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCaa                                                                                                         <  4:357421/148‑1 (MQ=255)                                        cccTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCaa                                                                                                         >  5:300740/1‑148 (MQ=255)                                                   gTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTcccc                                                                                              <  6:43052/148‑1 (MQ=255)                                                   gTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTcccc                                                                                              >  5:43052/1‑148 (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                                                         >  5:35783/1‑140 (MQ=255)                                                                                                   tCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAatat                                                         <  6:35783/140‑1 (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                                                                                 >  1:97646/1‑96 (MQ=255)                                                                                                                        ttACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCg                                                                                 <  2:97646/96‑1 (MQ=255)                                                                                                                             cgcgTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAATATGATGAGAATAATATTTCGATCATCATATTTTGCt                       <  5:148102/149‑1 (MQ=255)                                                                                                                             cgcgTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAATATGATGAGAATAATATTTCGATCATCATATTTTGCt                       <  6:295937/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)                                                                                                                                                     |                                                                                                                                                  CCAGGGTTCACGCAAGTGGATCTTCGCGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAATATGATGAGAATAATATTTCGATCATCATATTTTGCTTGGTGCGTGCCAGATGAGCGA  >  NC_002947/2242146‑2242435

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