Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	558,477	A→C	7.7%	E21D (GAA→GAC)	rplR →	50S ribosomal protein L18

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	558,477	0	A	C	7.7%	200.1 / 5.2	74	E21D (GAA→GAC)	rplR	50S ribosomal protein L18
Reads supporting (aligned to +/- strand):  ref base A (32/35);  new base C (3/3);  total (36/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 = 7.95e-01

GCCAGAGCCTTACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGT‑CCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGACGCGGC  >  NC_002947/558329‑558623                                                                                                                                                        |                                                                                                                                                   gCCAGAGCCTTACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGT‑CCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGaa                                                                                                                                                    >  5:1544/1‑149 (MQ=255)    agagCCTTACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGT‑CCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTc                                                                                                                                                 <  5:247978/149‑1 (MQ=255)    agagCCTTACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGT‑CCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTc                                                                                                                                                 <  5:244180/149‑1 (MQ=255)        ccTTACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGT‑CCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGACCTCGAAg                                                                                                                                             >  8:207778/1‑149 (MQ=255)          ttACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGT‑CCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTc                                                                                                                                           >  5:166840/1‑149 (MQ=255)          ttACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGT‑CCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTc                                                                                                                                           >  1:76225/1‑149 (MQ=255)            aCAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGT‑CCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTc                                                                                                                                           <  4:377613/147‑1 (MQ=255)                aggtaaaggtGTGCGTTACGCGGACGAAGTAGT‑CCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACTACAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTgc                                                                                                                                       <  1:310468/149‑1 (MQ=255)                 ggtaaaggtGTGCGTTACGCGGACGAAGTAGTGCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTgcgc                                                                                                                                     >  1:243543/1‑149 (MQ=255)                 ggtaaaggtGTGCGTTACGCGGACGAAGTAGT‑CCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGACTTCGAGTTCGTgcatc                                                                                                                                    <  8:335049/149‑4 (MQ=255)                 ggtaaaggtGTGCGTTACGCGGACGAAGTAGT‑CCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGcc                                                                                                                                    >  2:63815/1‑149 (MQ=255)                   taaaggtGTGCGTTACGCGGACGAAGTAGT‑CCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGACTTCGAGTTCGTgcatctt                                                                                                                                  <  8:5556/149‑6 (MQ=255)                    aaaggtGTGCGTTACGCGGACGAAGTAGT‑ACGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGAACGAC‑A‑AAAAAGTTATTAGACTGCGTCGCGCTCGGAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCtgt                                                                                                                                 >  7:242128/1‑149 (MQ=255)                      aggtGTGCGGTACGCCGACGAAGTAGT‑CCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAACCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGc                                                                                                                               <  4:298292/149‑1 (MQ=255)                             cGTTACGCGGACGAAGTAGT‑CCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTc                                                                                                                         <  5:29949/148‑1 (MQ=255)                              gTTACGCGGACGAAGTAGT‑CCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCg                                                                                                                       >  8:218140/1‑149 (MQ=255)                              gTTACGCGGACGAAGTAGT‑CCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCg                                                                                                                       >  7:283241/1‑149 (MQ=255)                              gTTACGCGGACGAAGTAGT‑CCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCg                                                                                                                       >  6:362480/1‑149 (MQ=255)                              gTTACGCGGACGAAGTAGT‑CCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCg                                                                                                                       <  3:241458/149‑1 (MQ=255)                                  cgcgGACGAAGTAGT‑CCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGctcc                                                                                                                   <  1:340785/149‑1 (MQ=255)                                  cgcgGACGAAGTAGT‑CCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGctcc                                                                                                                   >  2:194749/1‑149 (MQ=255)                                      gACGAAGTAGT‑CCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCgc                                                                                                               <  2:113325/149‑1 (MQ=255)                                         gAAGTAGT‑CCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCgcag                                                                                                             >  5:261516/1‑148 (MQ=255)                                           agtagt‑CCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCgcagca                                                                                                           >  6:349427/1‑148 (MQ=255)                                                          aagaagCCAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAgg                                                                                             <  1:63815/148‑1 (MQ=255)                                                           agaagGCAAGATGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTTCGCCCAGGt                                                                                            <  8:335622/148‑1 (MQ=255)                                                             aagCCAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGGCGCAAGATGCACGCACTCGAAGTCGTGCGCCTGTGCGTGTGCCGATGCTCTCCGCAGATCTACGCCCAGGTCAt                                                                                         >  4:292642/1‑149 (MQ=255)                                                                ccAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGAGTCGCGCTCGCAAAGCCCGTCTCAAGATGCACGCACTCGAAGTCGTGCGCCTGTGCGTGTTCAGCTCCTCGCAGCACATCAACGCCGAGGTCATTTc                                                                                      >  2:244314/1‑149 (MQ=255)                                                                caaagaagaagTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTc                                                                                      <  3:193621/147‑1 (MQ=255)                                                                 cAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGc                                                                                                                               >  7:105329/1‑107 (MQ=255)                                                                 cAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGc                                                                                                                               <  8:105329/107‑1 (MQ=255)                                                                  aagaagaagTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACTAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCgg                                                                                    >  6:127149/1‑149 (MQ=255)                                                                  aagaagaagTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCgg                                                                                    >  2:284057/1‑149 (MQ=255)                                                                    gaagaagTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGcc                                                                                  >  3:303613/1‑149 (MQ=255)                                                                    gaagaagTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGcc                                                                                  <  6:261516/149‑1 (MQ=255)                                                                       gaagTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGAc                                                                               >  5:124105/1‑149 (MQ=255)                                                                       gaagTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGAc                                                                               <  5:29649/149‑1 (MQ=255)                                                                        aagTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACg                                                                              >  4:203179/1‑149 (MQ=255)                                                                          gTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACgg                                                                             >  4:35065/1‑148 (MQ=255)                                                                           tAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGgca                                                                           <  2:77137/149‑1 (MQ=255)                                                                                ccTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAg                                                                      <  7:207778/149‑1 (MQ=255)                                                                                       actGACCGAC‑A‑AAAACGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTGCt                                                                 <  8:242128/145‑1 (MQ=255)                                                                                       aaTGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTg                                                                <  2:243543/148‑1 (MQ=255)                                                                                           aCCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCaa                                                            <  5:333567/148‑1 (MQ=255)                                                                                           aCCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAg                                                           <  7:218140/149‑1 (MQ=255)                                                                                                  a‑aaaaaGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAGCTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTc                                                     >  4:227694/1‑149 (MQ=255)                                                                                                  a‑aaaaaGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTc                                                     >  5:186030/1‑149 (MQ=255)                                                                                                    aaaaaGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTc                                                     >  1:304069/1‑148 (MQ=255)                                                                                                     aaaaGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGa                                                   >  4:86648/1‑149 (MQ=255)                                                                                                       aaGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGAcc                                                 <  5:349427/149‑1 (MQ=255)                                                                                                          ttATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTg                                              <  2:201650/149‑1 (MQ=255)                                                                                                           tATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGACTTCGAGTTCGTgc                                                                                                                                       >  2:47982/1‑59 (MQ=37)                                                                                                           tATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGACTTCGAGTTCGTgc                                                                                                                                       <  1:47982/59‑1 (MQ=37)                                                                                                             ttCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACTAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGa                                            <  5:127149/148‑1 (MQ=255)                                                                                                             ttCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACTAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGAc                                           <  3:35065/149‑1 (MQ=255)                                                                                                                     cGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCt                                   <  6:1544/149‑1 (MQ=255)                                                                                                                        cgcgCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCg                                <  5:362480/149‑1 (MQ=255)                                                                                                                             tCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTg                                                                >  8:314749/1‑112 (MQ=255)                                                                                                                             tCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTg                                                                <  7:314749/112‑1 (MQ=255)                                                                                                                               gCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATg                           <  5:249534/147‑1 (MQ=255)                                                                                                                                 aaaGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGgcgc                       >  7:138055/1‑149 (MQ=255)                                                                                                                                 aaaGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGgcgc                       >  6:254287/1‑149 (MQ=255)                                                                                                                                 aaaGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGgcgc                       >  5:8615/1‑149 (MQ=255)                                                                                                                                 aaaGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGgcgc                       >  8:168750/1‑149 (MQ=255)                                                                                                                                      aCGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTg                  <  8:283241/149‑1 (MQ=255)                                                                                                                                      aCGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTg                  >  4:73971/1‑149 (MQ=255)                                                                                                                                        gTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTgg                 >  6:275125/1‑148 (MQ=255)                                                                                                                                        gTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGc                <  6:166840/149‑1 (MQ=255)                                                                                                                                        gTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGc                <  2:76225/149‑1 (MQ=255)                                                                                                                                           tcAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCaa              >  4:269836/1‑148 (MQ=255)                                                                                                                                             aaGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAg                                       >  6:394555/1‑121 (MQ=255)                                                                                                                                             aaGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAg                                       <  5:394555/121‑1 (MQ=255)                                                                                                                                                gcgCACGACTTCGAGTTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGa        >  1:233919/3‑149 (MQ=255)                                                                                                                                                aTGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGa        <  3:86648/149‑1 (MQ=255)                                                                                                                                                      gAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGACGCGGc  <  1:244314/149‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                   GCCAGAGCCTTACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGT‑CCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGACGCGGC  >  NC_002947/558329‑558623

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