Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,350,111 C→T intergenic (‑4/+140) yciW ← / ← fabI putative oxidoreductase/enoyl‑[acyl‑carrier‑protein] reductase, NADH‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,350,1110CT98.9% 351.5 / ‑4.1 93intergenic (‑4/+140)yciW/fabIputative oxidoreductase/enoyl‑[acyl‑carrier‑protein] reductase, NADH‑dependent
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (39/53);  total (40/53)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

CTTTTTCAAGGATCACGTCGAGTAGAAAGGGATCGCTGACCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGCGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGTATAACAGAGATAACGGGCGGCAGAACGCCGCCCATCTTTACCAACAGAACGATTATTT  >  NC_000913/1349966‑1350256
                                                                                                                                                 |                                                                                                                                                 
cTTTTTCAAGTATCACGTCGAGTAGAAAGGGATCGCTGCCCTTTGCGGCTACCTGGACCAGAGGCACAACGTCATACTCCGCTGTCCTGGCTTGCGTTTCATGATACCAAGGGCTTTTTCCGGGGACGTGGCGTTGTTCCTTGGGTGCTcc                                                                                                                                              >  1:1022725/1‑151 (MQ=255)
 tttttCAAGGATCACGTCGAGTAGAAAGGGATCGCTGACCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCt                                                                                                                                             >  2:789227/1‑151 (MQ=255)
  ttttCAAGGATCACGTCGAGTAGAAAGGGATCGCTGACCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCt                                                                                                                                             <  1:4595/150‑1 (MQ=255)
    ttCAAGGATCACGTCGAGTAGAAAGGGATCGCTGACCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTgg                                                                                                                                          >  2:39275/1‑151 (MQ=255)
    ttCAAGGATCACGTCGAGTAGAAAGGGATCGCTGACCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTg                                                                                                                                           <  1:818011/150‑1 (MQ=255)
     tCAAGGATCACGTCGAGTAGAAAGGGATCGCTGACCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGt                                                                                                                                         >  1:94288/1‑151 (MQ=255)
      cAAGGATCACGTCGAGTAGAAAGGGATCGCTGACCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTc                                                                                                                                        <  1:958182/151‑1 (MQ=255)
       aaGGATCACGTCGAGTAGAAAGGGATCGCTGACCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCg                                                                                                                                       >  2:718835/1‑151 (MQ=255)
        aGGATCACGTCGAGTAGAAAGGGATCGCTGACCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGt                                                                                                                                      <  1:421377/151‑1 (MQ=255)
          gATCACGTCGAGTAGAAAGGGATCGCTGACCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTaa                                                                                                                                    <  1:886783/151‑1 (MQ=255)
             cACGTCGAGTAGAAAGGGATCGCTGACCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAgg                                                                                                                                 >  2:694936/1‑151 (MQ=255)
             cACGTCGAGTAGAAAGGGATCGATGACCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAgg                                                                                                                                 <  1:223824/151‑1 (MQ=255)
                gTCGAGTAGAAAGGGATCGCTGACCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGaaa                                                                                                                              <  2:169805/151‑1 (MQ=255)
                gTCGAGTAGAAAGGGATCGCTGACCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGaa                                                                                                                               >  2:91044/1‑150 (MQ=255)
                  cGAGTAGAAAGGGATCGCTGACCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAAt                                                                                                                             <  2:1107821/150‑1 (MQ=255)
                  cGAGTAGAAAGGGATCGCTGACCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATc                                                                                                                            >  2:518593/1‑151 (MQ=255)
                  cGAGTAGAAAGGGATCGCTGACCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATc                                                                                                                            >  1:104229/1‑151 (MQ=255)
                      tAGAAAGGGATCGCTGACCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGtt                                                                                                                         >  2:364816/1‑150 (MQ=255)
                       aGAAAGGGATCGCTGACCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTAt                                                                                                                       <  1:360873/151‑1 (MQ=255)
                        gAAAGGGATCGCTGACCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTTCTGTATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATc                                                                                                                      <  1:1022509/151‑1 (MQ=255)
                        gAAAGGGATCGCTGACCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATc                                                                                                                      >  2:1206996/1‑151 (MQ=255)
                        gAAAGGGATCGCTGACCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATc                                                                                                                      >  2:744786/1‑151 (MQ=255)
                          aaGGGATCGCTGACCTTTGCGGCTTCTGGGACCAGAGGCAGAACGTCATTCTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCt                                                                                                                    <  2:669963/151‑1 (MQ=255)
                            gggATCGCTGACCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGa                                                                                                                  <  2:1075723/151‑1 (MQ=255)
                                tCGCTGACCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGc                                                                                                               <  2:826532/150‑1 (MQ=255)
                                 cGCTGACCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCa                                                                                                              >  2:535885/1‑150 (MQ=255)
                                  gCTGACCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGCGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCa                                                                                                              >  2:247106/1‑149 (MQ=255)
                                   cTGACCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAgg                                                                                                           >  1:638603/1‑151 (MQ=255)
                                     gACCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAgg                                                                                                           <  1:310355/149‑1 (MQ=255)
                                     gACCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGc                                                                                                          <  2:634049/150‑1 (MQ=255)
                                      aCCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCgg                                                                                                        >  2:605102/1‑151 (MQ=255)
                                           tGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAgg                                                                                                   >  1:161794/1‑151 (MQ=255)
                                            gCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAggg                                                                                                  <  2:411880/151‑1 (MQ=255)
                                             cGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAgggg                                                                                                 >  1:78494/1‑151 (MQ=255)
                                               gCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAgggg                                                                                                 >  1:760716/1‑149 (MQ=255)
                                               gCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGgag                                                                                               <  2:1041207/151‑1 (MQ=255)
                                                  tCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGaa                                                                                             >  2:614881/1‑150 (MQ=255)
                                                    cGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGa                                                                                          <  2:224734/151‑1 (MQ=255)
                                                     gggACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGAc                                                                                         >  2:724491/1‑151 (MQ=255)
                                                        aCCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGa                                                                                      <  1:787395/151‑1 (MQ=255)
                                                        aCCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGa                                                                                      <  1:955374/151‑1 (MQ=255)
                                                         ccAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGAt                                                                                     >  1:908893/1‑151 (MQ=255)
                                                         ccAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGAt                                                                                     <  2:172781/151‑1 (MQ=255)
                                                          cAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATc                                                                                    >  2:426219/1‑151 (MQ=255)
                                                          cAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATc                                                                                    <  1:694900/151‑1 (MQ=255)
                                                           agagGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCg                                                                                   >  1:690657/1‑151 (MQ=255)
                                                            gagGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCgg                                                                                  <  1:9236/151‑1 (MQ=255)
                                                                cAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGa                                                                               <  1:877007/150‑1 (MQ=255)
                                                                cAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGAt                                                                              <  2:55339/151‑1 (MQ=255)
                                                                cAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGAt                                                                              >  1:636382/1‑151 (MQ=255)
                                                                      gTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTGTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGTCGGATCGGCGATAACaaa                                                                        <  1:1047608/151‑1 (MQ=255)
                                                                           aCTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATc                                                                    <  1:672446/150‑1 (MQ=255)
                                                                           aCTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCa                                                                   <  2:1165687/151‑1 (MQ=255)
                                                                            cTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATc                                                                    >  2:327380/1‑149 (MQ=255)
                                                                             tCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATc                                                                    >  1:467823/1‑148 (MQ=255)
                                                                              ccGTAGTGCTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAg                                                                  <  1:97217/149‑1 (MQ=255)
                                                                              ccGTAGTACTGGATTGCGTTTCATGATTCCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGaa                                                                <  2:1008378/151‑1 (MQ=255)
                                                                              ccGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAg                                                                  >  2:571756/1‑149 (MQ=255)
                                                                               cGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGaaa                                                               <  1:723312/151‑1 (MQ=255)
                                                                                gtagtaCTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAg                                                              <  2:360459/151‑1 (MQ=255)
                                                                                 tagtaCTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAgg                                                             >  2:1035636/1‑151 (MQ=255)
                                                                                   gtaCTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGta                                                           <  2:81343/151‑1 (MQ=255)
                                                                                    taCTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGGGAAAGACGGCTCGGCGATAACAAATATCAGcaa                                                               >  1:165879/1‑146 (MQ=255)
                                                                                     aCTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGtat                                                          <  1:518565/150‑1 (MQ=255)
                                                                                        ggATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGTATAAc                                                       >  2:816085/1‑150 (MQ=255)
                                                                                         gATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGTATAAc                                                       >  1:900380/1‑149 (MQ=255)
                                                                                          aTTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGTATAACaga                                                    <  1:451482/151‑1 (MQ=255)
                                                                                             gCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGTATAGCGgaga                                                  <  1:332554/150‑1 (MQ=255)
                                                                                             gCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGTATAACAgaga                                                  <  2:489812/150‑1 (MQ=255)
                                                                                              cGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTTGCGTTGTTCCATGGGTGCGCCTTTGTCGTAAAGGAAATCGTTATTCTGACGCAAGGCGGGAAGCGAAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGTATAACAGAGAt                                                 <  1:723635/150‑1 (MQ=255)
                                                                                              cGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGTATAACAGAGAt                                                 >  2:817023/1‑150 (MQ=255)
                                                                                                   cATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGTATAACAGAGATAACgg                                            <  1:148020/150‑1 (MQ=255)
                                                                                                    aTGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGTATAACAGAGATAACGggc                                          <  2:383018/151‑1 (MQ=255)
                                                                                                    aTGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGTATAACAGAGATAACGggc                                          >  1:200665/1‑151 (MQ=255)
                                                                                                    aTGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGTATAACAGAGATAACGggc                                          <  1:276367/151‑1 (MQ=255)
                                                                                                      gATACCAGTGGCTTTTGCCGGGGATGAGGCGTTGTTCCCTGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATAGGCGATAACAAATATCAGAAAGGTATAACAGAGATAACGggcgg                                        >  1:1013827/1‑151 (MQ=255)
                                                                                                         aCCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGTATAACAGAGATAACGGGCGGCa                                      <  2:812630/150‑1 (MQ=255)
                                                                                                            aGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGTATAACAGAGATAACGGGCGGCAGAAc                                  <  1:426189/151‑1 (MQ=255)
                                                                                                            aGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGTATAACAGAGATAACGGGCGGCAGAAc                                  <  1:663809/151‑1 (MQ=255)
                                                                                                             gTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGTATAACAGAGATAACGGGCGGCAGAAct                                 <  2:94281/151‑2 (MQ=255)
                                                                                                                    ttGCCGGTGATTTGGCGTTGTTCCATGGGTGCTCCTTGCTCGTAAAGGACATCGTTATCCTTACCCAAGGCGCGAAGGGGACAAAGACGTATCCGCTATAACAAATATCAGAAAGGTATAACAGCGATAACGGGCGGCATACCGCCGccc                           <  1:217825/150‑1 (MQ=255)
                                                                                                                       ccGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGTATAACAGAGATAACGGGCGGCAGAACGCCGCCCATCt                       <  2:950723/151‑1 (MQ=255)
                                                                                                                        cGGTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGTATAACAGAGATAACGGGCGGCAGAACGCCGCCCATCtt                      <  2:945220/151‑1 (MQ=255)
                                                                                                                         ggTGATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGTATAACAGAGATAACGGGCGGCAGAACGCCGCCCATCttt                     >  1:742495/1‑151 (MQ=255)
                                                                                                                            gATGTGGCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGTATAACAGAGATAACGGGCGGCAGAACGCCGCCCATCTTTAcc                  >  2:886680/1‑151 (MQ=255)
                                                                                                                                 ggCGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTAGCCTGACGCAAGGCGAGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGTATAACAGAGATAACGGGCGGCAGAACGCCGCCCATCTTTACCAACa              <  1:24437/150‑1 (MQ=255)
                                                                                                                                   cGTTGTTCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGTATAACAGAGATAACGGGCGGCAGAACGCCGCCCATCTTTACCAACAGaa           <  1:460546/151‑1 (MQ=255)
                                                                                                                                    gttgttCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGTATAACAGAGATAACGGGCGGCAGAACGCCGCCCATCTTTACCAACAGAAc          <  1:346817/151‑1 (MQ=255)
                                                                                                                                        ttCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGTATAACAGAGATAACGGGCGGCAGAACGCCGCCCATCTTTACCAACAGAACGatt      >  2:356382/1‑151 (MQ=255)
                                                                                                                                         tCCATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGTATAACAGAGATAACGGGCGGCAGAACGCCGCCCATCTTTACCAACAGAACGatt      >  1:188923/1‑150 (MQ=255)
                                                                                                                                          ccATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGTATAACAGAGATAACGGGCGGCAGAACGCCGACCATCTTTACCAACAGAACGattat    >  1:944841/1‑151 (MQ=255)
                                                                                                                                           cATGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGTATAACAGAGATAACGGGCGGCAGAACGCCGCCCATCTTTACCAACAGAACGattatt   <  1:279816/151‑1 (MQ=255)
                                                                                                                                            aTGGGTGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGTATAACAGAGATAACGGGCGGCAGAACGCCGCCCATCTTTACCAACAGAACGATTAttt  <  1:940207/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CTTTTTCAAGGATCACGTCGAGTAGAAAGGGATCGCTGACCTTTGCGGCTTCCGGGACCAGAGGCAGAACGTCATACTCCGTAGTACTGGATTGCGTTTCATGATACCAGTGGCTTTTGCCGGTGATGTGGCGTTGTTCCATGGGCGCTCCTTGGTCGTAAAGGAAATCGTTATCCTGACGCAAGGCGGGAAGGGGAGAAAGACGGATCGGCGATAACAAATATCAGAAAGGTATAACAGAGATAACGGGCGGCAGAACGCCGCCCATCTTTACCAACAGAACGATTATTT  >  NC_000913/1349966‑1350256

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