Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	344,032	A→G	15.1%	intergenic (+121/+68)	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,032	0	A	G	15.1%	132.4 / 16.7	60	intergenic (+121/+68)	PP_5440/PP_t72	hypothetical protein/tRNA‑Phe
Reads supporting (aligned to +/- strand):  ref base A (24/26);  new base G (5/4);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.47e-01

GGACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGAC  >  NC_002947/343902‑344178                                                                                                                                   |                                                                                                                                                    ggACGAGAAAGGTTTGGCCCGCAGACTCCAAATGGCTGCGGGTTTCACGAATGCATGGCGACCCTTCAGCTTGCGTTGTGACGTCTCGTGTTCTTCGACGCCCCCTGAAGTTAAAAATCACAGACAACAAAAAACACGCACTAc‑gcgcg                                                                                                                                  >  4:289213/1‑147 (MQ=255)   aCGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCgg                                                                                                                                   <  6:187598/146‑1 (MQ=255)   aCGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCgg                                                                                                                                   >  5:187598/1‑146 (MQ=255)       gTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCtt                                                                                                                            >  4:32960/1‑149 (MQ=255)        tAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGC‑CTAGTGCGGGTTTTttg                                                                                                                           >  6:349441/1‑149 (MQ=255)        tAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCttg                                                                                                                           >  4:276571/1‑149 (MQ=255)        tAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCttg                                                                                                                           >  1:276749/1‑149 (MQ=255)          aGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGGCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCttgt                                                                                                                          <  2:310590/148‑1 (MQ=255)           ggTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGC‑CTAGTGCGGGTTTTttgttg                                                                                                                        <  7:235940/149‑1 (MQ=39)           ggTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGC‑CTAGTGCGGGTTTTttgtt                                                                                                                         <  5:349441/148‑1 (MQ=39)            gTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTtgttgt                                                                                                                       <  1:227598/149‑1 (MQ=255)            gTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTtgttgt                                                                                                                       >  2:227598/1‑149 (MQ=255)               tGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCg                                                                                                                     <  3:126280/148‑1 (MQ=255)                  cccGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTc                                                                                                                 <  2:276749/149‑1 (MQ=255)                   ccGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGC‑CTAGTGCGGGtttttt                                                                                                                            <  2:225857/137‑1 (MQ=39)                   ccGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGC‑CTAGTGCGGGtttttt                                                                                                                            >  1:225857/1‑137 (MQ=39)                    cGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACACAAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCgg                                                                                                               >  3:211958/1‑149 (MQ=255)                           tGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTGTCTTGTTGTCGCTTCGGGTAACtt                                                                                                        >  3:73584/1‑149 (MQ=255)                            gCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACttt                                                                                                       >  3:177912/1‑149 (MQ=255)                             cTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGC‑CTAGTGCGGGtttttt                                                                                                                            <  7:68275/127‑1 (MQ=39)                             cTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGC‑CTAGTGCGGGtttttt                                                                                                                            >  8:68275/1‑127 (MQ=39)                              tCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCgg                                                                                                                                   <  2:261498/119‑1 (MQ=255)                              tCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCgg                                                                                                                                   >  1:261498/1‑119 (MQ=255)                                                 gATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGaa                                                                                   <  3:276571/148‑1 (MQ=255)                                                  aTTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAAccacca                                                                                             >  5:186690/1‑137 (MQ=255)                                                  aTTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAAccacca                                                                                             <  6:186690/137‑1 (MQ=255)                                                     gCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAggtgg                                                                              >  6:324674/1‑149 (MQ=255)                                                          gCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGcc                                                                          <  3:32960/148‑1 (MQ=255)                                                              gCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCaga                                                                      <  5:324674/148‑1 (MQ=255)                                                                cTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGAc                                                                   <  4:73584/149‑1 (MQ=255)                                                                 ttCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACg                                                                  >  6:38732/1‑149 (MQ=255)                                                                 ttCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACg                                                                  >  5:17471/1‑149 (MQ=255)                                                                    aGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGa                                                                <  4:177912/148‑1 (MQ=255)                                                                      gTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCttg                                                                                                                           <  7:59320/87‑1 (MQ=255)                                                                      gTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCttg                                                                                                                           >  8:59320/1‑87 (MQ=255)                                                                          cGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAAc                                                         <  8:209504/149‑1 (MQ=255)                                                                            ttGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAAccg                                                       >  7:200687/1‑149 (MQ=255)                                                                                    tctcGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAg                                                                                  >  8:326567/1‑114 (MQ=255)                                                                                    tctcGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAg                                                                                  <  7:326567/114‑1 (MQ=255)                                                                                      tcGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAggt                                                                                <  6:273768/114‑1 (MQ=255)                                                                                      tcGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAAc                                                         <  1:142132/137‑1 (MQ=255)                                                                                        gggctcttcGAAGCCGCCGTAAGTTAAAACTCCCAGACCACAACAACCCCGCACTAG‑GCGGGGTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCCGAAGCTTGAAGGTGGTGCCCAGCGACGGAGTCGAACCGCCGACACGa                                              <  4:211958/142‑1 (MQ=255)                                                                                             cttcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAaccac                                                                                               >  2:249140/1‑92 (MQ=255)                                                                                             cttcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAaccac                                                                                               <  1:249140/92‑1 (MQ=255)                                                                                                    gccgccTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCctc                               <  2:292809/149‑1 (MQ=255)                                                                                                     ccgccTGTAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCctct                              <  3:289213/149‑1 (MQ=255)                                                                                                       gccTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAAcc                                                        >  2:35206/1‑121 (MQ=255)                                                                                                       gccTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAAcc                                                        <  1:35206/121‑1 (MQ=255)                                                                                                       gccTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTc                                      <  3:119008/139‑1 (MQ=255)                                                                                                       gccTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTc                                      >  4:119008/1‑139 (MQ=255)                                                                                                               ttAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGAtttt                                       >  4:285392/1‑130 (MQ=255)                                                                                                               ttAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGAtttt                                       <  3:285392/130‑1 (MQ=255)                                                                                                                tAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGAc                   >  5:69518/1‑149 (MQ=255)                                                                                                                   aaaTCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGa                <  6:17471/149‑1 (MQ=255)                                                                                                                   aaaTCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGa                <  5:38732/149‑1 (MQ=255)                                                                                                                   aaaTCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGa                <  7:49463/149‑1 (MQ=255)                                                                                                                       cacaGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTa            >  2:210200/1‑149 (MQ=255)                                                                                                                           gACAACAAGAAACCCGC‑CTAGTGCGGGTTTTTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTg        >  2:358341/1‑149 (MQ=39)                                                                                                                              aacaaAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGc     >  3:130274/1‑149 (MQ=255)                                                                                                                                caaGAAACCCGC‑CTAGCGCGGGTTTTTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGcga   >  8:4919/1‑149 (MQ=255)                                                                                                                                 aaaaaaCCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTgggagct  >  5:241241/1‑145 (MQ=255)                                                                                                                                 aaaaaaCCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGcgac  >  1:169731/1‑149 (MQ=255)                                                                                                                                 aaaaaaCCCGCACTAG‑GCGGGTTTCTTGTTGGCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTgg       <  6:241241/144‑1 (MQ=255)                                                                                                                                   |                                                                                                                                                    GGACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGAC  >  NC_002947/343902‑344178

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