Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,982,841	A→C	5.6%	*120C (TGA→TGC)	PP_3510 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,982,841	0	A	C	5.6%	225.2 / 2.2	71	*120C (TGA→TGC)	PP_3510	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/36);  new base C (2/2);  total (33/38)
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.73e-01

AGGCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCAT  >  NC_002947/3982693‑3982988                                                                                                                                                     |                                                                                                                                                     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                                                                                            <  2:274979/149‑1 (MQ=255)                                                          aCCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTc                                                                                            <  6:218480/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                                                                >  7:195656/1‑149 (MQ=255)                                                                                      cAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTgg                                                                >  3:87337/1‑149 (MQ=255)                                                                                            agagGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGt                                                          <  8:195656/149‑1 (MQ=255)                                                                                             gagGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTg                                                         >  1:87658/1‑149 (MQ=255)                                                                                                gCCTGTGTACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATa                                                                                                                                   <  6:12716/72‑1 (MQ=255)                                                                                                gCCTGTGTACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATa                                                                                                                                   >  5:12716/1‑72 (MQ=255)                                                                                                gCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGt                                                          <  2:134836/145‑1 (MQ=255)                                                                                                gCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGt                                                          >  1:134836/1‑145 (MQ=255)                                                                                                gCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTgcg                                                       >  1:179642/1‑148 (MQ=255)                                                                                                gCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTgcg                                                       >  2:11944/1‑148 (MQ=255)                                                                                                  cTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTAt                                                                        <  3:292155/129‑1 (MQ=255)                                                                                                  cTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTAt                                                                        >  4:292155/1‑129 (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                         >  2:247820/1‑149 (MQ=255)                                                                                                                             gggAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCAc                         >  8:314208/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       >  6:282601/1‑149 (MQ=255)                                                                                                                                               gCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGgt       <  5:282601/149‑1 (MQ=255)                                                                                                                                                  ttGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGgtgt     >  7:72074/1‑148 (MQ=255)                                                                                                                                                  ttGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTc    <  6:409791/149‑1 (MQ=255)                                                                                                                                                    gACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCa   <  6:390893/148‑1 (MQ=255)                                                                                                                                                     aCCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCAt  >  3:168279/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                     AGGCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTG‑ACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTGTCAT  >  NC_002947/3982693‑3982988

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