Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 501,427 A→T 6.1% I271F (ATC→TTC)  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,4270AT 6.1% 199.7 / 2.8 66I271F (ATC→TTC) PP_0412polyamine ABC transportersubstrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (27/35);  new base T (2/2);  total (29/37)
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.06e-01

TCATGAGCTCTGCCTACAACGGCCGTATTGCCGCTGTGCAGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGACATGCCGACCACGCCTGAGAACATTGCCAACCAGGT  >  NC_002947/501284‑501572
                                                                                                                                                           |                                                                                                                                                 
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                                                                                                                                        >  3:218825/1‑136 (MQ=255)
                   cGGCCGTATTGCCGCTGTGCAGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGGGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGc                                                                                                                                        <  4:218825/136‑1 (MQ=255)
                     gCCGTATTGCAGCTGTGGAGTAAGAGAGCAACCTCAAGGTTGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGagcag                                                                                                                         <  1:369871/149‑1 (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                                                                                                                 >  1:24157/1‑139 (MQ=255)
                                       aGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGAt                                                                                                                 <  2:24157/139‑1 (MQ=255)
                                       aGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGa                                                                                                        >  3:312265/1‑148 (MQ=255)
                                       aGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGa                                                                                                        >  8:119855/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                                                                                                <  5:58524/146‑1 (MQ=255)
                                                 cAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCg                                                                                                >  6:58524/1‑146 (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                                                                          >  2:295279/1‑138 (MQ=255)
                                                                               cTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCaa                                                                          <  1:295279/138‑1 (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                                                                                >  6:367217/1‑111 (MQ=255)
                                                                                                    aaTTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCaa                                                                                <  5:367217/111‑1 (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                      <  2:18068/149‑1 (MQ=255)
                                                                                                                        aCGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGACATGCCGACCACGCCt                      <  8:345320/149‑1 (MQ=255)
                                                                                                                         cGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCaa                                                                          >  5:161251/1‑96 (MQ=255)
                                                                                                                         cGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCaa                                                                          <  6:161251/96‑1 (MQ=255)
                                                                                                                                 gAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGACATGCCGACCACGCCTGAGaa                 <  3:144241/149‑1 (MQ=255)
                                                                                                                                 gAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGCCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGACATGCCGACCACGCCTGAGaa                 <  7:70529/149‑1 (MQ=255)
                                                                                                                                  tgtagGC‑‑‑‑‑‑GATGAACTTCTTCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGAc                                     <  3:189889/125‑1 (MQ=255)
                                                                                                                                  tgtagGC‑‑‑‑‑‑GATGAACTTCTTCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGAc                                     <  3:189887/125‑1 (MQ=255)
                                                                                                                                  tgtagGC‑‑‑‑‑‑GATGAACTTCTTCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGAc                                     >  4:189887/4‑128 (MQ=255)
                                                                                                                                  tgtagGC‑‑‑‑‑‑GATGAACTTCTTCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGAc                                     >  4:189889/4‑128 (MQ=255)
                                                                                                                                  aaGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTga                                             <  7:315324/120‑1 (MQ=255)
                                                                                                                                  aaGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTga                                             >  8:315324/1‑120 (MQ=255)
                                                                                                                                       ggC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGACATGCCGACCACGCCTGAGAACATTg            <  5:105794/149‑1 (MQ=255)
                                                                                                                                        gC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGTAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGACATGCCGACCACGCCTGAGAACATTGc           >  2:342410/1‑149 (MQ=255)
                                                                                                                                                     aagTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGcc                                                                      <  6:85238/84‑1 (MQ=255)
                                                                                                                                                        ttCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGTAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGACATGCCGACCACGCCTGAGAACATTGCCAACCAGGt  <  1:342410/149‑1 (MQ=255)
                                                                                                                                                        ttCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGACATGc                                 <  1:348881/118‑1 (MQ=255)
                                                                                                                                                        ttCATCGCCAACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGACATGc                                 >  2:348881/1‑118 (MQ=255)
                                                                                                                                                           |                                                                                                                                                 
TCATGAGCTCTGCCTACAACGGCCGTATTGCCGCTGTGCAGAAAGAGAGCAACCTCAAGGTGGTGTGGAACGGCGGCATCTACGACTTCGATGCATGGGCAATTCCCAAAGGCGCCAAAAACGTC‑‑‑‑GAAGAGGC‑‑‑‑‑‑GAAGAAGTTCATCGCCTACAGCGTCAAGCCCGAGCAGCAGAAGATCTACTCCGAAAACATCGCCTATGGCCCGGCCAACTCCAAGGCCGTGAGCCTGCTGTCGGACGAAGTGAAGAAAGACATGCCGACCACGCCTGAGAACATTGCCAACCAGGT  >  NC_002947/501284‑501572

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: