Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,863,320	G→C	7.2%	A322P (GCG→CCG)	PP_1666 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,863,320	0	G	C	7.2%	168.0 / 2.0	55	A322P (GCG→CCG)	PP_1666	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (28/23);  new base C (2/2);  total (30/25)
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.97e-01

TGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGG  >  NC_002947/1863174‑1863467                                                                                                                                                   |                                                                                                                                                     tGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCCGGCACcgccgg                                                                                                                                                   >  2:843441/1‑150 (MQ=255) tGGCCCGCGTGCTTGAGCCCTATGGACCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACcggcgg                                                                                                                                                   <  1:714995/150‑1 (MQ=255)   gCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTg                                                                                                                                                 >  1:219141/1‑150 (MQ=255)          tGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCCGGCACCGCCGGTGcctgctc                                                                                                                                          <  1:843441/150‑6 (MQ=255)             ttGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGc                                                                                                                                       <  2:219141/150‑1 (MQ=255)                   ccTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAc                                                                                                                                 >  1:1297965/1‑150 (MQ=255)                    cTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAc                                                                                                                                 <  2:1091483/149‑1 (MQ=255)                    cTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAc                                                                                                                                 >  1:1091483/1‑149 (MQ=255)                    cTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCAGGTCGAGCAGGCACCGGCTGCGTGGGCCATGCCATTAcc                                                                                                                                >  2:307635/1‑150 (MQ=255)                            aatgtttGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTgg                                                                                                                        >  1:410187/6‑150 (MQ=255)                               gCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGACATGCCAGTACCGGCAGTGGAcg                                                                                                                     >  2:1038799/1‑150 (MQ=255)                                cTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGAcgc                                                                                                                    <  2:550470/150‑1 (MQ=255)                                     aGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAg                                                                                                               >  2:1132907/1‑150 (MQ=255)                                         ggcTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCGACGGGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCcgc                                                                                                           <  2:1328415/150‑1 (MQ=255)                                                   aCTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCcgcg                                                                                                 >  2:48340/1‑150 (MQ=255)                                                      aCCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCt                                                                                              >  1:910661/1‑150 (MQ=255)                                                           gACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGCCGGTGcc                                                                                                                                               >  1:950552/1‑96 (MQ=255)                                                           gACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGCCGGTGcc                                                                                                                                               <  2:950552/96‑1 (MQ=255)                                                             ccTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGcagc                                                                                       <  1:307635/150‑1 (MQ=255)                                                             ccTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGcagc                                                                                       <  2:910661/150‑1 (MQ=255)                                                              cTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGcagcc                                                                                      >  2:919780/1‑150 (MQ=255)                                                              cTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGcagcc                                                                                      <  2:1297965/150‑1 (MQ=255)                                                                aCGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGcagccag                                                                                    <  2:265065/150‑1 (MQ=255)                                                                         cGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGc                                                                           <  2:1423846/150‑1 (MQ=255)                                                                          gggCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAGGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGcc                                                                          >  1:803049/1‑150 (MQ=255)                                                                          gggCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGcc                                                                          <  1:1200924/150‑1 (MQ=255)                                                                            gCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGt                                                                        >  2:571810/1‑150 (MQ=255)                                                                            gCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCAGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGt                                                                        >  2:551154/1‑150 (MQ=255)                                                                                  gTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACg                                                                 >  1:430291/1‑150 (MQ=255)                                                                                       cAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCG‑TCGAACGGCTTg                                                            <  1:48340/150‑1 (MQ=255)                                                                                           tACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGtt                                                        >  1:895516/1‑150 (MQ=255)                                                                                                gCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGc                                                           >  1:1166798/1‑142 (MQ=255)                                                                                                gCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGc                                                           <  2:1166798/142‑1 (MQ=255)                                                                                                ctgcAGCTTGGCCCGGACACGCCGCCGGACGTCTCTGCCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGGGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGc                                                   <  2:373435/147‑1 (MQ=255)                                                                                                  ccAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCAAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGTCTTGCGTTTTCACTg                                                 >  2:227054/1‑150 (MQ=255)                                                                                                   cAGCTTGGCCAGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGGGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCtag                                                >  1:373435/1‑148 (MQ=255)                                                                                                       ttGGCCTGGCAAAGCTGGAGGGAGTGCCGGTCTATCATGCACCTGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCAGCGCCGCCCGAACCACCGCAGCTAGCACGCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAg                                            >  2:1020792/1‑150 (MQ=255)                                                                                                         ggcctggccAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGa                                          <  2:895516/150‑1 (MQ=255)                                                                                                          tccaggccAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAg                                         >  2:787882/5‑150 (MQ=255)                                                                                                           cctggccAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAgc                                        >  1:1024745/1‑150 (MQ=255)                                                                                                            ctggccAGGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAgcg                                       <  2:803049/150‑1 (MQ=255)                                                                                                                 cAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGc                                  >  2:1250179/1‑150 (MQ=255)                                                                                                                      tGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGtga                             <  1:919780/150‑1 (MQ=255)                                                                                                                      tGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGtga                             <  2:498909/150‑1 (MQ=255)                                                                                                                        cAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGtgact                           >  2:1068257/1‑150 (MQ=255)                                                                                                                         aGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCa                                                                                                        >  2:32268/1‑73 (MQ=255)                                                                                                                         aGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCa                                                                                                        <  1:32268/73‑1 (MQ=255)                                                                                                                          ggAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGtgactga                         >  2:1399729/1‑150 (MQ=255)                                                                                                                               tGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGc                    <  2:489954/150‑1 (MQ=255)                                                                                                                                  cGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTc                 <  1:1068257/150‑1 (MQ=255)                                                                                                                                   ggTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCg                <  1:253169/150‑1 (MQ=255)                                                                                                                                    gTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGc               >  2:1075077/1‑150 (MQ=255)                                                                                                                                         gCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTgg                                                >  2:842859/1‑112 (MQ=255)                                                                                                                                         gCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTgg                                                <  1:842859/112‑1 (MQ=255)                                                                                                                                          ccagCACCACCGGTGCCGGACATGCCAGTAACGGCATTTGGCGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTGCGCTGGTAATGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATc         >  2:1020715/4‑150 (MQ=255)                                                                                                                                             gCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGc                                  <  2:827904/122‑1 (MQ=255)                                                                                                                                             gCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGc                                  >  1:827904/1‑122 (MQ=255)                                                                                                                                                 cggcggTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGtttt                                                      <  1:1153352/98‑1 (MQ=255)                                                                                                                                                 cggcggTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGtttt                                                      >  2:1153352/1‑98 (MQ=255)                                                                                                                                                 cggcggTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTCCTAAGCAGCGCTTGCCGTGACTGACATGCGTCGCTGGATctggctgg  <  2:410187/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     TGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGG  >  NC_002947/1863174‑1863467

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