Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,242,296	G→A	9.5%	intergenic (+39/‑132)	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,296	0	G	A	9.5%	175.6 / 7.3	63	intergenic (+39/‑132)	PP_t30/PP_1978	tRNA‑Glu/TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (28/29);  new base A (3/3);  total (31/32)
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.68e-01

ACGCAAGTGGATCTTCGCGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAATATGATGAGAATAATATTTCGATCATCATATTTTGCTTGGTGCGTGCCAGATGAGCGA  >  NC_002947/2242155‑2242435                                                                                                                                                  |                                                                                                                                            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                                              <  2:15349/148‑1 (MQ=255)                                                                                             ccTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCTAAAAAATATGATGAGaataa                                              >  1:15349/1‑148 (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                       <  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)                                                                                                                                                  |                                                                                                                                            ACGCAAGTGGATCTTCGCGAAGGTTACGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAG‑‑‑TCCCCTAGGGGAC‑GCCACGTTTACCCGCGTTTTGCGGGACTTTAAAGGCTCATTCGCTTATTGAATGGGCCTTTTGTTTTTTCCCTTCTCAATAATTTTCCCCCGACAAACGGCCGCAGATTTTGCTTGCCAAAAAATATGATGAGAATAATATTTCGATCATCATATTTTGCTTGGTGCGTGCCAGATGAGCGA  >  NC_002947/2242155‑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: ‑