Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,982,835	G→T	6.2%	G118G (GGG→GGT)	PP_3510 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,982,835	0	G	T	6.2%	201.7 / 3.2	64	G118G (GGG→GGT)	PP_3510	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (29/31);  new base T (2/2);  total (31/33)
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.65e-01

AGGCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGT  >  NC_002947/3982693‑3982983                                                                                                                                               |                                                                                                                                                      aGGCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGa                                                                                                                                                 >  8:113470/1‑149 (MQ=255)          tGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCTGGGACTGGGGAAAACGCCTGCAGGGTCAATGCCCCTGcag                                                                                                                                         <  6:319290/148‑3 (MQ=255)             tGGCACTGGCCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGttt                                                                                                                                      >  7:328451/1‑148 (MQ=255)             tGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGTCAATGCCCCTGca                                                                                                                                          <  8:17777/144‑2 (MQ=255)             tGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCAGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGTCAATGCCCCTGca                                                                                                                                          >  7:17777/1‑143 (MQ=255)              ggCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCTCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGAAGGGGCATTGACCCTGCGGTTTG‑a                                                                                                                                   <  2:239982/149‑1 (MQ=255)               gCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑a                                                                                                                                   <  2:31321/148‑1 (MQ=255)                cACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑Ac                                                                                                                                  <  6:36318/148‑1 (MQ=255)                 aCTGACCTGCCAGGCATTGGACCACGCCGGAATACAACTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCAAGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGa                                                                                                                                >  1:274979/1‑149 (MQ=255)                  cTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGAt                                                                                                                               <  5:262839/149‑1 (MQ=255)                     gccTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATa                                                                                                                              <  8:328451/146‑1 (MQ=255)                     aCCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATa                                                                                                                              <  1:87763/147‑1 (MQ=255)                          ccAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGc                                                                                                                       >  5:409791/1‑149 (MQ=255)                            aGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCtt                                                                                                                     >  2:207514/1‑149 (MQ=255)                             ggCATTGGACCACGCCGGAATACAGGTGGACCATCGCGCTTTGCATGTACGCTATGAGAGCTAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGAGGTTTG‑ACGATAGTATCGCTTc                                                                                                                    <  6:393564/149‑1 (MQ=255)                                aTTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCa                                                                                                                 >  7:275966/1‑149 (MQ=255)                                aTTGGACCAAGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTAAGCTATGACAGAGAAGAGGAATGTGCACGTGTGCATGGCGGCGACTCCGGAAAACGCCTTCAGGGGCATTGACCATGCGGATTG‑ACGATAGTATCGCTTCGCa                                                                                                                 >  1:120377/1‑149 (MQ=255)                                   ggACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCAcc                                                                                                               <  2:73718/148‑1 (MQ=255)                                       cACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTAt                                                                                                           >  5:218480/1‑148 (MQ=255)                                            cGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTc                                                                                                     <  3:258205/149‑1 (MQ=255)                                                    aGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTgg                                                                                             <  7:336515/149‑1 (MQ=255)                                                    aGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTgg                                                                                             <  4:350556/149‑1 (MQ=255)                                                      cTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGcc                                                                                           <  7:118664/149‑1 (MQ=255)                                                         gACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCAAGGGTGCATGGCGGGGACTGGCGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACTATAGTATAGCTTCGCACCATATCAATTCGTGACTGGCCGt                                                                                         >  8:390695/1‑148 (MQ=255)                                                          aCCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTc                                                                                       <  6:218480/149‑1 (MQ=255)                                                          aCCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTc                                                                                       <  2:274979/149‑1 (MQ=255)                                                           ccATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGGGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTcc                                                                                      >  1:285110/1‑149 (MQ=255)                                                                ccGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTTGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATc                                                                                 <  7:113470/149‑1 (MQ=255)                                                                    ttGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGc                                                                              >  5:390893/1‑148 (MQ=255)                                                                          tGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGAt                                                                       <  6:169567/149‑1 (MQ=255)                                                                            tACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGTCTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGt                                                                     <  6:223060/149‑1 (MQ=255)                                                                               gCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTAtt                                                                  <  1:207514/149‑1 (MQ=255)                                                                               gCTATGACAGCGAAGAGGCCTGTGAACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTAt                                                                   >  1:52074/1‑148 (MQ=255)                                                                                   tGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAgg                                                              <  4:177159/149‑1 (MQ=255)                                                                                      cAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTgg                                                           >  3:87337/1‑149 (MQ=255)                                                                                      cAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTgg                                                           >  7:195656/1‑149 (MQ=255)                                                                                            agagGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGt                                                     <  8:195656/149‑1 (MQ=255)                                                                                             gagGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTg                                                    >  1:87658/1‑149 (MQ=255)                                                                                                gCCTGTGTACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATa                                                                                                                              >  5:12716/1‑72 (MQ=255)                                                                                                gCCTGTGTACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATa                                                                                                                              <  6:12716/72‑1 (MQ=255)                                                                                                gCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGt                                                     <  2:134836/145‑1 (MQ=255)                                                                                                gCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGt                                                     >  1:134836/1‑145 (MQ=255)                                                                                                gCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTgcg                                                  >  2:11944/1‑148 (MQ=255)                                                                                                gCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTgcg                                                  >  1:179642/1‑148 (MQ=255)                                                                                                  cTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTAt                                                                   >  4:292155/1‑129 (MQ=255)                                                                                                  cTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTAt                                                                   <  3:292155/129‑1 (MQ=255)                                                                                                     tgCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGTGCATTGACCCTGCGGTTTG‑AAGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCatat                                             <  7:390695/148‑1 (MQ=255)                                                                                                        aCGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGAATATAAca                                         <  2:285110/149‑1 (MQ=255)                                                                                                           ggTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGTGCTTTGACCCTGAGGTTTG‑ACGATAGTATCGCTTCGCACGCTATCAATTCGTGACTGGCCGTTCGCTATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCTTATATcagca                                      <  2:87658/149‑1 (MQ=255)                                                                                                            gTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAAcagcag                                     >  1:22294/1‑149 (MQ=255)                                                                                                                        gACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCg                         >  1:248618/1‑149 (MQ=255)                                                                                                                             gggAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCAc                    >  8:314208/1‑149 (MQ=255)                                                                                                                             gggAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCAc                    >  2:247820/1‑149 (MQ=255)                                                                                                                               gAAAAGGCATGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAAcagca                                      <  4:110284/129‑1 (MQ=255)                                                                                                                               gAAAACGCCTGCAGGGGCATTGACCCTGCGGTTAG‑ACGATAGTATCGCTTCGCACCCTATAAATTCGTGACTGGCCGTTACCGATCCGCCAAGGATGTATTAAGGTTGCATGGTGCGCATATAAcagca                                      >  3:110284/1‑129 (MQ=255)                                                                                                                                  aaCGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGc                >  3:39588/1‑148 (MQ=255)                                                                                                                                   aCGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCaa              <  8:163607/149‑1 (MQ=255)                                                                                                                                      accgCAGGGTCAATGCCCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGa            >  8:184972/4‑148 (MQ=255)                                                                                                                                         gCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATc                          >  8:128380/1‑131 (MQ=255)                                                                                                                                         gCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATc                          <  7:128380/131‑1 (MQ=255)                                                                                                                                          cAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTg                 >  6:386121/1‑139 (MQ=255)                                                                                                                                          cAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTg                 <  5:386121/139‑1 (MQ=255)                                                                                                                                          cAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACaa         <  2:22294/147‑1 (MQ=255)                                                                                                                                          cAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAAt        >  8:236624/1‑148 (MQ=255)                                                                                                                                            ggggCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTACCGATCCGCCAAGGATGTATTAAGGTGGCATTGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAg     >  1:313808/1‑149 (MQ=255)                                                                                                                                             gggCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAgg    <  8:275966/149‑1 (MQ=255)                                                                                                                                             gggCAATGACCCGGCGGTTTGCACCA‑CTTATCGCTTCGCACCCCATCCAATCGGGACTGGCCGTGCCCGCTCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAgg    <  2:120377/149‑1 (MQ=255)                                                                                                                                               gCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGgt  <  5:282601/149‑1 (MQ=255)                                                                                                                                               gCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGgt  >  6:282601/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                      AGGCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGT  >  NC_002947/3982693‑3982983

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