Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,621,622 G→A 7.6% intergenic (+20/‑147) deoD → / → yjjJ purine nucleoside phosphorylase 1; nicotinamide 1‑beta‑D‑riboside phosphorylase/putative protein kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,621,6220GA 7.6% 204.7 / 6.0 66intergenic (+20/‑147)deoD/yjjJpurine nucleoside phosphorylase 1; nicotinamide 1‑beta‑D‑riboside phosphorylase/putative protein kinase
Reads supporting (aligned to +/- strand):  ref base G (29/32);  new base A (2/3);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

CTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAAT  >  NC_000913/4621474‑4621770
                                                                                                                                                    |                                                                                                                                                    
cTGACCATCTGCACCGTATCTGACCACATCCGCATTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCg                                                                                                                                                      >  4:497163/1‑149 (MQ=255)
  gACCATCTGCACCGTATCTGACCACATCCGCACTCAGGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGAt                                                                                                                                                    >  4:366035/1‑149 (MQ=255)
   aCCATCTGCACCGTATCTGACCACATCCGCATTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGAtt                                                                                                                                                   <  3:497163/149‑1 (MQ=255)
   aCCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCAAtc                                                                                                                                                   <  4:140629/149‑2 (MQ=255)
    ccATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTg                                                                                                                                                  >  3:369282/1‑149 (MQ=255)
     cATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGc                                                                                                                                                 >  3:110500/1‑149 (MQ=255)
          gCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTtgt                                                                                                                                            <  4:208567/149‑1 (MQ=255)
           cACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTtgtg                                                                                                                                           <  2:331998/149‑1 (MQ=255)
            aCCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCtt                                                                                                                                              >  3:30992/1‑145 (MQ=255)
            aCCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCtt                                                                                                                                              <  4:30992/145‑1 (MQ=255)
            aCCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGa                                                                                                                                          >  1:63437/1‑149 (MQ=255)
                tATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGcc                                                                                                                                      <  4:150437/149‑1 (MQ=255)
                    tGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAg                                                                                                                                  <  4:535720/149‑1 (MQ=255)
                       ccACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCgg                                                                                                                               <  4:257420/149‑1 (MQ=255)
                        cacaTCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCggg                                                                                                                              >  3:217531/1‑149 (MQ=255)
                           aTCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGgaga                                                                                                                           <  4:344130/149‑1 (MQ=255)
                           aTCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGgaga                                                                                                                           >  3:460545/1‑149 (MQ=255)
                           aTCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGgaga                                                                                                                           <  3:457315/149‑1 (MQ=255)
                               gCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGc                                                                                                                       >  1:23788/1‑149 (MQ=255)
                                 aCTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTc                                                                                                                     <  3:85301/149‑1 (MQ=255)
                                       gAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCgg                                                                                                                                    >  2:109683/1‑128 (MQ=255)
                                       gAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCgg                                                                                                                                    <  1:109683/128‑1 (MQ=255)
                                        aGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCt                                                                                                                      <  1:482671/141‑1 (MQ=255)
                                        aGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCt                                                                                                                      >  2:482671/1‑141 (MQ=255)
                                         gCAGACCACTGCCGCTGAGAGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCTCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGGCTGCTCCGGCtttt                                                                                                             <  4:217531/149‑1 (MQ=255)
                                           aGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTa                                                                                                           <  3:349038/149‑1 (MQ=255)
                                             aCCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGt                                                                                                         >  4:31255/1‑149 (MQ=255)
                                                 cTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCt                                                                                                     <  4:460545/149‑1 (MQ=255)
                                                 cTGACGCTCAGCGGCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCt                                                                                                     >  3:116473/1‑149 (MQ=255)
                                                     cGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCtttt                                                                                                             >  1:151646/1‑137 (MQ=255)
                                                     cGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCtttt                                                                                                             <  2:151646/137‑1 (MQ=255)
                                                       cTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTCTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCAt                                                                                               <  3:514176/149‑1 (MQ=255)
                                                        tGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCAtt                                                                                              >  4:471837/1‑149 (MQ=255)
                                                           gCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCAttttt                                                                                           <  2:116440/149‑1 (MQ=255)
                                                           gCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCAttttt                                                                                           <  4:119201/149‑1 (MQ=255)
                                                            cGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTtc                                                                                          <  3:59798/149‑1 (MQ=255)
                                                                     aCCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCt                                                                                 <  1:390775/149‑1 (MQ=255)
                                                                        ttCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGa                                                                              >  1:436624/1‑149 (MQ=255)
                                                                           aaCGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTTTTCATTTTTCTCTCCAGCTTGAata                                                                           >  4:23544/1‑149 (MQ=255)
                                                                                   gATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGAtt                                                                                                                                                   <  4:549992/69‑1 (MQ=255)
                                                                                   gATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGAtt                                                                                                                                                   >  3:549992/1‑69 (MQ=255)
                                                                                   gATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCt                                                                   >  2:265327/1‑149 (MQ=255)
                                                                                    aTCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCta                                                                  <  4:504541/149‑1 (MQ=255)
                                                                                     tCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCtat                                                                 >  1:93609/1‑149 (MQ=255)
                                                                                             gCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCa                                                         <  3:254865/149‑1 (MQ=255)
                                                                                              cACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAg                                                        >  3:521525/1‑149 (MQ=255)
                                                                                                 tGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCCGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGa                                                     <  2:63437/149‑1 (MQ=255)
                                                                                                           tCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTtaa                                           <  4:116473/149‑1 (MQ=255)
                                                                                                                tgGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAAt                                      <  4:55423/149‑1 (MQ=255)
                                                                                                                tgGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAAt                                      >  3:55423/1‑149 (MQ=255)
                                                                                                                 gggCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATc                                     <  3:471837/149‑1 (MQ=255)
                                                                                                                 gggCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATc                                     >  1:527941/1‑149 (MQ=255)
                                                                                                                     gATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTAtt                                 >  3:231110/1‑149 (MQ=255)
                                                                                                                       tAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCa                               >  2:213383/1‑149 (MQ=255)
                                                                                                                                 ttGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTAtt                     <  2:406529/149‑1 (MQ=255)
                                                                                                                                    tgtTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCAc                  <  4:486356/149‑1 (MQ=255)
                                                                                                                                    tgtTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCAc                  >  3:424665/1‑149 (MQ=255)
                                                                                                                                      tttCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGa                >  1:278315/1‑149 (MQ=255)
                                                                                                                                         cGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATc             <  4:369282/149‑1 (MQ=255)
                                                                                                                                          ttcacAAGGCAATCGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCt            >  1:459344/5‑149 (MQ=255)
                                                                                                                                           tcacAAGGCAATCGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTa           <  2:459344/146‑1 (MQ=255)
                                                                                                                                              cAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTc        <  2:527941/149‑1 (MQ=255)
                                                                                                                                              cAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTc        >  2:236962/1‑149 (MQ=255)
                                                                                                                                               aaGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCa       >  4:548976/1‑149 (MQ=255)
                                                                                                                                               aaGGCAATCGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATaaa                                       >  1:484761/1‑117 (MQ=255)
                                                                                                                                               aaGGCAATCGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATaaa                                       <  2:484761/117‑1 (MQ=255)
                                                                                                                                                aGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCAt      >  3:145085/1‑149 (MQ=255)
                                                                                                                                                    gATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAAt  >  3:490920/1‑149 (MQ=255)
                                                                                                                                                    gATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAAt  <  2:278315/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAAT  >  NC_000913/4621474‑4621770

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