Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	345,710	C→A	7.2%	L713L (CTG→CTT)	PP_0287 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	345,710	0	C	A	7.2%	217.9 / 4.1	69	L713L (CTG→CTT)	PP_0287	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (32/32);  new base A (2/3);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

TCATACCAGTCCAGCACTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGCCGGGTCCGGCATG  >  NC_002947/345562‑345848                                                                                                                                                     |                                                                                                                                           tCATACCAGTCCAGCACTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCa                                                                                                                                           <  1:1016220/150‑1 (MQ=255) tCATACCAGTCCAGCACTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCa                                                                                                                                           <  2:1002432/150‑1 (MQ=255)             aGCACTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCaa                                                                                                                               <  1:292272/150‑1 (MQ=255)             aGCACTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGCCAAGCTGGCCGGCaa                                                                                                                               <  1:531791/150‑1 (MQ=255)              gCACTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAAc                                                                                                                              >  1:542250/1‑150 (MQ=255)                 cTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCtt                                                                                                                           >  1:380177/1‑150 (MQ=255)                  tGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTg                                                                                                                          <  2:1121879/150‑1 (MQ=255)                            aaCTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGt                                                                                                                >  2:952098/1‑150 (MQ=255)                              cTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGccc                                                                                                                       >  1:947587/1‑141 (MQ=255)                              cTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGccc                                                                                                                       <  2:947587/141‑1 (MQ=255)                                gCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGCTCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTg                                                                                                            <  2:221333/150‑1 (MQ=255)                                gCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTg                                                                                                            <  2:871415/150‑1 (MQ=255)                                           aTCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCGGTTCTGGTCGAggcggc                                                                                                 >  2:52226/1‑150 (MQ=255)                                            tCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCa                                                                                                >  2:425221/1‑150 (MQ=255)                                                ttGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGcc                                                                                            >  1:5832/1‑150 (MQ=255)                                                    aCAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCttt                                                                                         <  1:211838/149‑1 (MQ=255)                                                    aCAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCttt                                                                                         >  2:211838/1‑149 (MQ=255)                                                        ccccTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGcc                                                                                    >  1:600724/1‑150 (MQ=255)                                                        ccccTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGcc                                                                                    >  1:685930/1‑150 (MQ=255)                                                          ccTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGa                                                                                  >  2:276055/1‑150 (MQ=255)                                                                  gcgTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGc                                                                          <  1:502299/150‑1 (MQ=255)                                                                  gcgTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGc                                                                          <  1:106736/150‑1 (MQ=255)                                                                     tCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGa                                                                       <  1:383715/150‑1 (MQ=255)                                                                        ttCAGTTCCGGGCTCACTTTGGCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCa                                                                    >  2:542562/1‑150 (MQ=255)                                                                         tCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCCGATAAACGTCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCAc                                                                   >  1:347042/1‑150 (MQ=255)                                                                         tCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCAc                                                                   <  2:600724/150‑1 (MQ=255)                                                                             ttCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAg                                                                                 >  2:984719/1‑132 (MQ=255)                                                                             ttCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAg                                                                                 <  1:984719/132‑1 (MQ=255)                                                                             ttCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCa                                                               >  1:427771/1‑150 (MQ=255)                                                                                 gggCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAgg                                                                                                     >  1:670940/1‑108 (MQ=255)                                                                                 gggCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAgg                                                                                                     <  2:670940/108‑1 (MQ=255)                                                                                     tCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAggg                                                       <  1:952098/150‑1 (MQ=255)                                                                                       aCTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAggggg                                                     >  2:202878/1‑150 (MQ=255)                                                                                        cTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGa                                                    >  2:249313/1‑150 (MQ=255)                                                                                         tttGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGAc                                                   <  1:276055/150‑1 (MQ=255)                                                                                             tctcCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTATCAGGCTGTCGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGGCGAGGCGGCAGGCCTTTGCGCCGAGTCACATCGGACCGCGCCCGCGCGGGGGGTCTTCa                                               >  2:140629/1‑150 (MQ=255)                                                                                               tcCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGGCCACGGCAGCGCAGGGGGACTTCAac                                             >  2:258610/1‑150 (MQ=255)                                                                                                ccAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGAGGACTTCAaca                                            >  1:597201/1‑150 (MQ=255)                                                                                                 cAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGAc                                                   >  2:859837/1‑142 (MQ=255)                                                                                                 cAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGAc                                                   <  1:859837/142‑1 (MQ=255)                                                                                                 cAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAacac                                           >  1:95642/1‑150 (MQ=255)                                                                                                       ttacGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACAt                                     >  2:835397/4‑150 (MQ=255)                                                                                                            tCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAAcgt                                >  2:667591/1‑150 (MQ=255)                                                                                                            tCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAAcgt                                >  2:657586/1‑150 (MQ=255)                                                                                                            tCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGACTTTGCGCCTATGTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAAcgt                                <  2:542250/150‑1 (MQ=255)                                                                                                               aGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAAcgttcg                             <  1:873405/150‑1 (MQ=255)                                                                                                                gTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAAcgttcgt                            <  1:241466/150‑1 (MQ=255)                                                                                                                 tttttCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAAcgttcgtt                           <  1:175064/150‑1 (MQ=255)                                                                                                                 tttttCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAAcgttcgtt                           <  1:1002931/150‑1 (MQ=255)                                                                                                                 tttttCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAAcgttcgtt                           <  2:958226/150‑1 (MQ=255)                                                                                                                 tttttCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAAcgttcgtt                           <  2:347042/150‑1 (MQ=255)                                                                                                                   tttCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGa                         <  1:518604/150‑1 (MQ=255)                                                                                                                       aTCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTg                     >  2:203969/1‑150 (MQ=255)                                                                                                                        tCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTgg                    <  1:249313/150‑1 (MQ=255)                                                                                                                        tCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTgg                    <  1:542562/150‑1 (MQ=255)                                                                                                                        tCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTgg                    >  1:511074/1‑150 (MQ=255)                                                                                                                        tCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGAGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTgg                    <  1:202878/150‑1 (MQ=255)                                                                                                                              ttGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGcc              >  1:267557/1‑150 (MQ=255)                                                                                                                              ttGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGcc              >  2:840983/1‑150 (MQ=255)                                                                                                                              ttGTCTTTCAGGCGGTTGCCGGCCAGCTGGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGcc              <  1:195647/150‑1 (MQ=255)                                                                                                                                gTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGCCgg            >  1:405128/1‑150 (MQ=255)                                                                                                                                 tCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGCCggg           >  2:751159/1‑150 (MQ=255)                                                                                                                                         gCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGCCGGGTCCGGCAt   <  2:267557/150‑1 (MQ=255)                                                                                                                                         gCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGCCGGGTCCGGCAt   <  1:751159/150‑1 (MQ=255)                                                                                                                                          cGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGCCGGGTCCGGCATg  <  2:899717/150‑1 (MQ=255)                                                                                                                                            gTTGCCGCCAAGCTGGCCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGt                                                                                >  2:168309/1‑70 (MQ=37)                                                                                                                                            gTTGCCGCCAAGCTGGCCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGt                                                                                <  1:168309/70‑1 (MQ=37)                                                                                                                                            gTTGCCGCCAAGCTGGCCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCa                                                               >  2:655271/1‑87 (MQ=38)                                                                                                                                            gTTGCCGCCAAGCTGGCCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCa                                                               <  1:655271/87‑1 (MQ=38)                                                                                                                                                     |                                                                                                                                           TCATACCAGTCCAGCACTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGCCGGGTCCGGCATG  >  NC_002947/345562‑345848

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑