Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 501,419 A→T 6.2% K268M (AAG→ATG)  PP_0412 → polyamine ABC transportersubstrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947501,4190AT 6.2% 188.2 / 2.7 64K268M (AAG→ATG) PP_0412polyamine ABC transportersubstrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (27/33);  new base T (2/2);  total (29/35)
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.75e-01

TGGCGACGTGGTCATGAGCTCTGCCTACAACGGCCGTATTGCCGCTGTGCAGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGACATGCCGACCACGCCTGAGAACATTGC  >  NC_002947/501273‑501563
                                                                                                                                                              |                                                                                                                                                
tGGCGACGTTGTCATGAGCTCTGCCTACAACGGCCGTATTGCCGCTGTGCAGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑gaa                                                                                                                                                  <  2:44372/147‑1 (MQ=255)
tGGCGACGTGGTCATGAGCTCTGCCTACAACGGCCGTATTGCCGCTGTGCAGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑gaa                                                                                                                                                  <  8:90256/147‑1 (MQ=255)
tGGCGACGTGGTCATGAGCTCTGCCTACAACGGCCGTATTGCCGCTGTGCAGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑gaa                                                                                                                                                  <  2:50091/147‑1 (MQ=255)
  gCGACGTGGTCATGAGCTCTGCCTACAACGGCCGTATTGCCGCTGTGCAGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑gaagaa                                                                                                                                               <  3:139176/148‑1 (MQ=255)
  gCGACGTGGTCATGAGCTCTGCCTACAACGGCCGTATTGCCGCTGTGCAGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑gaa                                                                                                                                                  <  1:325241/145‑1 (MQ=255)
  gCGACGTGGTCATGAGCTCTGCCTACAACGGCCGTATTGCCGCTGTGCAGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑Gaagaag                                                                                                                                              >  8:328157/1‑149 (MQ=255)
   cGACGTGGTCATGAGCTCTGCCTACAACGGCCGTATTGCCGCTGTGCAGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGt                                                                                                                                             <  1:331093/149‑1 (MQ=255)
   cGACGTGGTCATGAGCTCTGCCTACAACGGCCGTATTGCCGCTGTGCAGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGt                                                                                                                                             >  3:154601/1‑149 (MQ=255)
           tCATGAGCTCTGCCTACAACGGCCGTATTGCCGCTGTGCAGAAAGAGAGCAACCTCTAGGTGGTGTGGCACGGCGGCATCTGCGACTGCGATGCATGGGCAATTCCGAAAGGCGCCACACACGTC‑‑‑‑GCAGCGGG‑‑‑‑‑‑GAAGGAGTTCATGGcc                                                                                                                                     >  6:311498/1‑149 (MQ=255)
             aTGAGCTCTGCCTACAACGGCCGTATTGCCGCTGTGCAGAAAGATAGCAACCTCAAGGTGTTTGGGAACGGAGGCATCTACGACTTCGATGCATTGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTa                                                                                                                                   <  6:361934/149‑1 (MQ=255)
             aTGAGCTCTGCCTACAACGGCCGTATTGCCGCTGTGCAGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTa                                                                                                                                   <  8:60284/149‑1 (MQ=255)
              tGAGCTCTGCCTACAACGGCCGTATTGCCGCTGTGCAGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTAc                                                                                                                                  <  6:322183/149‑1 (MQ=255)
               gAGCTCTGCCTACAACGGCCGTATTGCCGCTGTGCAGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTa                                                                                                                                   >  6:105794/1‑147 (MQ=255)
                              cGGCCGTATTGCCGCTGTGCAGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGGGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGc                                                                                                                               <  4:218825/136‑1 (MQ=255)
                              cGGCCGTATTGCCGCTGTGCAGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGGGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGc                                                                                                                               >  3:218825/1‑136 (MQ=255)
                                gCCGTATTGCAGCTGTGGAGTAAGAGAGCAACCTCAAGGTTGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGagcag                                                                                                                <  1:369871/149‑1 (MQ=255)
                                 ccGTATTGCCGCTGTGCAGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑Gaagaag                                                                                                                                              <  3:125056/118‑1 (MQ=255)
                                 ccGTATTGCCGCTGTGCAGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑Gaagaag                                                                                                                                              >  4:125056/1‑118 (MQ=255)
                                 ccGTATTGCCGCTGTGCAGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGagcagc                                                                                                               <  8:274569/149‑1 (MQ=255)
                                  cGTATTGCCGCTGTGCAGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAgcagca                                                                                                              <  7:328157/149‑1 (MQ=255)
                                      ttGCCGCTGTGCAGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGcagcag                                                                                                             >  7:345320/1‑146 (MQ=255)
                                        gCCGCTGTGCAGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGcagcag                                                                                                             >  2:15911/1‑144 (MQ=255)
                                         ccGCTGTGCAGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATc                                                                                                       <  5:142226/149‑1 (MQ=255)
                                          cGCTGTGCAGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATc                                                                                                       >  4:144241/1‑148 (MQ=255)
                                                  aGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGAt                                                                                                        <  2:24157/139‑1 (MQ=255)
                                                  aGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGAt                                                                                                        >  1:24157/1‑139 (MQ=255)
                                                  aGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGa                                                                                               >  8:119855/1‑148 (MQ=255)
                                                  aGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGa                                                                                               >  3:312265/1‑148 (MQ=255)
                                                      agagagCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACa                                                                                          <  7:118170/149‑1 (MQ=255)
                                                        agagCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATc                                                                                        >  6:143777/1‑149 (MQ=255)
                                                           gCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGc                                                                                      <  1:15911/148‑1 (MQ=255)
                                                            cAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCg                                                                                       >  6:58524/1‑146 (MQ=255)
                                                            cAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCg                                                                                       <  5:58524/146‑1 (MQ=255)
                                                                cTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTCGAAGGAAGAGGGCGGAGAGGAGAGAAAATCATCTCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCAGaaaa                                                                                            >  1:122062/1‑149 (MQ=255)
                                                                cTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATgg                                                                                <  5:311498/149‑1 (MQ=255)
                                                                   aaGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGccc                                                                             <  4:154601/149‑1 (MQ=255)
                                                                                      gcATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGt                                                           >  4:333335/1‑148 (MQ=255)
                                                                                          cTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCaa                                                                 <  1:295279/138‑1 (MQ=255)
                                                                                          cTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCaa                                                                 >  2:295279/1‑138 (MQ=255)
                                                                                             cGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCctgc                                                   >  3:154441/1‑149 (MQ=255)
                                                                                                cTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGt                                                <  5:143777/149‑1 (MQ=255)
                                                                                                 ttCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATc                                                                                        <  6:323321/108‑1 (MQ=255)
                                                                                                       gCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACg                                          >  4:218038/1‑148 (MQ=255)
                                                                                                            ggCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTg                                     <  2:122062/148‑1 (MQ=255)
                                                                                                              cAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTgaag                                  <  3:333335/149‑1 (MQ=255)
                                                                                                               aaTTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCaa                                                                       <  5:367217/111‑1 (MQ=255)
                                                                                                               aaTTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCaa                                                                       >  6:367217/1‑111 (MQ=255)
                                                                                                               aaTTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTgaag                                  >  6:6431/1‑148 (MQ=255)
                                                                                                                    ccAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGaagaaaga                             >  1:18068/1‑148 (MQ=255)
                                                                                                                        aGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGACATGc                        <  7:119855/149‑1 (MQ=255)
                                                                                                                            gcCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTa                                                                                   >  4:35027/1‑86 (MQ=255)
                                                                                                                            gcCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTa                                                                                   <  3:35027/86‑1 (MQ=255)
                                                                                                                            gcCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGACATGCCGAc                    >  8:70529/1‑149 (MQ=255)
                                                                                                                                aaaaCGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGACATGCCGACCACg                >  4:22891/1‑149 (MQ=255)
                                                                                                                                 aaaCGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGcc                                                             <  1:356383/103‑1 (MQ=255)
                                                                                                                                 aaaCGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGcc                                                             >  2:356383/1‑103 (MQ=255)
                                                                                                                                 aaaCGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGACATGCCGACCACGc               <  4:312265/149‑1 (MQ=255)
                                                                                                                                 aaaCGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGACATGCCGACCACGc               <  3:218038/149‑1 (MQ=255)
                                                                                                                                   aCGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGACATGCCGACCACGCCt             <  8:345320/149‑1 (MQ=255)
                                                                                                                                   aCGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGACATGCCGACCACGCCt             <  2:18068/149‑1 (MQ=255)
                                                                                                                                    cGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCaa                                                                 <  6:161251/96‑1 (MQ=255)
                                                                                                                                    cGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCaa                                                                 >  5:161251/1‑96 (MQ=255)
                                                                                                                                            gAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGACATGCCGACCACGCCTGAGaa        <  3:144241/149‑1 (MQ=255)
                                                                                                                                            gAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGCCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGACATGCCGACCACGCCTGAGaa        <  7:70529/149‑1 (MQ=255)
                                                                                                                                             tgtagGC‑‑‑‑‑‑GATGAACTTCTTCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGAc                            >  4:189889/4‑128 (MQ=255)
                                                                                                                                             tgtagGC‑‑‑‑‑‑GATGAACTTCTTCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGAc                            <  3:189887/125‑1 (MQ=255)
                                                                                                                                             tgtagGC‑‑‑‑‑‑GATGAACTTCTTCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGAc                            <  3:189889/125‑1 (MQ=255)
                                                                                                                                             tgtagGC‑‑‑‑‑‑GATGAACTTCTTCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGAc                            >  4:189887/4‑128 (MQ=255)
                                                                                                                                             aaGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTga                                    >  8:315324/1‑120 (MQ=255)
                                                                                                                                             aaGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTga                                    <  7:315324/120‑1 (MQ=255)
                                                                                                                                                  ggC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGACATGCCGACCACGCCTGAGAACATTg   <  5:105794/149‑1 (MQ=255)
                                                                                                                                                   gC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGTAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGACATGCCGACCACGCCTGAGAACATTGc  >  2:342410/1‑149 (MQ=255)
                                                                                                                                                              |                                                                                                                                                
TGGCGACGTGGTCATGAGCTCTGCCTACAACGGCCGTATTGCCGCTGTGCAGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGACATGCCGACCACGCCTGAGAACATTGC  >  NC_002947/501273‑501563

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: