Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,891,038	A→T	6.1%	A180A (GCT→GCA)	PP_1696 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,891,038	0	A	T	6.1%	200.7 / 3.3	65	A180A (GCT→GCA)	PP_1696	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (30/31);  new base T (2/2);  total (32/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 = 8.79e-01

GCCGAACACCACCGCGCTCAGCACCCCTCCCACGGTAAAGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGCACATGGTCGTTCAGTTCAAGCATGGCCACGCCCAGCAACGGGCCGGCGGCGTAGCCGACGTTGCTGAGGGTGTACTTGATGGCGAATACCTCGGCCCGC  >  NC_002947/1890891‑1891182                                                                                                                                                    |                                                                                                                                                 gCCGAACACCACCGCGCTCAGCACCCCTCCCACGGTAAAGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAg                                                                                                                                                 >  7:138434/1‑149 (MQ=255)  ccGAACACCACCGCGCTCAGCACCCCTCCCACGGTAAAGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGc                                                                                                                                                <  3:157067/149‑1 (MQ=255)     aaCACCACCGCGCTCAGCACCCCTCCCACGGTAAAGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAg                                                                                                                                             >  8:39838/1‑149 (MQ=255)      acacCACCGCGCTCAGCACCCCTCCCACGGTAAAGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATATGCCAgg                                                                                                                                            >  4:22159/1‑149 (MQ=255)      acacCACCGCGCTCAGCACCCCTCCCACGGTAAAGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATATGCCAgg                                                                                                                                            >  4:27179/1‑149 (MQ=255)      acacCACCGCGCTCAGCACCCCTCCCACGGTAAAGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAgg                                                                                                                                            <  5:188003/149‑1 (MQ=255)         ccaccGCGCTCAGCACCCCTCCCACGGTGAAGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGca                                                                                                                                          >  2:7315/1‑148 (MQ=255)          caccGCGCTCAGCACCCCTCCCACGGTAAAGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAgg                                                                                                                                            >  1:31434/1‑145 (MQ=255)          caccGCGCTCAGCACCCCTCCCACGGTAAAGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAgg                                                                                                                                            <  2:31434/145‑1 (MQ=255)          caccGCGCTCAGCACCCCTCCCACGGTAAAGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGcaca                                                                                                                                        >  4:360798/1‑149 (MQ=255)          caccGCGCTCAGCACCCCTCCCACGGTAAAGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGcaca                                                                                                                                        <  3:362143/149‑1 (MQ=255)          caccGCGCTCAGCACCCCTCCCACGGTAAAGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGcaca                                                                                                                                        >  3:270317/1‑149 (MQ=255)          caccGCGCTCAGCACCCCTCCCACGGTAAAGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGcaca                                                                                                                                        >  3:234836/1‑149 (MQ=255)                  tCAGCACCCCTCCCACGGTAAAGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGc                                                                                                                                           <  4:258689/138‑1 (MQ=255)                  tCAGCACCCCTCCCACGGTAAAGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGc                                                                                                                                           >  3:258689/1‑138 (MQ=255)                    aGCACCCCTCCCACGGTAAAGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATggccaggc                                                                                                                               >  8:82920/1‑148 (MQ=255)                        ccccTCCCACGGTAAAGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATggc                                                                                                                                    <  7:164217/139‑1 (MQ=255)                        ccccTCCCACGGTAAAGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATggc                                                                                                                                    >  8:164217/1‑139 (MQ=255)                         cccTCCCACGGTAAAGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATATGCCAGGCACATGGCCAGGCCAagca                                                                                                                         <  3:22159/149‑1 (MQ=255)                         cccTCCCACGGTAAAGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATATGCCAGGCACATGGCCAGGCCAagca                                                                                                                         <  3:27179/149‑1 (MQ=255)                           cTCCCACGGTAAAGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAagcagc                                                                                                                       <  7:71596/149‑1 (MQ=255)                            tCCCACGGTAAAGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAgcg                                                                                                                      >  4:293317/1‑149 (MQ=255)                             cccACGGTAAAGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGTGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAgcgc                                                                                                                     <  3:357930/149‑1 (MQ=255)                              ccACGGTAAAGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATggccagg                                                                                                                                >  5:58187/1‑137 (MQ=255)                              ccACGGTAAAGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATggccagg                                                                                                                                <  6:58187/137‑1 (MQ=255)                              ccACGGTAAAGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAgcg                                                                                                                      >  5:58569/1‑147 (MQ=255)                              ccACGGTAAAGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCa                                                                                                                    <  4:270317/149‑1 (MQ=255)                                  ggTAAAGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAgg                                                                                                                >  3:273435/1‑149 (MQ=255)                                      aaGCACACCAGCCGACGATCCCGCGTCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGa                                                                                                             <  5:227580/148‑1 (MQ=255)                                      aaGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGa                                                                                                             <  1:53689/148‑1 (MQ=255)                                      aaGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAg                                                                                                            <  1:7315/149‑1 (MQ=255)                                                   gacgaTCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGCACATgg                                                                                               <  7:39838/149‑1 (MQ=255)                                                     cgaTCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGCCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGCACATGGTc                                                                                             <  7:180231/149‑1 (MQ=255)                                                           cgcgCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAAc                                                                                                        >  6:132347/1‑132 (MQ=255)                                                           cgcgCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAAc                                                                                                        <  5:132347/132‑1 (MQ=255)                                                                cAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGCACATGGTCgttcagttc                                                                                    >  6:357339/1‑147 (MQ=255)                                                                 aGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGCACATgg                                                                                               <  4:79741/135‑1 (MQ=255)                                                                 aGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGCACATgg                                                                                               >  3:79741/1‑135 (MQ=255)                                                                  gCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGCACATGGTCGTTCAGTTCAAGc                                                                                <  4:234836/149‑1 (MQ=255)                                                                    cccAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGCACATGGTCGTTCAGTTCAAGCAt                                                                              >  2:98947/1‑149 (MQ=255)                                                                     ccAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGCACATGGTCGTTCAGTTCAAGCATg                                                                             <  8:138434/149‑1 (MQ=255)                                                                                 ggCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGCACATGGTCGTTCAGTTCAAGCATGGCCACGCCCAg                                                                  <  7:82920/148‑1 (MQ=255)                                                                                            gccgCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGcaca                                                                                                  >  1:4861/1‑105 (MQ=255)                                                                                            gccgCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGcaca                                                                                                  <  2:4861/105‑1 (MQ=255)                                                                                                  ggTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGCACATGGTCGTTCAGTTCAAGCATGGCCACGCCCAGCAACGGGCCGGCGGCGTa                                                <  4:273435/149‑1 (MQ=255)                                                                                                           ggCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGCACATGGTCGTTCAGTTCAAGCATGGCCACGCCCAGCAACGGGCCGGCGGCGTAGCCGACGtt                                       >  1:137604/1‑149 (MQ=255)                                                                                                             cgcgcTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGCACATGGTCGTTCAGTTCAAGCATGGCCACGCCCAGCAACGGGCCGGCGGCGTAGCCGACGTTGc                                     <  6:58569/149‑1 (MQ=255)                                                                                                                  tGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGCACATGGTCGTTCAGTTCAAGCATGGCCACGCCCAGCAACGGGCCGGCGGCGTAGCCGACGTTGCTGAgg                                >  6:31319/1‑149 (MQ=255)                                                                                                                    gCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGCACATGGTCGTTCAGTTCAAGCATGGCCACGCCCAGCAACGGGCCGGCGGCGTAGCCGACGTTGCTGAgg                                <  3:293317/147‑1 (MQ=255)                                                                                                                    gCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGCACATGGTCGTTCAGTTCAAGCATGGCCACGCCCAGCAACGGGCCGGCGGCGTAGCCGACGTTGCTGAgg                                <  7:115777/147‑1 (MQ=255)                                                                                                                      cTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGCACATGGTCGTTCAGTTCAAGCATGGCCACGCCCAGCAACGGGCCGGCGGCGTAGCCGACGTTGCTGAGGgtg                             >  5:323853/1‑148 (MQ=255)                                                                                                                        gCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGcac                                                                                                   <  7:25273/76‑1 (MQ=255)                                                                                                                        gCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGcac                                                                                                   >  8:25273/1‑76 (MQ=255)                                                                                                                            gCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGCACATGGTCGTTCAGTTCAAGCATGGCCACGCCCAGCAACGGGCCGGCGGCGTAGCCGACGTTGCTGAGGGTGTACTTGa                      <  5:357339/149‑1 (MQ=255)                                                                                                                              tGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGCACATGGTCgttcag                                                                                       <  3:246422/82‑1 (MQ=255)                                                                                                                              tGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGCACATGGTCgttcag                                                                                       >  4:246422/1‑82 (MQ=255)                                                                                                                                     aCCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGCACATGGTCGTTCAGTTCAAGCATGGCCACGCCCAGCAACGGGCCGGCGGCGTAGCCGACGTTGCTGAGGGAGTACTTGATGGCGAAt              <  6:269719/148‑1 (MQ=255)                                                                                                                                     aCCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGCACATGGTCGTTCAGTTCAAGCATGGCCACGCCCAGCAACGGGCAGGCGGCGTAGCCGACGTTGCTGAGGGTGTACTTGATGGCGAATa             >  6:353630/1‑149 (MQ=255)                                                                                                                                         aGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGCACATGGTCGTTCAGTTCAAGCATGGCCACGCCCAGCAACGGGCCGGCGGCGTAGCCGACGTTGCTGAGGGTGTACTTGATGGCGAATACCTc         >  5:85404/1‑149 (MQ=255)                                                                                                                                         aGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAAAGGCACATGGTCGTTCAGTTCAAGCATGGCCACGCCCAGCAACGGGCCGGCGGCGTAGCAGACGTTGCTTAGGGTGTACTTGATGGCGAATACCTc         >  8:343406/1‑149 (MQ=255)                                                                                                                                              ccAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGCACATGGTCGTTCAGTTCAAGCATGGCCACGCCCAGCAACGGGCCGGCGGCGTAGCCGACGTTGCTGAGGGTGTACTTGATGGCGAATACCTCGGccc    >  6:367926/1‑149 (MQ=255)                                                                                                                                               cAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGCACATGGTCGTTCAGTTCAAGCATGGCCACGCCCAGCAACGGGCCGGCGGCGTAGCCGACGTTGCTGAGGGTGTACTTGATGGCGAAt              <  7:261311/138‑1 (MQ=255)                                                                                                                                               cAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGCACATGGTCGTTCAGTTCAAGCATGGCCACGCCCAGCAACGGGCCGGCGGCGTAGCCGACGTTGCTGAGGGTGTACTTGATGGCGAAt              >  8:261311/1‑138 (MQ=255)                                                                                                                                                aaTAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGCACATGGTCGTTCAGTTCAAGCATGGCCACGCCCAGCAACGGGCCGGCGGCGTAGCCGACGTTGCTGAGGGTGTACTTGATGGCGAATACCTCGGCCCGc  >  2:205842/1‑149 (MQ=255)                                                                                                                                                aaTAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGCACATGGTCGTTCAGTTCAAGCATGGCCACGCCCAGCAACGGGCCGGCGGCGTAGCCGACGTTGCTGAGGGTGTACTTGATGGCGAATACCTCGGCCCGc  <  1:205842/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 GCCGAACACCACCGCGCTCAGCACCCCTCCCACGGTAAAGCACACCAGCCGACGATCCCGCGCCAGCCCCAGCGCCACTTGGCCGAAGCCGGCCGCAGGTTTGTCCGGCGCGCTGGCCTGCAGGCTGCGATCACCCAGGCGCCAATAAGCCAGGCACATGGCCAGGCCAAGCAGCGCAGAGGCGAGGAACGGCACATGGTCGTTCAGTTCAAGCATGGCCACGCCCAGCAACGGGCCGGCGGCGTAGCCGACGTTGCTGAGGGTGTACTTGATGGCGAATACCTCGGCCCGC  >  NC_002947/1890891‑1891182

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