Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,069,557	C→T	7.0%	S17S (TCG→TCA)	PP_5323 ←	M23/M37 family peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,069,557	0	C	T	7.0%	206.8 / 4.0	71	S17S (TCG→TCA)	PP_5323	M23/M37 family peptidase
Reads supporting (aligned to +/- strand):  ref base C (29/37);  new base T (2/3);  total (31/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.88e-01

TCGCGTTCAGCAGTACCTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTATTGCAGTTGGTCGGATGGGGCAACTCGGGCTTGG  >  NC_002947/6069415‑6069702                                                                                                                                               |                                                                                                                                                    tCGCGGTCAGCAGTACCTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCTAGGGCTTGGCGCCGGGGCCAGGGCGATCTGAATAGGAGATGACACCGGAATCGTCGGTTATCTTGTAGAGCGTCACGCCCTGGCTCGACGc                                                                                                                                                <  7:123839/148‑1 (MQ=255)   gcgTTCAGCAGTACCTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCGAGCCAGGGCGTGACGCtct                                                                                                                                             <  4:175855/149‑3 (MQ=255)   gcgTTCAGCAGTACCTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGa                                                                                                                                             <  1:109087/149‑1 (MQ=255)    cgTTCAGCAGTACCTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑g                                                                                                                                           <  1:156591/149‑1 (MQ=255)            agTACCTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAg                                                                                                                                   <  4:297567/149‑1 (MQ=255)             gTACCTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCCTCACGCCCTGGCTCGACGCTGA‑GGCCATGAg                                                                                                                                   >  7:300046/1‑148 (MQ=255)                ccTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAAt                                                                                                                               >  4:36221/1‑149 (MQ=255)                 cTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATg                                                                                                                              <  7:244491/149‑1 (MQ=255)                 cTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATg                                                                                                                              <  3:273157/149‑1 (MQ=255)                 cTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGAACTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATg                                                                                                                              >  6:288345/1‑149 (MQ=255)                   ggCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGTGCCAGGGAGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTAGTAGAGCGTCACGACCTGGCTCGACGCTGA‑GTCAATGAGAAAtta                                                                                                                             >  8:146393/1‑146 (MQ=255)                   ggCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCGAGCCAGGGCGTGACGCtctacaa                                                                                                                                         >  4:338670/1‑131 (MQ=255)                   ggCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCGAGCCAGGGCGTGACGCtctacaa                                                                                                                                         <  3:338670/137‑7 (MQ=255)                    gCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCGTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑cac                                                                                                                          <  6:322120/149‑1 (MQ=255)                      cTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAg                                                                                                                        <  7:155735/149‑1 (MQ=255)                             cGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAAcagca                                                                                                                 <  3:36221/149‑1 (MQ=255)                             cGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAAcagca                                                                                                                 <  4:214513/149‑1 (MQ=255)                                   gACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGgc                                                                                                           >  1:294686/1‑149 (MQ=255)                                   gACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGgc                                                                                                           >  2:219450/1‑149 (MQ=255)                                    aCCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGgcg                                                                                                          >  4:349323/1‑149 (MQ=255)                                     ccATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGCGCAATG‑CACAGAACCGCAGCAGgcgc                                                                                                         >  8:139352/1‑149 (MQ=255)                                        tGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGcgcgcg                                                                                                      <  3:327675/149‑1 (MQ=255)                                           gCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCggg                                                                                                    >  7:340210/1‑148 (MQ=255)                                            cTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAAcagcagcag                                                                                                             <  3:153701/138‑1 (MQ=255)                                            cTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAAcagcagcag                                                                                                             >  4:153701/1‑138 (MQ=255)                                                aGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCGAGCCAGGGCGTGACGCtctac                                                                                                                                           <  8:227146/106‑5 (MQ=38)                                                aGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATATTGTAGAGCGTCGAGCCAGGGCGTGACGCtctac                                                                                                                                           >  7:227146/1‑102 (MQ=37)                                                  cAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATGCCACAGAACAGCAGCAGGCGCGCGGGCATGATa                                                                                             <  1:37612/149‑1 (MQ=255)                                                    acacCAAGGGCTTGGCGCCGGGGCCAGGGCGATCTGAATATGAGATGACACCGCAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCTATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTg                                                                                          <  2:318236/149‑1 (MQ=255)                                                    acacCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTg                                                                                          <  1:219450/149‑1 (MQ=255)                                                        cAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGt                                                                                           >  4:54876/1‑144 (MQ=255)                                                        cAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGt                                                                                           <  3:54876/144‑1 (MQ=255)                                                        cAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCtt                                                                                      >  5:186300/1‑149 (MQ=255)                                                         aaGGGCTTGGCGCTGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCttt                                                                                     <  7:134250/149‑1 (MQ=255)                                                          aGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGAc                                                                                                                                                  <  4:259858/89‑1 (MQ=255)                                                          aGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGAc                                                                                                                                                  >  3:259858/1‑89 (MQ=255)                                                           gggCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCttttt                                                                                   <  3:349323/149‑1 (MQ=255)                                                               ttGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGa                                                                                 >  8:41227/1‑147 (MQ=255)                                                                  gcgcCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTGTTTGACAGt                                                                             >  5:258502/1‑148 (MQ=255)                                                                           ccAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGt                                                                   <  8:93153/149‑1 (MQ=255)                                                                             aGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCg                                                                 <  8:340210/149‑1 (MQ=255)                                                                               ggCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGGCg                                                                 >  8:334722/1‑147 (MQ=255)                                                                                gCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATcaaatcgg                                                               >  3:237554/1‑141 (MQ=255)                                                                                       gAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGaa                                                               >  2:172471/1‑141 (MQ=255)                                                                                       gAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGaa                                                               <  1:172471/141‑1 (MQ=255)                                                                                        aaTAGTAGATGACACCGGAATCTTCGGTGATCTTGTAGAGCGTCTCGCCCTGGCTCGACGCTGA‑TGCCATGAGCAATG‑CACAGAACAGCACCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTa                                                          <  7:146393/145‑1 (MQ=255)                                                                                         aTAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTagcagc                                                     >  1:10094/1‑149 (MQ=255)                                                                                              agaTGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGtt                                                <  5:288345/149‑1 (MQ=255)                                                                                               gaTGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGttt                                               <  6:186300/149‑1 (MQ=255)                                                                                                  gACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑g                                                                                                                                           >  6:64565/1‑55 (MQ=255)                                                                                                  gACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑g                                                                                                                                           <  5:64565/55‑1 (MQ=255)                                                                                                  gACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTa                                             <  6:305492/148‑1 (MQ=255)                                                                                                   acacCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACg                                           <  8:29004/149‑1 (MQ=255)                                                                                                      ccGGAATCGTCGGTGATCTTGTAGAGCGTCACTCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACagta                                        <  6:258502/149‑1 (MQ=255)                                                                                                       cGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGAAGCGGGTTTTACggtag                                       >  4:197271/1‑149 (MQ=255)                                                                                                       aggAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACggtag                                       >  2:323500/2‑149 (MQ=255)                                                                                                                     aTCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTATTGCAGTTgg                          <  2:10094/148‑1 (MQ=255)                                                                                                                         tGTAGAGCGTCACGCCCTGGCTCGACGCTGATGGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTAtt                                  >  1:38032/1‑137 (MQ=255)                                                                                                                         tGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTAtt                                  <  2:38032/136‑1 (MQ=255)                                                                                                                          gTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTATTGCAGTTGGTCGGAt                    >  6:200259/1‑149 (MQ=255)                                                                                                                            agagCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTATTGCAGTTGGTCGGATg                   <  7:242627/148‑1 (MQ=255)                                                                                                                               gCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTATTGCAGTTGGTCGGATgggg                <  5:200259/148‑1 (MQ=255)                                                                                                                                  tCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTATTGCAGTTGGTCGGATGGGGCAAc            >  6:331396/1‑149 (MQ=255)                                                                                                                                      gCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTATTGCAGTTGGTCGGATGGGGCAACTCgg        >  2:12583/1‑149 (MQ=255)                                                                                                                                        ccTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTATTGCAGTTGGTCGGATGGGGCAACTCGGGc      >  4:327150/1‑149 (MQ=255)                                                                                                                                          tGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTATTGCAGTTGGTCGGATGGGGCACCTCGGGCtt    >  5:36058/1‑149 (MQ=255)                                                                                                                                           ggCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTATTGCAGTTGGTCGGATGGGGCAACTCGGGCTTg   >  8:185341/1‑149 (MQ=255)                                                                                                                                            gCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGACCAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTATTGCAGTTGGTCGGATGGGGCAACTCGGGCTTg   <  7:139352/148‑1 (MQ=255)                                                                                                                                            gCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTATTGCAGTTGGTCGGATGGGGCAACTCGGGCTTgg  <  2:294686/149‑1 (MQ=255)                                                                                                                                            gCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTATTGCAGTTGGTCGGATGGGGCAACTCGGGCTTgg  <  7:334722/149‑1 (MQ=255)                                                                                                                                              tCGAAGCTGA‑GGCCATTAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCAGGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCATAGCAGCGGGTTTTACGGTAGGTATTGCAGTTGGTAGGATg                   <  7:111545/130‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                    TCGCGTTCAGCAGTACCTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTATTGCAGTTGGTCGGATGGGGCAACTCGGGCTTGG  >  NC_002947/6069415‑6069702

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