Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,914,874	C→T	7.8%	intergenic (‑88/‑158)	PP_3454 ← / → PP_3455	transcriptional regulator/multidrug RND transporter membrane fusion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,914,874	0	C	T	7.8%	187.3 / 3.5	66	intergenic (‑88/‑158)	PP_3454/PP_3455	transcriptional regulator/multidrug RND transporter membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base C (27/33);  new base T (2/3);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.62e-01

GTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTAGCATATACACGGTCATCATC  >  NC_002947/3914727‑3915020                                                                                                                                                    |                                                                                                                                                   gTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCa                                                                                                                                                   >  6:299337/1‑149 (MQ=255)   aGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCAcc                                                                                                                                                 >  6:327624/1‑149 (MQ=255)     cTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTGAACCCGCTCCTACGGGGACCTGCAGCAGCGGGTTTACCcg                                                                                                                                               <  5:3258/149‑1 (MQ=25)      tGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCcgc                                                                                                                                              >  8:329547/1‑149 (MQ=255)      tGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCcgc                                                                                                                                              >  5:120691/1‑149 (MQ=255)      tGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCcgc                                                                                                                                              >  1:378441/1‑149 (MQ=255)      nggCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCACCGTAGGAGCGGGTTCACCcgc                                                                                                                                              >  2:109996/2‑149 (MQ=255)        gCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGGTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTGAACCcgctc                                                                                                                                            <  3:295854/149‑3 (MQ=255)           aTCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                                         <  6:54412/149‑1 (MQ=255)            tCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAagag                                                                                                                                        <  3:71068/149‑1 (MQ=255)                 tCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAg                                                                                                                                   >  8:4/1‑149 (MQ=255)                      cAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAg                                                                                                                              <  5:299337/149‑1 (MQ=255)                        gggCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGc                                                                                                                            <  5:327624/149‑1 (MQ=255)                              tGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGGTCACCCGCGAAGCGGCCAGCACAGGCCCAGGc                                                                                                                      >  4:322143/1‑149 (MQ=255)                                  gtcTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGaa                                                                                                                  >  4:351573/1‑149 (MQ=255)                                  gtcTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGa                                                                                                                   <  6:92653/148‑1 (MQ=255)                                   tcTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCGGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAat                                                                                                                 <  1:109996/149‑1 (MQ=255)                                         aCCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGaa                                                                                                           <  6:120691/149‑1 (MQ=255)                                           cAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAAcc                                                                                                         <  8:158411/149‑1 (MQ=255)                                            aGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGc                                                                                                                                      >  6:139837/1‑119 (MQ=255)                                              gAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCaa                                                                                                      >  3:14085/1‑149 (MQ=255)                                              gAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCaa                                                                                                      >  6:279372/1‑149 (MQ=255)                                               aaaaTGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAAc                                                                                                     <  7:257642/149‑1 (MQ=255)                                                   tGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACt                                                                                                 <  7:329547/149‑1 (MQ=255)                                                               gTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCAt                                                                                     >  4:82415/1‑149 (MQ=255)                                                                   ttGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCAttctt                                                                                 <  5:34776/149‑1 (MQ=255)                                                                     gggCGTCTGGGTTTAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCAttctt                                                                                 <  7:165714/147‑1 (MQ=255)                                                                              ggTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACttttt                                                                       >  6:38989/1‑148 (MQ=255)                                                                              ggTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACttttt                                                                       <  7:246178/148‑1 (MQ=255)                                                                               gTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTca                                                                     >  2:174805/1‑149 (MQ=255)                                                                                     gTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCACCCACTCGACAGGATCATTCTTCACACTTTTTCACACTCg                                                               >  1:57757/1‑149 (MQ=255)                                                                                      tAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGc                                                              >  2:296188/1‑149 (MQ=255)                                                                                           cTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGTAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAAGCACTCGCCAGGATCATTGTTCACACTTTTTCACACTCTCCGTa                                                          >  5:186927/1‑148 (MQ=38)                                                                                           cTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTAc                                                         <  4:247670/149‑1 (MQ=255)                                                                                                 cGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCCCACTTTTTCACACTCGCCCTACAGGTAt                                                   <  1:214481/149‑1 (MQ=255)                                                                                                  gAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTAtt                                                  <  6:25602/149‑1 (MQ=255)                                                                                                    gcgcTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTaa                                                <  2:201104/149‑1 (MQ=255)                                                                                                       cTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGCAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACCCTTTTTCACACTCGCCCTACAGGTATTAACAg                                             >  3:2950/1‑149 (MQ=255)                                                                                                          cgcgGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAgc                                            <  5:38989/147‑1 (MQ=255)                                                                                                             gggTAAACCCGTTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATACAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAg                                                       <  3:105468/133‑1 (MQ=255)                                                                                                              ggTGAACCCGCTCCTACGGGGACCTGCAGCAGCGGGTTTACCCGCGAAGCGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTAtt                                                  <  4:226516/137‑1 (MQ=25)                                                                                                              ggTGAACCCGCTCCTACGGGGACCTGCAGCAGCGGGTTTACCCGCGAAGCGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTAtt                                                  >  3:226516/1‑137 (MQ=25)                                                                                                              ggTGAACCCGCTCCTACGGGGACCTGCAGCAGCGGGTTTACCCGCGAAGCGCTCGCCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACa                                              <  5:221871/141‑1 (MQ=18)                                                                                                              ggTGAACCCGCTCCTACGGGGACCTGCAGCAGCGGGTTTACCCGCGAAGCGCTCGCCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACa                                              >  6:221871/1‑141 (MQ=18)                                                                                                                tAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATACTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAgg                                    >  2:77688/1‑149 (MQ=255)                                                                                                                   aCCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACCCTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAg                                 >  5:151266/1‑149 (MQ=255)                                                                                                                    cccGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGc                                <  3:351573/149‑1 (MQ=255)                                                                                                                       gctgctGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTc                             <  1:296188/149‑1 (MQ=255)                                                                                                                        ctgctgCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCt                            <  2:378441/149‑1 (MQ=255)                                                                                                                          gctgcAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGcc                                                             <  6:42516/114‑1 (MQ=255)                                                                                                                          gctgcAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGcc                                                             >  5:42516/1‑114 (MQ=255)                                                                                                                             gcAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCt                        >  6:136237/1‑148 (MQ=255)                                                                                                                                 gTCCCCGTAGCAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATATAACACAACCACTCGACAGCATCATTCTTCACACTTTTTAACACACGCCCTACAGGTATTAACAGCCGCCAAGGCAGCGTCTGCCATAGCa                   >  4:311111/1‑149 (MQ=38)                                                                                                                                   ccccGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTAGCata                 <  4:170273/149‑1 (MQ=255)                                                                                                                                     ccGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCGCTACAGGTATTAACAGCCGCCCAGGCAGCGTc                             <  5:110390/135‑1 (MQ=255)                                                                                                                                     ccGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTc                             >  6:110390/1‑135 (MQ=255)                                                                                                                                      cGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTAGCATATac              <  7:293828/149‑1 (MQ=255)                                                                                                                                        tAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCACCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTAGCATATacac            <  3:82415/149‑1 (MQ=255)                                                                                                                                          ggAGCGGTTTCACGCGCGAAGAGGCCTGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACTCTTTTTCACACTCGACCTACAGGTTTTAACATCCGCCCAGGCAGCGTCTGCGTTAGCATATACACgg          <  3:311111/149‑1 (MQ=255)                                                                                                                                          ggAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTAGCATATACACg           <  1:174805/148‑1 (MQ=255)                                                                                                                                           gAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTAGCATATACACGGt         >  1:289241/1‑149 (MQ=255)                                                                                                                                            aGCTGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTAGCATATACACGGtc        <  2:57757/149‑1 (MQ=255)                                                                                                                                            aGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTAGCATATACACGGtc        <  7:4/149‑1 (MQ=255)                                                                                                                                            aGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACAATCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTAGCATATACACGGtc        <  5:279372/149‑1 (MQ=255)                                                                                                                                             gCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTAGCATATACACGGtca       >  6:360863/1‑149 (MQ=255)                                                                                                                                                  ttCACCCGCGAAGAGGCCAGCACAGTCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTAGCATATACACGGTCatcatc  <  6:186927/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTAGCATATACACGGTCATCATC  >  NC_002947/3914727‑3915020

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