Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,654,138	G→C	10.8%	A877A (GCC→GCG)	PP_1449 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,654,138	0	G	C	10.8%	190.2 / 9.3	64	A877A (GCC→GCG)	PP_1449	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/30);  new base C (2/5);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

ACCTCGACCTTTGCTGCGGACAACTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTTGTTCATGTTCATGGGTATCCTGAGCTGCACGGTTGTCGATGCGTTTGGCGGTCA  >  NC_002947/1653992‑1654286                                                                                                                                                   |                                                                                                                                                      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                                                                                                              <  5:85984/149‑1 (MQ=255)                                        aTTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAg                                                                                                              <  6:302464/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                                                     <  5:64042/122‑1 (MQ=255)                                                                                                                           gggCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTtgtt                                                     >  6:64042/1‑122 (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)                                                                                                                                                 ggggCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTTGTTCATGTTCATGGGTATCCTGAGCTGCACGGTTGTCGATGCGTTTGGCGGt    >  1:193861/1‑149 (MQ=255)                                                                                                                                                  gggCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAgtcgtc                                                                <  2:311248/88‑1 (MQ=255)                                                                                                                                                  gggCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAgtcgtc                                                                >  1:311248/1‑88 (MQ=255)                                                                                                                                                   ggCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTTGTTCATGTTCATGGGTATCCTGAGCTGCACGGTTGTCGATGCGTTTGGCGGTCa  <  5:160633/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                      ACCTCGACCTTTGCTGCGGACAACTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTTGTTCATGTTCATGGGTATCCTGAGCTGCACGGTTGTCGATGCGTTTGGCGGTCA  >  NC_002947/1653992‑1654286

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