Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,395,455 G→A 11.8% A464T (GCA→ACA)  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,4550GA11.8% 263.3 / 19.5 92A464T (GCA→ACA) nnrbifunctional 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 G (39/42);  new base A (7/4);  total (46/46)
Fisher's exact test for biased strand distribution p-value = 5.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAA  >  NC_000913/4395309‑4395600
                                                                                                                                                  |                                                                                                                                                 
gAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCa                                                                                                                                                 >  3:642090/1‑149 (MQ=255)
       ggTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCtgt                                                                                                                                           <  2:484395/148‑1 (MQ=255)
       ggTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCtgt                                                                                                                                           >  1:484395/1‑148 (MQ=255)
         tCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTc                                                                                                                                        <  1:631374/149‑1 (MQ=255)
           gccgccCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTcgc                                                                                                                                      >  3:197526/1‑149 (MQ=255)
           gccgccCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTcgc                                                                                                                                      <  2:423333/149‑1 (MQ=255)
              gccCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCa                                                                                                                                   <  1:237146/149‑1 (MQ=255)
              gccCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCa                                                                                                                                   <  4:622030/149‑1 (MQ=255)
              gccCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCa                                                                                                                                   <  1:430506/149‑1 (MQ=255)
                ccATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACg                                                                                                                                 >  2:216205/1‑149 (MQ=255)
                  aTCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGt                                                                                                                               <  4:405481/149‑1 (MQ=255)
                   tCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTg                                                                                                                              >  2:30583/1‑149 (MQ=255)
                    ccTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGc                                                                                                                             <  1:669174/149‑1 (MQ=255)
                        aCGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCa                                                                                                                         >  1:254264/1‑149 (MQ=255)
                        aCGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCa                                                                                                                         >  1:9989/1‑149 (MQ=255)
                        aCGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCa                                                                                                                         >  2:544141/1‑149 (MQ=255)
                          gCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCGCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGc                                                                                                                       >  4:617124/1‑149 (MQ=255)
                            tttAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGc                                                                                                                             <  1:229541/141‑1 (MQ=255)
                            tttAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGc                                                                                                                             >  2:229541/1‑141 (MQ=255)
                              tAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGAc                                                                                                                   <  4:642090/149‑1 (MQ=255)
                              tAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGAc                                                                                                                   <  1:216205/149‑1 (MQ=255)
                                 gCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCtgt                                                                                                                                           >  4:432486/1‑122 (MQ=255)
                                 gCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCtgt                                                                                                                                           <  3:432486/122‑1 (MQ=255)
                                 gCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTa                                                                                                                >  1:16586/1‑149 (MQ=255)
                                 gCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTa                                                                                                                >  3:628064/1‑149 (MQ=255)
                                    ttattGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCt                                                                                                                                             >  2:450769/1‑117 (MQ=255)
                                    ttattGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCt                                                                                                                                             <  1:450769/117‑1 (MQ=255)
                                       ttGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGCACAGGCtgca                                                                                                                                          <  1:648549/117‑3 (MQ=255)
                                       ttGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGCACAGGCtgca                                                                                                                                          >  2:648549/1‑115 (MQ=255)
                                              cGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGtt                                                                                                   <  3:443998/149‑1 (MQ=255)
                                                gAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGAc                                                                                                                   >  4:337941/1‑131 (MQ=255)
                                                gAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGAc                                                                                                                   <  3:337941/131‑1 (MQ=255)
                                                gAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGtttt                                                                                                 >  4:93566/1‑149 (MQ=255)
                                                gAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGCACAGGCtgcatc                                                                                                                                        >  4:440828/1‑106 (MQ=255)
                                                gAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGCACAGGCtgcatc                                                                                                                                        <  3:440828/110‑5 (MQ=255)
                                                   aTGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGa                                                                                              <  1:544141/149‑1 (MQ=255)
                                                    tGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTg                                                                                                >  3:290382/1‑146 (MQ=255)
                                                    tGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTg                                                                                                <  4:290382/146‑1 (MQ=255)
                                                        ggCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAAcgcg                                                                                         <  4:508551/149‑1 (MQ=255)
                                                           aTGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCgg                                                                                      <  2:669757/149‑1 (MQ=255)
                                                             ggCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGTCGCGTTTTGGAAcgcg                                                                                         >  1:654697/1‑144 (MQ=255)
                                                             ggCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAAcgcg                                                                                         <  2:654697/144‑1 (MQ=255)
                                                                        gcATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAAcc                                                                         >  3:334456/1‑149 (MQ=255)
                                                                        gcATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAAcc                                                                         >  3:253271/1‑149 (MQ=255)
                                                                          aTGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGa                                                                       >  4:65313/1‑149 (MQ=255)
                                                                           tGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGAt                                                                      <  3:236441/149‑1 (MQ=255)
                                                                                       tctctGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACg                                                          >  4:69759/1‑149 (MQ=255)
                                                                                         tctGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCt                                                        >  1:10922/1‑149 (MQ=255)
                                                                                              tattatTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGc                                                   >  2:452026/1‑149 (MQ=255)
                                                                                              tattatTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGc                                                   >  3:528653/1‑149 (MQ=255)
                                                                                               attattGGCGCATTGCTTGGGCAAAAACTGTCGCTGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCg                                                  >  2:667742/1‑149 (MQ=255)
                                                                                                  attGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTAt                                               >  2:677374/1‑149 (MQ=255)
                                                                                                  attGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTAt                                               <  3:617124/149‑1 (MQ=255)
                                                                                                     ggCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGt                                            <  4:528011/149‑1 (MQ=255)
                                                                                                     ggCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGt                                            <  2:9989/149‑1 (MQ=255)
                                                                                                      gcgcATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGtt                                           <  1:530308/149‑1 (MQ=255)
                                                                                                        gcATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCa                                                            >  1:180970/1‑130 (MQ=255)
                                                                                                        gcATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCa                                                            <  2:180970/130‑1 (MQ=255)
                                                                                                        gcATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTaa                                         >  2:15624/1‑149 (MQ=255)
                                                                                                        gcATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTaa                                         <  1:281399/149‑1 (MQ=255)
                                                                                                         cATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAAc                                        <  4:226766/149‑1 (MQ=255)
                                                                                                              cTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGa                                   <  1:610204/149‑1 (MQ=255)
                                                                                                              cTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGa                                   <  2:254264/149‑1 (MQ=255)
                                                                                                               ttGGGCAAAAACTGTCGCTGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGaa                                  <  1:667742/149‑1 (MQ=255)
                                                                                                                 gggCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGt                                >  1:664231/1‑149 (MQ=255)
                                                                                                                   gCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGa                              >  3:175855/1‑149 (MQ=255)
                                                                                                                   gCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGa                              <  2:10922/149‑1 (MQ=255)
                                                                                                                   gCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGa                              >  2:112932/1‑149 (MQ=255)
                                                                                                                     aaaaaCTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACt                            >  4:165237/1‑149 (MQ=255)
                                                                                                                     aaaaaCTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACt                            >  2:559723/1‑149 (MQ=255)
                                                                                                                      aaaaCTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTg                           >  2:345454/1‑149 (MQ=255)
                                                                                                                       aaaCTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGa                          >  3:557296/1‑149 (MQ=255)
                                                                                                                               gCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAAcc                  >  4:341058/1‑149 (MQ=255)
                                                                                                                                ccGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCa                 <  4:197526/149‑1 (MQ=255)
                                                                                                                                   tATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCatga              <  3:69759/149‑1 (MQ=255)
                                                                                                                                   tATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCatga              >  1:452096/1‑149 (MQ=255)
                                                                                                                                   tATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCatga              <  1:575110/149‑1 (MQ=255)
                                                                                                                                       atgCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGaa          <  1:26129/149‑1 (MQ=255)
                                                                                                                                       acacAGCCTGCACAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGaa          >  1:608355/4‑149 (MQ=255)
                                                                                                                                       acacAGCCTGCACAGGCTGCGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGaa          >  4:265024/4‑149 (MQ=255)
                                                                                                                                         acAGCCTGCACAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCatga              >  3:63994/2‑143 (MQ=255)
                                                                                                                                         acAGCCTGCACAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCatga              <  4:63994/142‑1 (MQ=255)
                                                                                                                                         acAGCCTGCACAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATc        >  3:244921/2‑149 (MQ=255)
                                                                                                                                         acAGCCTGCACAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATc        >  1:156691/2‑149 (MQ=255)
                                                                                                                                           aGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGa      <  4:628064/149‑1 (MQ=255)
                                                                                                                                            gCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGGATCGAg     <  2:452096/149‑1 (MQ=255)
                                                                                                                                            gCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATc        <  2:33764/146‑1 (MQ=255)
                                                                                                                                            gCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATc        >  1:33764/1‑146 (MQ=255)
                                                                                                                                            gCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAg     <  2:449296/149‑1 (MQ=255)
                                                                                                                                            gCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAg     <  3:341058/149‑1 (MQ=255)
                                                                                                                                            gCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAg     <  1:15624/149‑1 (MQ=255)
                                                                                                                                            gCCTGCACAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAg     <  4:244921/149‑1 (MQ=255)
                                                                                                                                               tgtgCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACCGATAAAAACCATGATGAATCGAGTaa  >  2:114906/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GAACCGTGGTCGCCGCCCATCCTGACGCTTTAGGCATTATTGATGCCGGAAATGCAGGCATGGCGAGCGGCGGCATGGGCGATGTGCTCTCTGGTATTATTGGCGCATTGCTTGGGCAAAAACTGTCGCCGTATGATGCAGCCTGTGCAGGCTGTGTCGCGCACGGTGCGGCAGCTGACGTACTGGCGGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAA  >  NC_000913/4395309‑4395600

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