Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 311,012 T→C S109G (AGC→GGC)  ecpR ← putative transcriptional regulator for the ecp operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913311,0120TC99.0% 396.9 / NA 100S109G (AGC→GGC) ecpRputative transcriptional regulator for the ecp operon
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (51/48);  minor base A (1/0);  total (52/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGGCCGTCATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCTATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGTTTATAAATTCAGATAAAGAATCCGTACTGAGATTGATAAACACGAAAT  >  NC_000913/310864‑311158
                                                                                                                                                    |                                                                                                                                                  
tGGGCCGTCATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCAt                                                                                                                                                  >  2:1269462/1‑151 (MQ=255)
 gggCCGTCATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCAtt                                                                                                                                                 <  4:1483933/151‑1 (MQ=255)
 gggCCGTCATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCAt                                                                                                                                                  <  2:1559005/150‑1 (MQ=255)
 gggCCGTCATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCAt                                                                                                                                                  <  4:1233974/150‑1 (MQ=255)
  ggCCGTCATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTg                                                                                                                                                <  4:438967/151‑1 (MQ=255)
   gCCGTCATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGt                                                                                                                                               <  2:1183093/151‑1 (MQ=255)
      gTCATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAAt                                                                                                                                            <  4:1039343/151‑1 (MQ=255)
        cATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATa                                                                                                                                           <  2:927004/150‑1 (MQ=255)
         aTGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATAcc                                                                                                                                         >  2:1372024/1‑151 (MQ=255)
          tGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATAcc                                                                                                                                         >  4:201330/1‑150 (MQ=255)
           gCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCaa                                                                                                                                       <  4:1679812/151‑1 (MQ=255)
           gCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCa                                                                                                                                        >  2:1504090/1‑150 (MQ=255)
           gCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCa                                                                                                                                        <  2:930424/150‑1 (MQ=255)
            cGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCGGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCaat                                                                                                                                      <  3:400056/151‑1 (MQ=255)
            cGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCaat                                                                                                                                      >  3:21441/1‑151 (MQ=255)
               aTAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCaataat                                                                                                                                   >  2:414460/1‑151 (MQ=255)
                tAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAAtt                                                                                                                                  >  4:652150/1‑151 (MQ=255)
                  atttccatttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTg                                                                                                                                <  1:478141/151‑1 (MQ=255)
                    ttccatttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCt                                                                                                                              >  3:467319/1‑151 (MQ=255)
                     tccatttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCt                                                                                                                              <  4:904604/150‑1 (MQ=255)
                     tccatttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCt                                                                                                                              <  2:1123974/150‑1 (MQ=255)
                     tccatttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCt                                                                                                                              <  2:1592546/150‑1 (MQ=255)
                      ccatttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTaa                                                                                                                            <  3:319353/151‑1 (MQ=255)
                      ccatttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTa                                                                                                                             <  2:734561/150‑1 (MQ=255)
                       catttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTa                                                                                                                             >  4:1438466/1‑149 (MQ=255)
                        atttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGc                                                                                                                         >  3:14208/1‑151 (MQ=255)
                         tttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAgg                                                                                                                          <  2:440765/150‑1 (MQ=255)
                          ttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGc                                                                                                                         <  4:1131418/150‑1 (MQ=255)
                          ttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTAAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGcc                                                                                                                        >  4:1157987/1‑151 (MQ=255)
                            cccGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCctc                                                                                                                      >  4:773537/1‑151 (MQ=255)
                             ccGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCctct                                                                                                                     >  4:56346/1‑151 (MQ=255)
                              cGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTg                                                                                                                    <  4:1659968/151‑1 (MQ=255)
                                 tCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACt                                                                                                                 >  2:169088/1‑151 (MQ=255)
                                    gTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTc                                                                                                              >  4:1153221/1‑151 (MQ=255)
                                      gATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTg                                                                                                            >  4:1291399/1‑151 (MQ=255)
                                          ttATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGcc                                                                                                         <  4:1510630/150‑1 (MQ=255)
                                          ttATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCg                                                                                                        >  4:559497/1‑151 (MQ=255)
                                           tATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCg                                                                                                        <  3:9380/150‑1 (MQ=255)
                                            aTCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGc                                                                                                       >  4:838504/1‑150 (MQ=255)
                                              ctttcttAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGcaat                                                                                                    >  4:923438/1‑151 (MQ=255)
                                               tttcttAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGcaatc                                                                                                   >  4:586813/1‑151 (MQ=255)
                                                  cttAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGcaatca                                                                                                  <  2:1062592/149‑1 (MQ=255)
                                                    tAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGcaatcaat                                                                                                <  2:1447649/149‑1 (MQ=255)
                                                    tAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATAc                                                                                              >  2:637931/1‑151 (MQ=255)
                                                     aaTATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACg                                                                                             <  2:206980/151‑1 (MQ=255)
                                                      atatCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACg                                                                                             >  4:730122/1‑150 (MQ=255)
                                                          cTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTg                                                                                        >  3:305477/1‑151 (MQ=255)
                                                              tctcAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTc                                                                                    >  2:1031279/1‑151 (MQ=255)
                                                               ctcAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTcc                                                                                   >  2:1572534/1‑151 (MQ=255)
                                                                 cAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTcccc                                                                                 <  4:1453238/151‑1 (MQ=255)
                                                                 cAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTcccc                                                                                 >  2:1456053/1‑151 (MQ=255)
                                                                 cAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTcccc                                                                                 >  3:262539/1‑151 (MQ=255)
                                                                  aGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCt                                                                                >  4:1429475/1‑151 (MQ=255)
                                                                  aGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCt                                                                                >  2:1446487/1‑151 (MQ=255)
                                                                     aaCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTa                                                                             <  2:1310253/151‑1 (MQ=255)
                                                                          gCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACcgc                                                                         <  1:449543/150‑1 (MQ=255)
                                                                           ccATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATTTTGCAATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCgg                                                                       >  1:57153/1‑151 (MQ=255)
                                                                           ccATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCgg                                                                       <  4:640779/151‑1 (MQ=255)
                                                                            cATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACTCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCgg                                                                       >  4:1124567/1‑150 (MQ=255)
                                                                             aTTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGAt                                                                     <  1:196874/151‑1 (MQ=255)
                                                                             aTTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGAt                                                                     >  2:555743/1‑151 (MQ=255)
                                                                             aTTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAAGCAATACGACCTGTTTCCCCTTTACCGCGGAt                                                                     <  2:938917/151‑1 (MQ=255)
                                                                               tAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATa                                                                    >  4:86444/1‑150 (MQ=255)
                                                                                aaTCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATaa                                                                   <  2:288226/150‑1 (MQ=255)
                                                                                aaTCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAAc                                                                  >  4:954784/1‑151 (MQ=255)
                                                                                aaTCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAAc                                                                  <  2:530027/151‑1 (MQ=255)
                                                                                   cacaTAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACcat                                                               >  2:667957/1‑151 (MQ=255)
                                                                                    acaTAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACcatt                                                              <  2:493708/151‑1 (MQ=255)
                                                                                            cAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATc                                                       >  1:39614/1‑150 (MQ=255)
                                                                                              gTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACg                                                    >  2:772672/1‑151 (MQ=255)
                                                                                                tCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGt                                                  >  2:527213/1‑151 (MQ=255)
                                                                                                   tttCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGt                                                  <  1:438426/148‑1 (MQ=255)
                                                                                                     tCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGTTtat                                              <  4:415055/150‑1 (MQ=255)
                                                                                                      cTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGTTtata                                             >  4:386865/1‑150 (MQ=255)
                                                                                                       tAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGTTTATaa                                            >  2:1584998/1‑150 (MQ=255)
                                                                                                        aaTATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGTTTATaa                                            <  3:201821/149‑1 (MQ=255)
                                                                                                                gACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGTTTATAAATTCAGATaa                                  <  2:1406898/151‑1 (MQ=255)
                                                                                                                gACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGTTTATAAATTCAGATa                                   >  2:1680913/1‑150 (MQ=255)
                                                                                                                gACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGTTTATAAATTCAGATa                                   >  2:1612452/1‑150 (MQ=255)
                                                                                                                gACTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGTTTATAAATTCAGATa                                   >  2:1239306/1‑150 (MQ=255)
                                                                                                                 aCTCAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGTTTATAAATTCAGATaa                                  >  1:199593/1‑150 (MQ=255)
                                                                                                                    cAGTCCAGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGTTTATAAATTCAGATAAAGa                               <  1:82402/150‑1 (MQ=255)
                                                                                                                          aGCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGTTTATAAATTCAGATAAAGAATCCGTa                        <  2:411550/151‑1 (MQ=255)
                                                                                                                           gCGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGTTTATAAATTCAGATAAAGAATCCGTa                        >  4:669183/1‑150 (MQ=255)
                                                                                                                            cGTATACCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGTTTATAAATTCAGATAAAGAATCCGTACt                      >  1:476471/1‑151 (MQ=255)
                                                                                                                              tataCCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGTTTATAAATTCAGATAAAGAATCCGTACTg                     >  2:621975/1‑150 (MQ=255)
                                                                                                                              tataCCACGCCCCTAATATTGCCATTGGAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGTTTATAAATTCAGATAAAGAATCCGTACTg                     >  4:1618083/1‑150 (MQ=255)
                                                                                                                               ataCCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGTTTATAAATTCAGATAAAGAATCCGTACTga                    <  2:1547108/150‑1 (MQ=255)
                                                                                                                               ataCCACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGTTTATAAATTCAGATAAAGAATCCGTACTg                     >  4:1399154/1‑149 (MQ=255)
                                                                                                                                   cACGCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGTTTATAAATTCAGATAAAGAATCCGTACTGAGATTg               <  2:520391/151‑1 (MQ=255)
                                                                                                                                      gCCCCTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGTTTATAAATTCAGATAAAGAATCCGTACTGAGATTGATa            <  2:1425434/151‑1 (MQ=255)
                                                                                                                                        cccTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGTTTATAAATTCAGATAAAGAATCCGTACTGAGATTGATaaa          <  4:508246/151‑1 (MQ=255)
                                                                                                                                         ccTAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGCCTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGTTTATAAATTCAGATAAAGAATCCGTACTGAGATTGATa            <  4:1271058/148‑1 (MQ=255)
                                                                                                                                           tAATATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGTTTATAAATTCAGATAAAGAATCCGTACTGAGATTGATAAacac       >  3:474776/1‑151 (MQ=255)
                                                                                                                                            aaTATTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGTTTATAAATTCAGATAAAGAATCCGTACTGAGATTGATAAaca        <  2:1542403/149‑1 (MQ=255)
                                                                                                                                             atatTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGTTTATAAATTCAGATAAAGAATCCGTACTGAGATTGATAAACACg      >  2:1046808/1‑150 (MQ=255)
                                                                                                                                               atTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGTTTATAAATTCAGATAAAGAATCCGTACTGAGATTGATAAACACGaaa   <  4:67932/151‑1 (MQ=255)
                                                                                                                                               atTGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGCTTCCCCTTTACCGCGGATAACCATTCATTATCACGGTTTATAAATTCAGATAAAGAATCCGTACTGAGATTGATAAACACGaaa   <  4:184689/151‑1 (MQ=255)
                                                                                                                                                ttGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGTTTATAAATTCAGATAAAGAATCCGTACTGAGATTGATAAACACGAAAt  >  2:1238802/1‑151 (MQ=255)
                                                                                                                                                 tGCCATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGTTTATAAATTCAGATAAAGAATCCGTACTGAGATTGATAAACACGAAAt  <  4:1505966/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGGGCCGTCATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACCACGCCCCTAATATTGCTATTGTAATACCAATAATTTGCTAAGGCCTCTGACTTTCTGGCCGCAATCAATACGACCTGTTTCCCCTTTACCGCGGATAACCATTCATTATCACGGTTTATAAATTCAGATAAAGAATCCGTACTGAGATTGATAAACACGAAAT  >  NC_000913/310864‑311158

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: