Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,735,139	Δ1 bp	13.6%	coding (36/3093 nt)	PP_3302 →	RND family transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,735,139	0	A	.	13.6%	185.8 / 26.7	66	coding (36/3093 nt)	PP_3302	RND family transporter
Reads supporting (aligned to +/- strand):  ref base A (29/28);  new base . (2/7);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.45e-01

GCACGCTGGCGATCGGATCGTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTCAAGCAGATGACGATTATCACTGCCTGGCCCGGCGCCACTGCG  >  NC_002947/3735002‑3735286                                                                                                                                          |                                                                                                                                                      gCACGCTGGCGATCGGATCGTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑gcgagc                                                                                                                                            >  8:312682/1‑148 (MQ=255)  cACGCTGGCGATCGGATCGTGGCGCTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGC‑GTACTgcgagcgc                                                                                                                                          <  5:353837/149‑2 (MQ=255)  cACGCTGGCGATCGGATCGTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGC‑GTACTgcgagcg                                                                                                                                           <  3:22231/148‑1 (MQ=255)  cACGCTGGCGATCGGATCGTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCGCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGC‑GTACTgcgagcgc                                                                                                                                          <  1:146383/149‑2 (MQ=255)   aCGGTGGCGATCGGATCGTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGC‑GTACTgcgagcgcg                                                                                                                                         <  3:17609/149‑3 (MQ=255)   aCGCTGGCGATCGGATCGTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGC‑GTACTgcgagcgcg                                                                                                                                         <  1:240483/149‑3 (MQ=255)   aCGCTGGCGATCGGATCGTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGC‑GTACTgcgagcgcg                                                                                                                                         <  3:44623/149‑3 (MQ=255)     gCTGGCGATCGGATCGTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGC‑GTACTgcgagcgcgga                                                                                                                                       <  2:27501/149‑5 (MQ=255)        ggCGATCGGATCGTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGatcaat                                                                                                                                     >  2:347214/1‑148 (MQ=255)        ggCGATCGGATCGTCGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGatcaatc                                                                                                                                    <  7:173261/149‑1 (MQ=255)         gCGATCGGATCGTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGatcaatca                                                                                                                                   >  5:247583/1‑149 (MQ=255)         gCGATCGGATCGTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGatcaatc                                                                                                                                    <  7:234075/148‑1 (MQ=255)                gatcgTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCtt                                                                                                                            >  2:42587/1‑149 (MQ=255)                 atcgTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTc                                                                                                                           <  3:340249/149‑1 (MQ=255)                     tGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTCC‑GCGAGCGATCAATCACCCTCTTCCtca                                                                                                                       >  4:26386/1‑147 (MQ=255)                     tGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTcc                                                                                                                          <  3:337319/146‑1 (MQ=255)                         aCTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATtgt                                                                                                                   >  7:313992/1‑149 (MQ=255)                                  gcATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGa                                                                                                                       >  7:101205/1‑136 (MQ=255)                                  gcATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGa                                                                                                                       <  8:101205/136‑1 (MQ=255)                                  gcATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑c                                                                                                         >  8:178901/1‑149 (MQ=255)                                         ctgcGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGcc                                                                                                  >  6:36673/1‑149 (MQ=255)                                               gATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACg                                                                                            >  5:125758/1‑149 (MQ=255)                                                      tCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGTCCt                                                                                     >  1:367267/1‑149 (MQ=255)                                                        gAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCttt                                                                                   <  4:45213/149‑1 (MQ=255)                                                         aGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCtttt                                                                                  >  7:23616/1‑149 (MQ=255)                                                           gTCATGGAACAGGCTCAAACAGCTCACCTCATGGGGCACCAACAATGAGCCAGCTGCGTTTCAACTTATCAGCGCTCGCAGTAC‑GCGAGCTATCAATCACCCTCTACCAGATAGTGCTGATCC‑CCTTCGCCGACACGCTGCCCttttt                                                                                 >  8:265809/1‑148 (MQ=255)                                                              aTGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCaa                                                                              <  7:200225/148‑1 (MQ=255)                                                              aTGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAg                                                                             >  2:107824/1‑149 (MQ=255)                                                                 gAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCt                                                                           <  6:247583/148‑1 (MQ=255)                                                                 gAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTg                                                                          <  2:299062/149‑1 (MQ=255)                                                                 gAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTg                                                                          >  3:256912/1‑149 (MQ=255)                                                                  aaCAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTgg                                                                         <  7:209936/149‑1 (MQ=255)                                                                     aGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGcc                                                                      <  7:280898/149‑1 (MQ=255)                                                                       gCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGt                                                                    >  8:170579/1‑149 (MQ=255)                                                                       gCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGt                                                                    >  6:154289/1‑149 (MQ=255)                                                                         tCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTg                                                                   >  7:286421/1‑148 (MQ=255)                                                                          cAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTg                                                                   >  6:19452/1‑147 (MQ=255)                                                                            aaCAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTCC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCag                                                                >  6:164810/1‑148 (MQ=255)                                                                             aCAGCTCAGGTCATTGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑gcgagc                                                                                                                                            <  2:365271/72‑1 (MQ=255)                                                                             aCAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑gcgagc                                                                                                                                            >  1:365271/1‑72 (MQ=255)                                                                             aCAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCt                                                                                                       >  5:281309/1‑108 (MQ=255)                                                                             aCAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCt                                                                                                       <  6:281309/108‑1 (MQ=255)                                                                             aCAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAACGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCagaa                                                              <  1:158362/149‑1 (MQ=255)                                                                              cAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCagaag                                                             <  7:178901/149‑1 (MQ=255)                                                                                    aGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATcc                                                         <  7:363191/147‑1 (MQ=255)                                                                                      gTCATGGGGCCCCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGccc                                                     >  6:49703/1‑149 (MQ=255)                                                                                      gTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGccc                                                     <  8:23616/149‑1 (MQ=255)                                                                                        cATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCtt                                                   >  3:348818/1‑149 (MQ=255)                                                                                         aTGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCtt                                                   <  5:154289/148‑1 (MQ=255)                                                                                          tGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCGCCCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTc                                                  <  1:107824/149‑1 (MQ=255)                                                                                          tGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCa                                                 <  5:19452/149‑1 (MQ=255)                                                                                            gggCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCAcc                                               <  8:313992/149‑1 (MQ=255)                                                                                              gCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGt                                             <  1:42587/149‑1 (MQ=255)                                                                                               cACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTc                                            <  6:125758/149‑1 (MQ=255)                                                                                                  caaccaTGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTCAAg                                         <  7:312682/149‑1 (MQ=255)                                                                                                      caTGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTCAAGCAGa                                     <  5:36673/149‑1 (MQ=255)                                                                                                      caTGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTCAAGCAGa                                     <  4:348818/149‑1 (MQ=255)                                                                                                      caTGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTCAAGCAGa                                     <  8:286421/149‑1 (MQ=255)                                                                                                                         aaCTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTCAAGCAGATGACGATAATCACTGCCTg                  >  5:65161/1‑149 (MQ=255)                                                                                                                          aCTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTCAAGCAGATGACGATTATCACTGCCTgg                 >  6:131397/1‑149 (MQ=255)                                                                                                                            ttATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTCAAGCAGATGACGATTATCACTGCCTGGcc               >  2:198841/1‑149 (MQ=255)                                                                                                                                 atcgcTCGC‑GTACTGCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTCAAGCAGATGACGATTATCACTGCCTGGCCCGgcg          >  6:110941/3‑149 (MQ=255)                                                                                                                                  tcgcTCGC‑GTACTGCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTCAAGCAGATGACGATTATCACTGCCTGGCCCGgcgc         >  6:302645/2‑149 (MQ=255)                                                                                                                                        gCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTCAAGCAGATGACGATTATCACTGCCTGGCCCGGCGCCACTgc   >  3:7470/1‑149 (MQ=255)                                                                                                                                        gCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTCAAGCAGATGACGATTATCACTGCCTGGCCCGGCGCCACTgc   >  1:169060/1‑149 (MQ=255)                                                                                                                                         cAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTCAAGCAGATGACGATTATCACTGCCTGGCCCGGCGCCACTgcg  <  5:49703/149‑1 (MQ=255)                                                                                                                                          |                                                                                                                                                      GCACGCTGGCGATCGGATCGTGGCACTGGGTGCGCATCTGCTGCGCGATGGCGTCGAGGTCATGGAACAGGCTCAAACAGCTCAGGTCATGGGGCACCAACCATGAGCCAGCTGCGTTTCAACTTATCCGCGCTCGCAGTAC‑GCGAGCGATCAATCACCCTCTTCCTGATTGTGCTGATCG‑CCTTCGCCGGCACGCTGGCCTTTTTCAAGCTGGGCCGTGCAGAAGATCCGCCCTTCACCGTCAAGCAGATGACGATTATCACTGCCTGGCCCGGCGCCACTGCG  >  NC_002947/3735002‑3735286

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