Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 503,837:1 +C 10.4% intergenic (+40/‑111) PP_0414 → / → rpe polyamine ABC transporter permease/ribulose‑5‑phosphate 3‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947503,8371.C10.4% 286.1 / 28.7 77intergenic (+40/‑111)PP_0414/rpepolyamine ABC transporter permease/ribulose‑5‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base . (32/37);  new base C (2/6);  total (34/43)
Fisher's exact test for biased strand distribution p-value = 2.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.03e-01

CATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCCATTCTCTCTGC  >  NC_002947/503689‑503985
                                                                                                                                                      |                                                                                                                                                     
cATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCg                                                                                                                                                        <  3:176255/149‑1 (MQ=255)
  ttGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCCGGCTTTAGCCgcg                                                                                                                                                      >  6:74417/1‑146 (MQ=255)
    gCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTTCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGaa                                                                                                                                                    <  3:209249/148‑1 (MQ=255)
    gCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCCGGCTTTAGCCgcgaa                                                                                                                                                    <  6:217945/148‑5 (MQ=255)
     ccGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTACCTGTGCCGGCCTCTTCCGGCTTTAGCCGCGAaga                                                                                                                                                 <  5:23897/149‑1 (MQ=38)
      cgcgGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTACCTGTGCCGGCCTCTTCCGGCTTTAGCCGCGAagag                                                                                                                                                <  8:325381/149‑1 (MQ=38)
      cgcgGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTACCTGTGCCGGCCTCTTCCGGCTTTAGCCGCGAagag                                                                                                                                                <  6:30594/149‑1 (MQ=38)
        cggcggcgACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCTGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGc                                                                                                                                              <  7:50066/149‑1 (MQ=255)
        cggcggcgACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCCGGCTTTAGCCGCGAAGAGGc                                                                                                                                              <  6:173165/149‑1 (MQ=39)
           cggcgACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGTCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCgg                                                                                                                                           <  2:250468/149‑1 (MQ=255)
           cggcgACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCgg                                                                                                                                           <  2:241917/149‑1 (MQ=255)
           cggcgACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCgg                                                                                                                                           <  2:229035/149‑1 (MQ=255)
               gACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGca                                                                                                                                         >  6:250371/1‑147 (MQ=255)
                 ccTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCgg                                                                                                                                                      <  2:136181/133‑1 (MQ=255)
                 ccTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGCCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAgg                                                                                                                                     >  1:134700/1‑149 (MQ=255)
                  cTTGCTGATTGCGTTCTCGGTTTTTTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGTCTTGATGATGGCTTGCCTGTGCCGGCCTCTTCCGGCTTTAGCCGCGAAGAGGCCGGTACAggc                                                                                                                                    <  5:378768/149‑2 (MQ=25)
                  cTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGt                                                                                                                                    <  4:222758/149‑1 (MQ=255)
                   ttGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGtt                                                                                                                                   >  5:252947/1‑149 (MQ=255)
                   ttGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGtt                                                                                                                                   >  2:38022/1‑149 (MQ=255)
                     tcTGATTGCGTTCTCTGGGTTGTTGCTGATGACCCTGTTGGGGTGGCGTGGGCTTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATTATGGCTTGCCTGTACCGGCCTCTTCGCGGCAAAAGCCGGAAGATGCCGGCACAGGtt                                                                                                                                   <  6:47808/146‑1 (MQ=255)
                          ttGGGTTCTCGGTGTTGGTGCTGCTGACCCTGGAGTGGTGGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTg                                                                                                                             <  8:278820/148‑1 (MQ=255)
                          ttGCTTTCTGGGTGTTGTTGCTGCTGACCGTGGAGTGGGTGCGTGGGCGTAGCGAGTAGCTGATGTGTCAACAGCCGGTTTGATGATGGCTTGCCTGTACGGGTCTCTGCGCGGTTAAAGCGGGAAGAGGCCGGCACAGTTTGCGTTg                                                                                                                             <  2:134700/148‑1 (MQ=255)
                             cGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTc                                                                                                                          <  3:190998/148‑1 (MQ=255)
                             cGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTc                                                                                                                          <  7:162848/148‑1 (MQ=255)
                              gTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGa                                                                                                                        >  2:60621/1‑149 (MQ=255)
                               ttCTCGTTGGTGTTGCTGCTGACCCTGGCGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGac                                                                                                                       <  7:141407/149‑1 (MQ=255)
                               ttCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGCAGGCCGGCACAGGTTGCGTTGGTCGa                                                                                                                        <  1:4878/149‑1 (MQ=255)
                                    ggTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCaa                                                                                                                  >  7:305766/1‑149 (MQ=255)
                                     gtgtGGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACCGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAAt                                                                                                                 <  4:8303/149‑1 (MQ=255)
                                      tgttgtTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAAtt                                                                                                                <  1:301661/149‑1 (MQ=255)
                                          gttgCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGAt                                                                                                             >  2:353720/1‑148 (MQ=255)
                                           ttgCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATa                                                                                                            <  6:233363/148‑1 (MQ=255)
                                           ttgCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATa                                                                                                            <  8:305766/148‑1 (MQ=255)
                                                 ctgACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATc                                                                                                      <  5:250371/148‑1 (MQ=255)
                                                    aCCCTGGAGTGGTTGCGTGGGCGTAGCGTGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGcc                                                                                                  <  8:257272/149‑1 (MQ=255)
                                                     cccTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGcc                                                                                                  <  1:60621/148‑1 (MQ=255)
                                                      ccTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCtt                                                                                                <  4:348592/149‑1 (MQ=255)
                                                       cTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTGCCGGCCTCTTCCGGCTTTAGCCGCGAaga                                                                                                                                                 <  2:245391/99‑1 (MQ=25)
                                                       cTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTGCCGGCCTCTTCCGGCTTTAGCCGCGAaga                                                                                                                                                 >  1:245391/1‑99 (MQ=25)
                                                        tGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCa                                                                                              >  1:145396/1‑149 (MQ=255)
                                                          gAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCAt                                                                                             <  1:353720/148‑1 (MQ=255)
                                                           aGTGGTTGCGTGGGCGTAGTGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCg                                                                                           <  1:239260/149‑1 (MQ=255)
                                                           aGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCg                                                                                           <  7:7315/149‑1 (MQ=255)
                                                           aGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCg                                                                                           >  6:69992/1‑149 (MQ=255)
                                                                    gTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCt                                                                                   >  3:171758/1‑148 (MQ=255)
                                                                     tGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCg                                                                                           <  1:62092/139‑1 (MQ=255)
                                                                     tGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCg                                                                                           >  2:62092/1‑139 (MQ=255)
                                                                     tGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGa                                                                                 <  1:38022/149‑1 (MQ=255)
                                                                          gTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTa                                                                             >  1:345988/1‑148 (MQ=255)
                                                                          gTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTa                                                                             >  8:263086/1‑148 (MQ=255)
                                                                          gTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTAc                                                                            <  5:69992/149‑1 (MQ=255)
                                                                          gTAGCGAGAAGCTGAGGACTCAACAGACGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTTATACCAATCAGCCTTCATCGCTTTCCCTGAGTTa                                                                             >  6:82080/1‑148 (MQ=255)
                                                                              cGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATg                                                                        >  1:24153/1‑149 (MQ=255)
                                                                               gagaAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATgc                                                                       >  5:281168/1‑149 (MQ=255)
                                                                               gagaAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATgc                                                                       >  1:51107/1‑149 (MQ=255)
                                                                                agaAGCTGAGGACTCATCTGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACTAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATgcg                                                                      <  1:284018/149‑1 (MQ=255)
                                                                                agaAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACGCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATgc                                                                       >  1:85932/1‑149 (MQ=255)
                                                                                   aGCTGATGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATgcgcgc                                                                   <  6:39070/149‑1 (MQ=255)
                                                                                         ggACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGt                                                             <  7:263086/149‑1 (MQ=255)
                                                                                             tCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAAt                                                         >  7:50336/1‑149 (MQ=255)
                                                                                             tCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAAt                                                         >  4:117288/1‑149 (MQ=255)
                                                                                               aaCAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATgcgcgc                                                                   >  7:10851/1‑137 (MQ=255)
                                                                                               aaCAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATgcgcgc                                                                   <  8:10851/137‑1 (MQ=255)
                                                                                               aaCAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAAttt                                                       >  5:1994/1‑149 (MQ=255)
                                                                                               aaCAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAAttt                                                       <  2:345988/149‑1 (MQ=255)
                                                                                                            atgatgGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCACTGAGTTACAATGCGAGCCACCCTGAATTTGCCAACAGGTgcg                                          >  2:103600/1‑149 (MQ=255)
                                                                                                              gatgGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTGACAATGCGCGCCACCGTGAATTTGCCAACCGGTgcgc                                         >  6:310753/1‑148 (MQ=255)
                                                                                                                 ggCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACAGCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATg                                     >  1:259964/1‑149 (MQ=255)
                                                                                                                   cTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCa                                   <  2:85932/149‑1 (MQ=255)
                                                                                                                      gCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGcc                                >  4:86237/1‑149 (MQ=255)
                                                                                                                      gCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGcc                                >  4:322667/1‑149 (MQ=255)
                                                                                                                        cTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGccc                               >  3:160803/1‑148 (MQ=255)
                                                                                                                        cTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCt                              <  4:171758/149‑1 (MQ=255)
                                                                                                                                    tctTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGcc                   >  1:103667/1‑149 (MQ=255)
                                                                                                                                        cgcgGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGccc                  <  2:24153/146‑1 (MQ=255)
                                                                                                                                         gcgGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTc              >  5:380793/1‑149 (MQ=255)
                                                                                                                                         ccgGCTTTAGCCGCGAAGAGGCCGGTACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCt               >  3:275329/2‑149 (MQ=255)
                                                                                                                                            gCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCCAt           <  2:145396/149‑1 (MQ=255)
                                                                                                                                                  gCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCt               <  8:256513/139‑1 (MQ=255)
                                                                                                                                                  gCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCt               >  7:256513/1‑139 (MQ=255)
                                                                                                                                                   ccGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCCATTctctct    >  5:288346/1‑149 (MQ=255)
                                                                                                                                                     ggAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCCATTCTCTCTGc  >  7:302524/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCCATTCTCTCTGC  >  NC_002947/503689‑503985

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: