Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	558,478	C→T	6.7%	L22F (CTC→TTC)	rplR →	50S ribosomal protein L18

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	558,478	0	C	T	6.7%	219.8 / 3.7	74	L22F (CTC→TTC)	rplR	50S ribosomal protein L18
Reads supporting (aligned to +/- strand):  ref base C (34/35);  new base T (2/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 = 9.44e-01

AGAGCCTTACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGT‑CCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGACGCGGC  >  NC_002947/558332‑558623                                                                                                                                                      |                                                                                                                                                  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                                                                                                                   >  2:194749/1‑149 (MQ=255)                               cgcgGACGAAGTAGT‑CCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGctcc                                                                                                                   <  1:340785/149‑1 (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                                                                                                                                       <  1:47982/59‑1 (MQ=37)                                                                                                        tATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGACTTCGAGTTCGTgc                                                                                                                                       >  2:47982/1‑59 (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                       >  6:254287/1‑149 (MQ=255)                                                                                                                              aaaGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGgcgc                       >  7:138055/1‑149 (MQ=255)                                                                                                                              aaaGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGgcgc                       >  5:8615/1‑149 (MQ=255)                                                                                                                              aaaGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGgcgc                       >  8:168750/1‑149 (MQ=255)                                                                                                                                   aCGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTg                  >  4:73971/1‑149 (MQ=255)                                                                                                                                   aCGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTg                  <  8:283241/149‑1 (MQ=255)                                                                                                                                     gTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTgg                 >  6:275125/1‑148 (MQ=255)                                                                                                                                     gTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGc                <  2:76225/149‑1 (MQ=255)                                                                                                                                     gTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGc                <  6:166840/149‑1 (MQ=255)                                                                                                                                        tcAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCaa              >  4:269836/1‑148 (MQ=255)                                                                                                                                          aaGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAg                                       <  5:394555/121‑1 (MQ=255)                                                                                                                                          aaGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAg                                       >  6:394555/1‑121 (MQ=255)                                                                                                                                             gcgCACGACTTCGAGTTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGa        >  1:233919/3‑149 (MQ=255)                                                                                                                                             aTGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGa        <  3:86648/149‑1 (MQ=255)                                                                                                                                                   gAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGACGCGGc  <  1:244314/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                  AGAGCCTTACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGT‑CCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGAC‑A‑AAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGACGCGGC  >  NC_002947/558332‑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: ‑