Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,318,939	C→A	7.7%	N400K (AAC→AAA)	PP_3789 →	efflux transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,318,939	0	C	A	7.7%	209.8 / 6.6	65	N400K (AAC→AAA)	PP_3789	efflux transporter
Reads supporting (aligned to +/- strand):  ref base C (33/27);  new base A (2/3);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCA‑CTGAT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGA  >  NC_002947/4318805‑4319082                                                                                                                                             |                                                                                                                                                 gcggcAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCA‑CTGAT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGCGCAAAGTTGCCCggggtca                                                                                                                                    <  4:234636/149‑4 (MQ=255)   ggcAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCA‑CTGAT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGCGCAAAGTTGCCCggg                                                                                                                                        <  2:73863/143‑1 (MQ=255)   ggcAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCA‑CTGAT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGCGCAAAGTTGCCCggg                                                                                                                                        >  1:73863/1‑143 (MQ=255)   ggcAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCA‑CTGAT‑‑CGC‑‑C‑TGCGCGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATc                                                                                                                                 <  1:63885/149‑1 (MQ=255)    gcAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCA‑CTGAT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCt                                                                                                                                >  3:37044/1‑149 (MQ=255)     cAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCA‑CTGAT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCtt                                                                                                                               <  6:96601/149‑1 (MQ=255)      aGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCA‑CTGAT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCttt                                                                                                                              <  1:17709/149‑1 (MQ=255)        tCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCA‑CTGAT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTaa                                                                                                                            <  5:267047/149‑1 (MQ=255)          cggccGGCGTTGTTGGGTATCATACTATCGCTGGGCGGCCCCCCTGGTTTGTGGCTCTCACCTGATCCCGCGACATGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCt                                                                                                                                <  1:106494/146‑1 (MQ=255)               ggCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCA‑CTGAT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCa                                                                                                                     <  2:294726/149‑1 (MQ=255)                     gtgtGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCA‑CTGAT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTAcgc                                                                                                               >  2:147434/1‑149 (MQ=255)                       gtgGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCA‑CTGAT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATc                                                                                                                                 <  1:96624/129‑1 (MQ=255)                       gtgGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCA‑CTGAT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATc                                                                                                                                 >  2:96624/1‑129 (MQ=255)                            ttCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCA‑CTGAT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTgg                                                                                                        <  4:37044/149‑1 (MQ=255)                                tACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCA‑CTGAT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGgtcgt                                                                                                    >  4:123898/1‑149 (MQ=255)                                 aCTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCA‑CTGAT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGgtcgt                                                                                                    >  5:225056/1‑148 (MQ=255)                                 aCTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCA‑CTGAT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGgtcgt                                                                                                    >  2:196874/1‑148 (MQ=255)                                 aCTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCA‑CTGAT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGgtcgt                                                                                                    <  1:196874/148‑1 (MQ=255)                                  cTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCA‑CTGAT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCggg                                                                                                                                        <  2:208117/112‑1 (MQ=255)                                  cTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCA‑CTGAT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCggg                                                                                                                                        >  1:208117/1‑112 (MQ=255)                                         tcgcGGACCCCCAGGGTTTGTGCTTCTCA‑CTGAT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTgg                                                                                                        <  5:168143/136‑1 (MQ=255)                                         tcgcGGACCCCCAGGGTTTGTGCTTCTCA‑CTGAT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTgg                                                                                                        >  6:168143/1‑136 (MQ=255)                                                 ccccAGGGTTTGTGCTTCTCA‑CTGAT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCAc                                                                                   >  6:173916/1‑149 (MQ=255)                                                 ccccAGGGTTTGTGCTTCTCA‑CTGAT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCAc                                                                                   <  3:312343/149‑1 (MQ=255)                                                  cccAGGGTTTGTGCTTCTCA‑CTGAT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGgtcgt                                                                                                    <  3:100319/131‑1 (MQ=255)                                                  cccAGGGTTTGTGCTTCTCA‑CTGAT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGgtcgt                                                                                                    >  4:100319/1‑131 (MQ=255)                                                         tttGTGCTTCTCA‑CTGAT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTgg                                                                                                        >  1:115252/1‑120 (MQ=255)                                                         tttGTGCTTCTCA‑CTGAT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTgg                                                                                                        <  2:115252/120‑1 (MQ=255)                                                                     a‑CTGAT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGtct                                                                                            >  8:294061/1‑120 (MQ=255)                                                                     a‑CTGAT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGtct                                                                                            <  7:294061/120‑1 (MQ=255)                                                                          aT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCt                                                             >  5:116509/1‑147 (MQ=255)                                                                                c‑‑c‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCt                                                        >  3:83310/1‑148 (MQ=255)                                                                                c‑‑c‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCt                                                        >  4:223747/1‑148 (MQ=255)                                                                                   c‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTc                                                                                                                      <  4:279412/85‑1 (MQ=255)                                                                                       cGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATg                                                    <  8:6881/148‑1 (MQ=255)                                                                                         tGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAg                                                 >  3:167220/1‑149 (MQ=255)                                                                                           cTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGtt                                               <  1:147434/149‑1 (MQ=255)                                                                                           cTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGtt                                               <  6:225056/149‑1 (MQ=255)                                                                                              gCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTtata                                            >  4:173860/1‑149 (MQ=255)                                                                                               cGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTtat                                           <  4:83310/149‑1 (MQ=255)                                                                                                  gTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATaaa                                        >  4:51951/1‑149 (MQ=255)                                                                                                     gcagcCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTAcgc                                                                                                               >  4:45860/1‑75 (MQ=255)                                                                                                     gcagcCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATAc                                     >  1:233692/1‑149 (MQ=255)                                                                                                         ccTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCGCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATg                                 >  3:220998/1‑149 (MQ=255)                                                                                                               tgCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCt                            <  2:233692/148‑1 (MQ=255)                                                                                                               tgCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTg                           <  4:295042/149‑1 (MQ=255)                                                                                                                gCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTg                           >  3:73657/1‑148 (MQ=255)                                                                                                                 ccATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTgg                          >  4:123061/1‑148 (MQ=255)                                                                                                                       gcCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCACGCGGCTGGAAATGAc                   >  7:106420/1‑149 (MQ=255)                                                                                                                         ccTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACt                  >  4:271251/1‑148 (MQ=255)                                                                                                                          cTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTAc                >  4:192268/1‑149 (MQ=255)                                                                                                                             gACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGTTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTAcc               <  1:255698/147‑1 (MQ=255)                                                                                                                             gACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTAcc               <  3:223747/147‑1 (MQ=255)                                                                                                                             gACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCtt             >  5:227175/1‑149 (MQ=255)                                                                                                                               cTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTg            >  4:217407/1‑148 (MQ=255)                                                                                                                               cTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTg            >  4:212956/1‑148 (MQ=255)                                                                                                                               cTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTg            >  5:135796/1‑148 (MQ=255)                                                                                                                               atctccccGCGCAAAGTTGCCCGGGGAGGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATg                     >  1:14845/5‑140 (MQ=255)                                                                                                                               atctccccGCGCAAAGTTGCCCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATg                     <  2:14845/135‑1 (MQ=255)                                                                                                                                  aCCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCACACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGtat         <  3:123898/148‑1 (MQ=255)                                                                                                                                   ccccGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACACTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATAc       >  6:34817/1‑149 (MQ=255)                                                                                                                                    cccGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATAcc      >  6:12366/1‑149 (MQ=255)                                                                                                                                    cccGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATAcc      >  4:259452/1‑149 (MQ=255)                                                                                                                                    cccGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATAcc      <  4:73657/149‑1 (MQ=255)                                                                                                                                        ggCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGa  <  6:116509/149‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                 GCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCA‑CTGAT‑‑CGC‑‑C‑TGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGA‑GATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGA  >  NC_002947/4318805‑4319082

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