Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,395,454 T→C 5.4% C463C (TGT→TGC nnr → bifunctional NAD(P)H‑hydrate repair enzyme; C‑terminal domain ADP‑dependent (S)‑NAD(P)H‑hydrate dehydratase and N‑terminal domain NAD(P)H‑hydrate epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,395,4540TC 5.4% 416.2 / 8.9 129C463C (TGT→TGCnnrbifunctional NAD(P)H‑hydrate repair enzyme; C‑terminal domain ADP‑dependent (S)‑NAD(P)H‑hydrate dehydratase and N‑terminal domain NAD(P)H‑hydrate epimerase
Reads supporting (aligned to +/- strand):  ref base T (63/59);  new base C (2/5);  total (65/64)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 131 total aligned reads displayed.

CGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAA  >  NC_000913/4395307‑4395600
                                                                                                                                                   |                                                                                                                                                  
cGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCtgtg                                                                                                                                                   >  2:1520647/1‑149 (MQ=255)
    aCCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAgg                                                                                                                                               <  4:41070/149‑1 (MQ=255)
     ccGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGc                                                                                                                                              <  4:1403994/149‑1 (MQ=255)
       gTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCtg                                                                                                                                            >  4:1525695/1‑149 (MQ=255)
        tGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATACAGCCTGCACAGGCtgc                                                                                                                                           <  3:647824/149‑2 (MQ=255)
         ggTCGCCGCCCATCCTCACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGCACAGGCtgca                                                                                                                                          <  1:1178399/149‑3 (MQ=255)
          gTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTgtgt                                                                                                                                         >  1:1553089/1‑149 (MQ=255)
          gTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTgtgt                                                                                                                                         >  4:314407/1‑149 (MQ=255)
          gTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGCACAGGCtgcat                                                                                                                                         <  2:756339/149‑4 (MQ=255)
            cgccgcCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTcg                                                                                                                                       <  3:1568657/149‑1 (MQ=255)
                gccCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCa                                                                                                                                   >  4:474386/1‑149 (MQ=255)
                   cATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACgg                                                                                                                                >  2:275160/1‑149 (MQ=255)
                         gACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCtgt                                                                                                                                           >  2:745504/1‑132 (MQ=255)
                         gACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGc                                                                                                                          <  1:374857/149‑1 (MQ=255)
                          aCGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCa                                                                                                                         <  2:1329394/149‑1 (MQ=255)
                                tAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGAc                                                                                                                   <  2:544285/149‑1 (MQ=255)
                                tAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGAc                                                                                                                   <  4:561968/149‑1 (MQ=255)
                                   gCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTa                                                                                                                >  1:817511/1‑149 (MQ=255)
                                    cATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCt                                                                                                                      >  3:1051909/1‑142 (MQ=255)
                                    cATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTAc                                                                                                               <  2:817511/149‑1 (MQ=255)
                                     attattGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACt                                                                                                              >  4:787672/1‑149 (MQ=255)
                                     attattGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACt                                                                                                              >  1:584860/1‑149 (MQ=255)
                                     attattGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACt                                                                                                              <  3:497675/149‑1 (MQ=255)
                                      ttattGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTg                                                                                                             >  2:957060/1‑149 (MQ=255)
                                        attGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTggc                                                                                                           >  1:638202/1‑149 (MQ=255)
                                        attGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTggc                                                                                                           >  2:200883/1‑149 (MQ=255)
                                               ccGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGt                                                                                                    <  4:1203330/149‑1 (MQ=255)
                                                   aaaTGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTAGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTg                                                                                                <  2:528536/149‑1 (MQ=255)
                                                       gCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAAc                                                                                            <  2:1189638/149‑1 (MQ=255)
                                                        cAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAAcg                                                                                           >  1:1094414/1‑149 (MQ=255)
                                                                gCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGGGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGAt                                                                                   <  3:1525695/149‑1 (MQ=255)
                                                                  gAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGc                                                                                 >  1:1238314/1‑149 (MQ=255)
                                                                  gAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGc                                                                                 >  1:1176138/1‑149 (MQ=255)
                                                                    gcggcgGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTg                                                                               <  3:557359/149‑1 (MQ=255)
                                                                    gcggcgGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTg                                                                               <  2:38763/149‑1 (MQ=255)
                                                                               ggCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGAtct                                                                    >  4:622136/1‑149 (MQ=255)
                                                                                gCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGAtctc                                                                   <  3:635194/149‑1 (MQ=255)
                                                                                 cGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATctct                                                                  <  1:534669/149‑1 (MQ=255)
                                                                                  gATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCtt                                                                 >  3:223953/1‑149 (MQ=255)
                                                                                  gATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCtt                                                                 >  3:675734/1‑149 (MQ=255)
                                                                                   aTGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCttt                                                                <  2:1452234/149‑1 (MQ=255)
                                                                                    tgtgCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCtttt                                                               <  2:1176138/149‑1 (MQ=255)
                                                                                      tgCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTcc                                                             >  4:848974/1‑149 (MQ=255)
                                                                                       gCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCgg                                                                                      <  3:15101/123‑1 (MQ=255)
                                                                                         tctctGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCaa                                                                           >  3:1481734/1‑132 (MQ=255)
                                                                                         tctctGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCaa                                                                           <  4:1481734/132‑1 (MQ=255)
                                                                                         tctctGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACg                                                          >  2:1374331/1‑149 (MQ=255)
                                                                                             tGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTAc                                                      >  1:202524/1‑149 (MQ=255)
                                                                                             tGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTAc                                                      >  4:305708/1‑149 (MQ=255)
                                                                                               gTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAg                                                    >  4:931576/1‑149 (MQ=255)
                                                                                                tattatTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGc                                                   >  1:399533/1‑149 (MQ=255)
                                                                                                tattatTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGc                                                   >  2:48010/1‑149 (MQ=255)
                                                                                                   tattGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTa                                                <  1:1444909/149‑1 (MQ=255)
                                                                                                     ttGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTAtt                                              <  3:787672/149‑1 (MQ=255)
                                                                                                      tGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTg                                             <  2:1553089/149‑1 (MQ=255)
                                                                                                        gcgcATTTCTTGGGCAAAACCTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGGGGCGCGTTTTGGAACGCGCGGGATGCTGGGAACCGATCGCTTTTCCCCGCTACAGCGTATTGtt                                           <  2:1067724/149‑1 (MQ=255)
                                                                                                        gcgcATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGtt                                           >  3:1213615/1‑149 (MQ=255)
                                                                                                          gcATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTaa                                         >  1:196606/1‑149 (MQ=255)
                                                                                                          gcATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTaa                                         <  2:202524/149‑1 (MQ=255)
                                                                                                          gcATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTaa                                         >  2:150373/1‑149 (MQ=255)
                                                                                                            aTTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGc                                                                                                                                                  >  4:258803/1‑42 (MQ=255)
                                                                                                            aTTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGc                                                                                                                                                  <  3:258803/42‑1 (MQ=255)
                                                                                                               gCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCgg                                    >  2:641385/1‑149 (MQ=255)
                                                                                                               gCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCgg                                    >  2:433712/1‑149 (MQ=255)
                                                                                                                cTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGa                                   >  3:1426505/1‑149 (MQ=255)
                                                                                                                 ttGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTTCAGGCTGTGTCGCGCACGTTGCGGCACATTACGTAGTGGCGTCGCGTTGTGGAACGCGCGGTAAGCTGGCAACCGATCCCTTTTCCACGCTACAGCGTAGTGTTAACCCGGaa                                  >  2:393488/1‑149 (MQ=255)
                                                                                                                     gCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGa                              >  4:244211/1‑149 (MQ=255)
                                                                                                                     gCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGa                              <  2:1372753/149‑1 (MQ=255)
                                                                                                                     gCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGa                              <  3:314407/149‑1 (MQ=255)
                                                                                                                     gCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACCGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGa                              >  3:819620/1‑149 (MQ=255)
                                                                                                                      cAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGGAGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCCGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGAc                             <  2:432730/149‑1 (MQ=255)
                                                                                                                       aaaaaCTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGa                                                                                              <  3:519159/83‑1 (MQ=255)
                                                                                                                       aaaaaCTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGTCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACt                            >  1:1099977/1‑149 (MQ=255)
                                                                                                                       aaaaaCTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGt                                            >  4:926530/1‑133 (MQ=255)
                                                                                                                       aaaaaCTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGt                                            <  3:926530/133‑1 (MQ=255)
                                                                                                                       aaaaaCTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACt                            >  1:17391/1‑149 (MQ=255)
                                                                                                                       aaaaaCTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACt                            >  1:90369/1‑149 (MQ=255)
                                                                                                                       aaaaaCTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACt                            >  3:39749/1‑149 (MQ=255)
                                                                                                                        aaaaCTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTg                           <  4:223953/149‑1 (MQ=255)
                                                                                                                         aaaCTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGa                          <  4:1167143/149‑1 (MQ=255)
                                                                                                                          aaCTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGAt                         >  3:791204/1‑149 (MQ=255)
                                                                                                                           aCTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATa                        >  4:37399/1‑149 (MQ=255)
                                                                                                                               tCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATaaaaa                    <  2:17391/149‑1 (MQ=255)
                                                                                                                                 gCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGTCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAAcc                  >  2:49735/1‑149 (MQ=255)
                                                                                                                                   cGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCat                <  3:153375/149‑1 (MQ=255)
                                                                                                                                   cGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCat                >  3:728387/1‑149 (MQ=255)
                                                                                                                                      atgatgCAGCCTGTGCAGGCTGTGACGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTAGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCatgat             <  3:1378969/149‑1 (MQ=255)
                                                                                                                                        gacacAGCCTGCACAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCAtgatga           >  1:655426/5‑149 (MQ=255)
                                                                                                                                           gCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATc        >  4:1307707/1‑149 (MQ=255)
                                                                                                                                           gCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATc        <  1:827024/149‑1 (MQ=255)
                                                                                                                                           acAGCCTGCACAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGaa                                  <  2:919001/122‑1 (MQ=255)
                                                                                                                                             aGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGa      <  3:244211/149‑1 (MQ=255)
                                                                                                                                             aGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACCGATAAAAACCATGATGAATCGa      <  2:578237/149‑1 (MQ=255)
                                                                                                                                              gCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCAtgatga           >  3:114915/1‑143 (MQ=255)
                                                                                                                                              gCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCAtgatga           <  4:114915/143‑1 (MQ=255)
                                                                                                                                              gCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAg     <  3:474386/149‑1 (MQ=255)
                                                                                                                                              gCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAg     <  3:850300/149‑1 (MQ=255)
                                                                                                                                              gCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAg     <  1:139746/149‑1 (MQ=255)
                                                                                                                                                cTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTa   >  3:688386/1‑149 (MQ=255)
                                                                                                                                                 tgtgCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTaa  <  4:404589/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGGAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAA  >  NC_000913/4395307‑4395600

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