Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	345,708	G→C	7.2%	A714G (GCC→GGC)	PP_0287 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	345,708	0	G	C	7.2%	218.0 / 7.5	69	A714G (GCC→GGC)	PP_0287	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/32);  new base C (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 = 6.82e-01

TCGTCATACCAGTCCAGCACTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGCCGGGTCCGGCATG  >  NC_002947/345559‑345848                                                                                                                                                      |                                                                                                                                             tCGTCATACCAGTCCAGCACTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCgg                                                                                                                                              >  2:106736/1‑150 (MQ=255) tCGTCATACCAGTCCAGCACTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCgg                                                                                                                                              >  2:175064/1‑150 (MQ=255) tCGTCATACCAGTCCAGCACTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCgg                                                                                                                                              <  2:25225/150‑1 (MQ=255)    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                                                                                         >  2:211838/1‑149 (MQ=255)                                                       aCAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCttt                                                                                         <  1:211838/149‑1 (MQ=255)                                                           ccccTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGcc                                                                                    >  1:685930/1‑150 (MQ=255)                                                           ccccTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGcc                                                                                    >  1:600724/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                                                                                 <  1:984719/132‑1 (MQ=255)                                                                                ttCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAg                                                                                 >  2:984719/1‑132 (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                           <  2:347042/150‑1 (MQ=255)                                                                                                                    tttttCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAAcgttcgtt                           <  1:1002931/150‑1 (MQ=255)                                                                                                                    tttttCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAAcgttcgtt                           <  2:958226/150‑1 (MQ=255)                                                                                                                    tttttCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAAcgttcgtt                           <  1:175064/150‑1 (MQ=255)                                                                                                                      tttCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGa                         <  1:518604/150‑1 (MQ=255)                                                                                                                          aTCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTg                     >  2:203969/1‑150 (MQ=255)                                                                                                                           tCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTgg                    >  1:511074/1‑150 (MQ=255)                                                                                                                           tCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTgg                    <  1:249313/150‑1 (MQ=255)                                                                                                                           tCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTgg                    <  1:542562/150‑1 (MQ=255)                                                                                                                           tCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGAGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTgg                    <  1:202878/150‑1 (MQ=255)                                                                                                                                 ttGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGcc              >  2:840983/1‑150 (MQ=255)                                                                                                                                 ttGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGcc              >  1:267557/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                                                                                <  1:168309/70‑1 (MQ=37)                                                                                                                                               gTTGCCGCCAAGCTGGCCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGt                                                                                >  2:168309/1‑70 (MQ=37)                                                                                                                                               gTTGCCGCCAAGCTGGCCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCa                                                               <  1:655271/87‑1 (MQ=38)                                                                                                                                               gTTGCCGCCAAGCTGGCCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCa                                                               >  2:655271/1‑87 (MQ=38)                                                                                                                                                      |                                                                                                                                             TCGTCATACCAGTCCAGCACTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGCCGGGTCCGGCATG  >  NC_002947/345559‑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: ‑