Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	344,041	Δ1 bp	15.4%	intergenic (+130/+59)	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,041	0	A	.	15.4%	181.6 / 29.6	65	intergenic (+130/+59)	PP_5440/PP_t72	hypothetical protein/tRNA‑Phe
Reads supporting (aligned to +/- strand):  ref base A (27/28);  new base . (6/4);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

GGACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGG  >  NC_002947/343902‑344182                                                                                                                                            |                                                                                                                                               ggACGAGAAAGGTTTGGCCCGCAGACTCCAAATGGCTGCGGGTTTCACGAATGCATGGCGACCCTTCAGCTTGCGTTGTGACGTCTCGTGTTCTTCGACGCCCCCTGAAGTTAAAAATCACAGACAACAAAAAACACGCACTAc‑gcgcg                                                                                                                                      >  4:289213/1‑147 (MQ=255)   aCGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCgg                                                                                                                                       >  5:187598/1‑146 (MQ=255)   aCGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCgg                                                                                                                                       <  6:187598/146‑1 (MQ=255)       gTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCtt                                                                                                                                >  4:32960/1‑149 (MQ=255)        tAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGC‑CTAGTGCGGGTTTTttg                                                                                                                               >  6:349441/1‑149 (MQ=255)        tAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCttg                                                                                                                               >  1:276749/1‑149 (MQ=255)        tAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCttg                                                                                                                               >  4:276571/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                                                                                                                           >  2:227598/1‑149 (MQ=255)            gTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTtgttgt                                                                                                                           <  1:227598/149‑1 (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                                                                                                                                >  8:68275/1‑127 (MQ=39)                             cTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGC‑CTAGTGCGGGtttttt                                                                                                                                <  7:68275/127‑1 (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                                                                      >  5:17471/1‑149 (MQ=255)                                                                 ttCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACg                                                                      >  6:38732/1‑149 (MQ=255)                                                                    aGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGa                                                                    <  4:177912/148‑1 (MQ=255)                                                                      gTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCttg                                                                                                                               >  8:59320/1‑87 (MQ=255)                                                                      gTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCttg                                                                                                                               <  7:59320/87‑1 (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                                                                                                   <  1:249140/92‑1 (MQ=255)                                                                                             cttcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAaccac                                                                                                   >  2:249140/1‑92 (MQ=255)                                                                                                    gccgccTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCctc                                   <  2:292809/149‑1 (MQ=255)                                                                                                     ccgccTGTAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCctct                                  <  3:289213/149‑1 (MQ=255)                                                                                                       gccTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAAcc                                                            <  1:35206/121‑1 (MQ=255)                                                                                                       gccTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAAcc                                                            >  2:35206/1‑121 (MQ=255)                                                                                                       gccTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTc                                          >  4:119008/1‑139 (MQ=255)                                                                                                       gccTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTc                                          <  3:119008/139‑1 (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                    <  7:49463/149‑1 (MQ=255)                                                                                                                   aaaTCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGa                    <  6:17471/149‑1 (MQ=255)                                                                                                                   aaaTCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGa                    <  5:38732/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)                                                                                                                                 aagaaaCCCGC‑CTAGTGCGGGTTTTTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGcgac      >  2:325355/4‑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)                                                                                                                                     aaCCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAgg  >  3:163733/1‑149 (MQ=255)                                                                                                                                            |                                                                                                                                               GGACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAG‑GCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGG  >  NC_002947/343902‑344182

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