Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	159,007	C→G	6.1%	R461P (CGC→CCC) ‡	PP_0149 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	159,007	0	C	G	6.1%	207.9 / 4.3	66	R461P (CGC→CCC) ‡	PP_0149	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (35/27);  new base G (2/2);  total (37/29)
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.82e-01

CTTGGCGAACCAGCTACGGCGCTGCTCGGCATTGAGGGCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTGGATTGGCGGTGCCGCTGCGCAGCAGCTTCAACGACT  >  NC_002947/158865‑159151                                                                                                                                               |                                                                                                                                                 cTTGGCGAACCAGCTACGGCGCTGCTCGGCATTGAGGGCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACg                                                                                                                                            >  6:186310/1‑149 (MQ=255)     gCGAACCAGCTACGGCGCTGCTCGGCATTGAGGGCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTc                                                                                                                                        <  5:144502/149‑1 (MQ=255)      cGAACCAGCTACGGCGCTGCTCGGCATTGAGGGCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGTTACGCCTGGGCTACCGCTTcg                                                                                                                                       <  4:399324/149‑1 (MQ=255)                  ggCGCTGCTCGGCATTGAGGGCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATc                                                                                                                            >  3:113531/1‑148 (MQ=255)                     gctgctCGGCATTGAGGTCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTGAGCACCTCGAACATGGCCTGGTACCGGTCTCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGc                                                                                                                        <  5:285208/149‑1 (MQ=255)                       tgctCGGCATTGAGGGCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGTTACGCCTGGGCTACCGCTTcga                                                                                                                                      >  6:69347/1‑132 (MQ=255)                       tgctCGGCATTGAGGGCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGTTACGCCTGGGCTACCGCTTcga                                                                                                                                      <  5:69347/133‑2 (MQ=255)                        gctCGGCATTGAGGGCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTAACGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCAc                                                                                                                      >  3:240617/1‑148 (MQ=255)                            ggCATTGAGGGCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTg                                                                                                                  <  3:268835/148‑1 (MQ=255)                                     ggCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGc                                                                                                          >  1:118532/1‑147 (MQ=255)                                     ggCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCa                                                                                                        <  8:8304/149‑1 (MQ=255)                                      gccagccaGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCt                                                                                                                   >  8:369444/1‑137 (MQ=255)                                      gccagccaGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCt                                                                                                                   <  7:369444/137‑1 (MQ=255)                                           ccagTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTcccc                                                                                                  >  5:63925/1‑149 (MQ=255)                                             agTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTg                                                                                                <  6:63925/149‑1 (MQ=255)                                               tCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGcc                                                                                              <  5:310451/149‑1 (MQ=255)                                                 ttCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCt                                                                                            >  4:109439/1‑149 (MQ=255)                                                    aaCCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCAGCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGa                                                                                         <  4:369247/149‑1 (MQ=255)                                                    aaCCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGa                                                                                         <  3:109439/149‑1 (MQ=255)                                                       cAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCa                                                                                      >  3:24187/1‑149 (MQ=255)                                                        agagCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCAc                                                                                     <  4:113531/149‑1 (MQ=255)                                                          agCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACgg                                                                                   >  5:402625/1‑149 (MQ=255)                                                            cttgcttgATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGc                                                                                  >  6:313701/1‑148 (MQ=255)                                                                 ttgATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTg                                                                                                           >  2:55562/1‑118 (MQ=255)                                                                 ttgATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTg                                                                                                           <  1:55562/118‑1 (MQ=255)                                                                 ttgATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAAACGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCa                                                                             >  1:116675/1‑148 (MQ=255)                                                                     ttGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGTAAACCg                                                                        >  6:146712/1‑149 (MQ=255)                                                                        ggaaCTGCCTTCTCCCCCCCTCGCACATGGCCGGGTACCGGTCCCCGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTg                                                                     <  3:26047/147‑1 (MQ=255)                                                                        gTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTg                                                                     >  5:154628/1‑149 (MQ=255)                                                                             tGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCat                                                                >  7:100006/1‑149 (MQ=255)                                                                               cTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCatga                                                              <  8:57030/149‑1 (MQ=255)                                                                               cTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCatga                                                              >  6:390155/1‑149 (MQ=255)                                                                               cTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCatg                                                               <  6:402625/148‑1 (MQ=255)                                                                                   ctcAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGcc                                                                                 <  1:21323/126‑1 (MQ=255)                                                                                   ctcAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGcc                                                                                 >  2:21323/1‑126 (MQ=255)                                                                                   ctcAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCt                                                          >  5:411446/1‑149 (MQ=255)                                                                                       gCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTc                                                       >  6:92395/1‑148 (MQ=255)                                                                                           cTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGt                                                  >  1:156265/1‑149 (MQ=255)                                                                                               aaCATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGa                                                                                         >  8:161209/1‑106 (MQ=255)                                                                                               aaCATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGa                                                                                         <  7:161209/106‑1 (MQ=255)                                                                                               aaCATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGtt                                              <  2:306710/149‑1 (MQ=255)                                                                                                  aTGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAg                                           >  2:289802/1‑149 (MQ=255)                                                                                                  aTGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAg                                           <  8:100006/149‑1 (MQ=255)                                                                                                       cTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTg                                      >  7:165348/1‑149 (MQ=255)                                                                                                          gTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTgg                                     <  5:186310/147‑1 (MQ=255)                                                                                                            aCCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAg                                                    <  7:328177/130‑1 (MQ=255)                                                                                                            aCCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAg                                                    >  8:328177/1‑130 (MQ=255)                                                                                                             ccGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTGGATTg                                 <  2:156265/148‑1 (MQ=255)                                                                                                                gTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGt                                               >  6:2553/1‑131 (MQ=255)                                                                                                                   gcgGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTGGATTGGCGGTGc                          <  6:154628/149‑1 (MQ=255)                                                                                                                      gTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGc                                                                                                                        >  3:90244/1‑52 (MQ=255)                                                                                                                      gTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGc                                                                                                                        <  4:90244/52‑1 (MQ=255)                                                                                                                       tAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTGGATTGGCGGTGCCGCt                      >  5:57731/1‑149 (MQ=255)                                                                                                                        aaGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTGGATTGGCGGTGCCGCTg                     >  1:373355/1‑149 (MQ=255)                                                                                                                          gAATCGAAGCGGTAGACCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATTATGATGTCCAGGTCGTTGAGGAGTTGATTGGCGGTTCCGCt                      >  6:286808/1‑146 (MQ=255)                                                                                                                           aaTCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTGGATTGGCGGTGCCGCTgcg                   >  2:66090/1‑148 (MQ=255)                                                                                                                           aaTCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTGGATTGGCGGTGCCGCTgcg                   >  3:312269/1‑148 (MQ=255)                                                                                                                            cgcGAAGCGTTACGCCTGGGCTACCGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTGGATTGGCGGTGCCGCTgcg                   >  4:37869/3‑147 (MQ=255)                                                                                                                             tCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTGGATTGGCGGTGCCGCTGCgcag                >  4:258914/1‑149 (MQ=255)                                                                                                                              cGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTGGATTGGCGGTGCCGCTGCgcag                >  8:47666/1‑148 (MQ=255)                                                                                                                                   cGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTGGATTGGCGGTGCCGCTGCGCAGCAGCtt          <  5:313701/149‑1 (MQ=255)                                                                                                                                      tAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGTTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGAAACCCGTTTTGCATGATGCTGTCCAGGTCGTTGAGGCGTGGATTGGCGGTGCCGCTGCGCAGCAGCTTaaa       <  5:286808/149‑1 (MQ=255)                                                                                                                                        gCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTGGATTGGCGGTGCCGCTGCGCAGCAGCTTCAAc      <  7:47666/148‑1 (MQ=255)                                                                                                                                           cAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCATGCCCTTGAGCACGCCCAGAAAACCGTTGGGCATGATGCTGTCAAGGGCGTTGAGGCGTGAATTGGCGGTCCCGCTGCCAGGAAGCTTCAACGAcg  >  4:207292/1‑148 (MQ=255)                                                                                                                                            aGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTGGATTGGCGGTGCCGCTGCGCAGCAGCTTCa        >  6:418557/1‑142 (MQ=255)                                                                                                                                            aGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTGGATTGGCGGTGCCGCTGCGCAGCAGCTTCa        <  5:418557/142‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                 CTTGGCGAACCAGCTACGGCGCTGCTCGGCATTGAGGGCCAGCCAGTCTTCAACCAGAGCTTGCTTGATTGGTAACTGCTTTCTCAGCACCTCGAACATGGCCTGGTACCGGTCGCGGTAAGAATCGAAGCGGTAGCCCAGGCGTAACGCTTCGCGCGGGTTATCAAGCACGCTGGTATCTGCCAGTCCCCTGCCCTTGAGCACGGCCAGCAAACCGTTGGGCATGATGCTGTCCAGGTCGTTGAGGCGTGGATTGGCGGTGCCGCTGCGCAGCAGCTTCAACGACT  >  NC_002947/158865‑159151

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