Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,654,134	G→C	10.2%	R879G (CGC→GGC)	PP_1449 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,654,134	0	G	C	10.2%	198.1 / 7.5	66	R879G (CGC→GGC)	PP_1449	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/32);  new base C (2/5);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.92e-01

CCTTGACCTCGACCTTTGCTGCGGACAACTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTTGTTCATGTTCATGGGTATCCTGAGCTGCACGGTTGTCGATGCGTTT  >  NC_002947/1653987‑1654278                                                                                                                                                    |                                                                                                                                                  ccTTGACCTCGACCTTTGCTGCTGACAACTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAgcgc                                                                                                                                                  <  5:47692/149‑1 (MQ=255) ccTTGACCTCGACCTTTGCTGCGGACAACTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAgcgc                                                                                                                                                  >  2:77065/1‑149 (MQ=255) ccTTGACCTCGACCTTTGCTGCGGACAACTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAgcgc                                                                                                                                                  >  6:346542/1‑149 (MQ=255)     gACCTCGACCTTTGCTGCGGACAACTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCgg                                                                                                                                                <  5:327764/147‑1 (MQ=255)      aCCTCGACCTTTGCTGCGGACAACTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCgggg                                                                                                                                              <  6:364063/148‑1 (MQ=255)          cGACCTTTGCTGCGGACAACTCGCTGACCTGAGCAATTCCACTGCTCTGGTTTTCCAGGCGCTTGAGTTGTTCGACCGTCATGTCGGCCCCCACGCTGGTAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCCCCGCGCTTcc                                                                                                                                         <  8:169199/149‑1 (MQ=255)           gACCTTTGCTGCGGACACCTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCCCCGCGCTTcca                                                                                                                                        <  4:314556/149‑2 (MQ=255)             ccTTTGCTGCGGACAACTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCCCCGCGCTTccagc                                                                                                                                      <  1:36389/149‑4 (MQ=255)                 ttctgcGGACAACTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCCCCGCGCTTccagcaa                                                                                                                                    <  1:369796/145‑6 (MQ=255)                 tgctgcGGACACCTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGGCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCAccc                                                                                                                                   <  3:82732/148‑1 (MQ=255)                 tgctgcGGACAACTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCAcc                                                                                                                                    <  1:179977/147‑1 (MQ=255)                    tgcGGACAACTCGCTGACCTGAGCAATTCCACGGCGCAGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCa                                                                                                                               >  6:150720/1‑149 (MQ=255)                      cGGACAACTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCa                                                                                                                               <  1:325246/147‑1 (MQ=255)                           aaCTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCg                                                                                                                        >  1:88969/1‑149 (MQ=255)                           aaCTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCg                                                                                                                        >  4:209286/1‑149 (MQ=255)                            aCTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGGTTGAGGTGTTCGACCTTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCGTCCCCCAGCACCAGGCGt                                                                                                                       >  8:164320/1‑149 (MQ=255)                             cTCGCTGACCTGAGCAATTCCCCGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTGTTCGGGCTAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGtt                                                                                                                      <  1:77917/149‑1 (MQ=255)                             cTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCTTGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGtt                                                                                                                      >  5:303086/1‑149 (MQ=255)                             cTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGtt                                                                                                                      <  5:346542/149‑1 (MQ=255)                              tCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTc                                                                                                                     >  5:302464/1‑149 (MQ=255)                              tCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTc                                                                                                                     >  5:202400/1‑149 (MQ=255)                                        gAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCgg                                                                                                           >  1:105670/1‑149 (MQ=255)                                         aGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCggg                                                                                                          >  6:203176/1‑149 (MQ=255)                                         aGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCg                                                                                                            <  6:20571/147‑1 (MQ=255)                                         aGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTAGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGGCTTCGGTCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCgg                                                                                                           <  5:67682/148‑1 (MQ=255)                                           cAATTCCTCGGCGCTGGTTCTCCAGGCGCTTGAGGTGTTCGACCGTCATGGCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGATGCTTGCTGGAAGCGCGGGGCTTCCTCACGCAGCACCAGGCGTGCGATAGgggggggc                                                                                                        <  8:290746/149‑3 (MQ=255)                                           cAATTCCACGGCGCTGGTGTTCCAGGCGCGTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGGGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTc                                                                                                        <  8:167927/149‑1 (MQ=255)                                            aagtCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTAGGCCCCCACGCTGGGAGGTTCCATGTCGTCTTCGGGCGAAGCTTGCTGGAAGCCCCGCGCTTcca                                                                                                                                        <  6:104648/113‑2 (MQ=255)                                            aaTTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCATCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGt                                                                                                         >  3:368471/1‑147 (MQ=255)                                            aaTTCCACGGCGCTGGTTTTCAAGGCGGTTGAGGTGTTCGACCGTCATGTCGGCCCCCCCGCTGGGAGGTGCCATGGCGTCGTCGGGCGAAGATTGCTGGAAGCCCCGCGCTTcca                                                                                                                                        >  5:104648/1‑115 (MQ=255)                                             aTTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAg                                                                                                      <  6:302464/149‑1 (MQ=255)                                             aTTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAg                                                                                                      <  5:85984/149‑1 (MQ=255)                                               tCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGt                                                                                                    >  7:94992/1‑149 (MQ=255)                                                     gcgcTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCgg                                                                                               >  2:225331/1‑148 (MQ=255)                                                      cgcTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTa                                                                                             >  8:277443/1‑149 (MQ=255)                                                         tGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTc                                                                                          >  4:252112/1‑149 (MQ=255)                                                               tCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGAGCGGGTCAGGTCGCGGTATTCGATAc                                                                                     <  1:10461/148‑1 (MQ=255)                                                                      gcTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGccc                                                                             <  2:105670/149‑1 (MQ=255)                                                                                ttCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTAccc                                                                    <  6:9010/148‑1 (MQ=255)                                                                                ttCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTAccc                                                                    >  4:259618/1‑148 (MQ=255)                                                                                ttCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCg                                                                   >  3:78052/1‑149 (MQ=255)                                                                                   gACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCTGTATTCGATACTGGCGCCCAGGCTACCCGAc                                                                 <  3:209286/148‑1 (MQ=255)                                                                                   gACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAc                                                                 <  2:88969/148‑1 (MQ=255)                                                                                    aCCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCAGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAc                                                                 >  6:160633/1‑147 (MQ=255)                                                                                      cGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAgt                                                            >  2:204895/1‑149 (MQ=255)                                                                                         cATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAgtcgt                                                         >  3:77659/1‑149 (MQ=255)                                                                                               ggCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGc                                                   >  2:45192/1‑149 (MQ=255)                                                                                                 cccccACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCtt                                                 <  6:303086/149‑1 (MQ=255)                                                                                                   cccACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTtg                                               <  3:252112/149‑1 (MQ=255)                                                                                                    ccACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTtgt                                              <  1:77065/149‑1 (MQ=255)                                                                                                     cACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTtgtt                                             >  5:51403/1‑149 (MQ=255)                                                                                                     cACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTtgt                                              >  5:37864/1‑148 (MQ=255)                                                                                                        gCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTTGTCACGGCTTtgttcat                                          >  8:67128/1‑149 (MQ=255)                                                                                                            ggAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTtgttcatgtt                                       <  3:259618/148‑1 (MQ=255)                                                                                                                gTGCCATGGCGTCTTCTGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTTGTTCATGTTCATgg                                  <  6:37864/149‑1 (MQ=255)                                                                                                                gTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTTGTTCATGTTCATgg                                  <  1:225331/149‑1 (MQ=255)                                                                                                                    cATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGACACAGTCGTCACGGCTTTGTTCATGTTCATGGGTa                               <  1:45192/149‑1 (MQ=255)                                                                                                                          gTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTTGTTCATGTTCATGGGTATCCTGAg                        <  5:150720/149‑1 (MQ=255)                                                                                                                                gggCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTtgtt                                             >  6:64042/1‑122 (MQ=255)                                                                                                                                gggCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTtgtt                                             <  5:64042/122‑1 (MQ=255)                                                                                                                                  gCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTTGTTCATGTTCATGGGTATCCTGAGCTGCACgg                >  6:140834/1‑149 (MQ=255)                                                                                                                                   cGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTTGTTCATGTTCATGGGTATCCTGAGCTGCACGGt               >  3:344415/1‑149 (MQ=255)                                                                                                                                   cGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGC‑TTGTTCATGTTCATGGGTATCCTGAGCTGCACGGtt              >  7:30997/1‑149 (MQ=255)                                                                                                                                   cGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGC‑TTGTTCATGTTCATGGGTATCCTGAGCTGCACGGtt              <  8:30997/149‑1 (MQ=255)                                                                                                                                    gAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTTGTTCATGTTCATGGGTATCCTGAGCTGCACGGtt              <  4:344415/149‑1 (MQ=255)                                                                                                                                            gaggAAGCCCCGCGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTTGTTCATGTTCATGGGTATCCTGAGCTGCACGGTTGTCGATGc      >  7:332864/3‑149 (MQ=255)                                                                                                                                               gAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTTGTTCATGTTCATGGGTATCCTGAGCTGCACGGTTGTCGATGCGt    <  7:277443/148‑1 (MQ=255)                                                                                                                                                aaGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTTGTTCATGTTCATGGGTATCCTGAGCTGCACGGTTGTCGATGCGttt  <  6:202400/149‑1 (MQ=255)                                                                                                                                                aaGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTTGTTCATGTTCATGGGTATCCTGAGCTGCACGGTTGTCGATGCGttt  <  2:222266/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  CCTTGACCTCGACCTTTGCTGCGGACAACTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTTGTTCATGTTCATGGGTATCCTGAGCTGCACGGTTGTCGATGCGTTT  >  NC_002947/1653987‑1654278

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