Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 27,186 G→A 5.1% V304M (GTG→ATG)  ispH → 4‑hydroxy‑3‑methylbut‑2‑enyl diphosphate reductase, 4Fe‑4S protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091327,1860GA 5.1% 258.5 / 2.7 78V304M (GTG→ATG) ispH4‑hydroxy‑3‑methylbut‑2‑enyl diphosphate reductase, 4Fe‑4S protein
Reads supporting (aligned to +/- strand):  ref base G (39/35);  new base A (2/2);  total (41/37)
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.51e-01

GTGAAAGAGGTTAAATGCGTCGGCGTGACTGCGGGCGCATCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGAC  >  NC_000913/27042‑27334
                                                                                                                                                |                                                                                                                                                    
gTGAAAGAGGTTAAATGCGTCGGCGTGACTGCGGGCGCATCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGcc                                                                                                                                                  >  1:276590/1‑149 (MQ=255)
gTGAAAGAGGTTAAATGCGTCGGCGTGACTGCGGGCGCATCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGcc                                                                                                                                                  >  3:409015/1‑149 (MQ=255)
gTGAAAGAGGTTAAATGCGTCGGCGTGACTGCGGGCGCATCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGcc                                                                                                                                                  >  1:518853/1‑149 (MQ=255)
gTGAAAGAGGTTAAATGCGTCGGCGTGACTGCGGGCGCATCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGcc                                                                                                                                                  <  4:33747/149‑1 (MQ=255)
   aaaGAGGTTAAATGCGTCGGCGTGACTGCGGGCGCATCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGaa                                                                                                                                               >  3:283404/1‑149 (MQ=255)
   aaaGAGGTTAAATGCGTCGGCGTGACTGCGGGCGCATCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGaa                                                                                                                                               >  1:612934/1‑149 (MQ=255)
    aaGAGGTTAAATGCGTCGGCGTGACTGCGGGCGCATCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGaaa                                                                                                                                              >  4:580049/1‑149 (MQ=255)
           tAAATGCGTCGGCGTGACTGCGGGCGCATCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGc                                                                                                                                       <  4:283404/149‑1 (MQ=255)
           tAAATGCGTCGGCGTGACTGCGGGCGCATCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGc                                                                                                                                       <  4:650016/149‑1 (MQ=255)
           tAAATGCGTCGGCGTGACTGCGGGCGCATCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGc                                                                                                                                       <  4:325865/149‑1 (MQ=255)
              aTGCGTCGGCGTGACTGCGGGCGCATCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCgtg                                                                                                                                    >  1:447278/1‑149 (MQ=255)
               tGCGTCGGCGTGACTGCGGGCGCATCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCgtgt                                                                                                                                   >  3:665213/1‑149 (MQ=255)
                  gTCGGCGTGACTGCGGGCGCATCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGa                                                                                                                                <  3:790595/149‑1 (MQ=255)
                   tCGGCGTGACTGCGGGCGCATCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGat                                                                                                                               <  2:350156/149‑1 (MQ=255)
                      gCGTGACTGCGGGCGCATCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGATGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATAtt                                                                                                                            >  4:567500/1‑149 (MQ=255)
                       cGTGACTGCGGGCGCATCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACATTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGGTGCGTGTCGATATTc                                                                                                                           <  2:114521/149‑1 (MQ=255)
                         tGACTGCGGGCGCATCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGt                                                                                                                         <  3:304730/149‑1 (MQ=255)
                         tGACTGCGGGCGCATCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGt                                                                                                                         >  1:794577/1‑149 (MQ=255)
                             tGCGGGCGCATCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAg                                                                                                                     >  4:240244/1‑149 (MQ=255)
                               cGGGCGCATCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTc                                                                                                                   <  2:57004/149‑1 (MQ=255)
                               cGGGCGCATCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTc                                                                                                                   <  3:580049/149‑1 (MQ=255)
                                    gcATCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTa                                                                                                              >  3:204588/1‑149 (MQ=255)
                                     cATCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTaa                                                                                                             <  2:777234/149‑1 (MQ=255)
                                     cATCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTaa                                                                                                             <  1:375450/149‑1 (MQ=255)
                                      aTCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAg                                                                                                            >  3:253426/1‑149 (MQ=255)
                                      aTCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAg                                                                                                            >  3:479442/1‑149 (MQ=255)
                                      aTCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAg                                                                                                            <  4:288076/149‑1 (MQ=255)
                                      aTCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAg                                                                                                            >  1:191590/1‑149 (MQ=255)
                                       tCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGt                                                                                                           >  4:408839/1‑149 (MQ=255)
                                          gCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCAt                                                                                                        <  3:500389/149‑1 (MQ=255)
                                          gCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCAt                                                                                                        >  3:124673/1‑149 (MQ=255)
                                           cTCCGGCTATTCTGGTGCAGAATGTTGTGGCACGTTTGCAGCGGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCAtt                                                                                                       >  4:672792/1‑149 (MQ=255)
                                           cTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCAtt                                                                                                       <  2:518853/149‑1 (MQ=255)
                                                gATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTagcag                                                                                                  >  2:830959/1‑149 (MQ=255)
                                                       tGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGt                                                                                           >  3:528836/1‑149 (MQ=255)
                                                          tGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAATGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTAt                                                                                        >  4:9313/1‑149 (MQ=255)
                                                              gAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGa                                                                                    >  3:683859/1‑149 (MQ=255)
                                                                      tggCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCg                                                                            <  2:428597/149‑1 (MQ=255)
                                                                       ggCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAATGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTa                                                                                                      <  3:850903/122‑1 (MQ=255)
                                                                       ggCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAATGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTa                                                                                                      >  4:850903/1‑122 (MQ=255)
                                                                         cACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTGGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTc                                                                         >  3:806320/1‑149 (MQ=255)
                                                                          aCGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCt                                                                        >  3:153241/1‑149 (MQ=255)
                                                                          aCGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCt                                                                        <  1:116206/149‑1 (MQ=255)
                                                                           cGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCtt                                                                       <  1:524995/149‑1 (MQ=255)
                                                                            gTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTg                                                                      <  4:204588/149‑1 (MQ=255)
                                                                               tGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTa                                                                   <  2:794577/149‑1 (MQ=255)
                                                                                  agcagcTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCg                                                                >  1:872156/1‑149 (MQ=255)
                                                                                         gggCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCAtttt                                                         >  1:504992/1‑149 (MQ=255)
                                                                                            cGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTa                                                      <  3:664286/149‑1 (MQ=255)
                                                                                            cGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTa                                                      <  4:244502/149‑1 (MQ=255)
                                                                                            cGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTa                                                      <  2:612934/149‑1 (MQ=255)
                                                                                             ggtggtGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTAt                                                     >  1:429972/1‑149 (MQ=255)
                                                                                             ggtggtGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTAt                                                     <  1:508462/149‑1 (MQ=255)
                                                                                             ggtggtGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTAt                                                     <  2:578644/149‑1 (MQ=255)
                                                                                                 gtgAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGgag                                                 >  2:572676/1‑149 (MQ=255)
                                                                                                   gAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGaa                                               >  4:166958/1‑149 (MQ=255)
                                                                                                   gAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGaa                                               >  2:780265/1‑149 (MQ=255)
                                                                                                   gAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGaa                                               <  4:386938/149‑1 (MQ=255)
                                                                                                     aGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGcc                                                                                                                                                  >  4:605052/1‑48 (MQ=255)
                                                                                                     aGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGcc                                                                                                                                                  <  3:605052/48‑1 (MQ=255)
                                                                                                      gCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATg                                                                               >  1:113745/1‑114 (MQ=255)
                                                                                                      gCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATg                                                                               <  2:113745/114‑1 (MQ=255)
                                                                                                      gCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAAc                                            >  4:134266/1‑149 (MQ=255)
                                                                                                      gCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAAc                                            >  2:717293/1‑149 (MQ=255)
                                                                                                        cATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACAt                                          <  1:55160/149‑1 (MQ=255)
                                                                                                         aTTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATg                                         >  4:837590/1‑149 (MQ=255)
                                                                                                          ttCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGc                                        <  2:276590/149‑1 (MQ=255)
                                                                                                           tCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCg                                       >  1:334762/1‑149 (MQ=255)
                                                                                                              gCTGGTAGGCCGTGAAGACAACATTGTTTTCGAAGTGACCAAAGAGCTGCGTGTCGCTATTCGTGAAGTCGATTATGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGTGttt                                    <  3:269943/149‑1 (MQ=255)
                                                                                                                tGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTAc                                  <  4:60113/149‑1 (MQ=255)
                                                                                                                           gaagaaAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTc                       <  2:504992/149‑1 (MQ=255)
                                                                                                                              gaaAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGAt                    >  4:528846/1‑149 (MQ=255)
                                                                                                                              gaaAACATTGTTTTCGAAATGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGAt                    <  3:9313/149‑1 (MQ=255)
                                                                                                                               aaaaCATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATa                   <  4:683859/149‑1 (MQ=255)
                                                                                                                                      tGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGAccc            >  4:12471/1‑149 (MQ=255)
                                                                                                                                        ttttCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCg          >  4:321028/1‑149 (MQ=255)
                                                                                                                                            cGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCAt      <  4:528836/149‑1 (MQ=255)
                                                                                                                                              aaGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTg    >  3:604776/1‑149 (MQ=255)
                                                                                                                                                gTGGCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGAc  <  4:722545/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GTGAAAGAGGTTAAATGCGTCGGCGTGACTGCGGGCGCATCGGCTCCGGATATTCTGGTGCAGAATGTGGTGGCACGTTTGCAGCAGCTGGGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGAC  >  NC_000913/27042‑27334

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