Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,069,547	C→A	6.9%	V21L (GTG→TTG) ‡	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,547	0	C	A	6.9%	210.7 / 4.3	67	V21L (GTG→TTG) ‡	PP_5323	M23/M37 family peptidase
Reads supporting (aligned to +/- strand):  ref base C (27/35);  new base A (2/3);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

GTCGAAGGTTTTCGCGTTCAGCAGTACCTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTATTGCAGTTGGTCGGATGGGGCAAC  >  NC_002947/6069404‑6069692                                                                                                                                                |                                                                                                                                                    gTCGAAGGTTTTCGCGTTCAGCAGTACCTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCt                                                                                                                                                >  8:155735/1‑149 (MQ=255)            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)                                       gCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGccc                                                                                                                                                 <  1:283885/110‑1 (MQ=255)                                       gCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGccc                                                                                                                                                 >  2:283885/1‑110 (MQ=255)                                        cGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAAcagca                                                                                                       <  3:36221/149‑1 (MQ=255)                                        cGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAAcagca                                                                                                       <  4:214513/149‑1 (MQ=255)                                              gACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGgc                                                                                                 >  2:219450/1‑149 (MQ=255)                                              gACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGgc                                                                                                 >  1:294686/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                                                                                                   >  4:153701/1‑138 (MQ=255)                                                       cTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAAcagcagcag                                                                                                   <  3:153701/138‑1 (MQ=255)                                                        tCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGc                                                                                                                                             <  8:207037/97‑1 (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)                                                                                                                                                |                                                                                                                                                    GTCGAAGGTTTTCGCGTTCAGCAGTACCTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGA‑GGCCATGAGCAATG‑CACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTATTGCAGTTGGTCGGATGGGGCAAC  >  NC_002947/6069404‑6069692

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