Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,029,177 T→G K599Q (AAG→CAG)  fadB ← fused 3‑hydroxybutyryl‑CoA epimerase/delta(3)‑cis‑delta(2)‑trans‑enoyl‑CoA isomerase/enoyl‑CoA hydratase/3‑hydroxyacyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,029,1770TG98.5% 229.9 / ‑3.9 65K599Q (AAG→CAG) fadBfused 3‑hydroxybutyryl‑CoA epimerase/delta(3)‑cis‑delta(2)‑trans‑enoyl‑CoA isomerase/enoyl‑CoA hydratase/3‑hydroxyacyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (34/30);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

TAATGCCTTCCTCCAGACAGCGCACCACTTCGTTGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTTACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGCGCGTCGATGGCATCGCGGTAATCTTTCTGCATCCGCTGCGGGAAGCCTGCTGCCATGACAGCCTGAGCGTG  >  NC_000913/4029030‑4029321
                                                                                                                                                   |                                                                                                                                                
tAATGCCTTCCTCCAGACAGCGCACCACTTCGTTGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGAcc                                                                                                                                               >  2:313685/1‑151 (MQ=255)
  aTGCCTTCCTCCAGACAGCGCACCACTTCGTTGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCtt                                                                                                                                             <  2:921252/151‑1 (MQ=255)
     ccTTCCTCCAGACAGCGCACCACTTCGTTGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTg                                                                                                                                           <  1:572695/150‑1 (MQ=255)
     ccTTCCTCCAGACAGCGCACCACTTCGTTGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGc                                                                                                                                          >  1:97490/1‑151 (MQ=255)
          ctccAGACAGCGCACCACTTCGTTGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGtct                                                                                                                                     >  2:338441/1‑151 (MQ=255)
             cAGACAGCGCACCACTTCGTTGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGtcttct                                                                                                                                  <  1:917427/151‑1 (MQ=255)
              aGACAGCGCACCACTTCGTTGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTcttctt                                                                                                                                 <  1:896525/151‑1 (MQ=255)
               gACAGCGCACCACTTCGTTGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCttt                                                                                                                                >  1:377450/1‑151 (MQ=255)
                 cAGCGCACCACTTCGTTGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTtat                                                                                                                              >  2:452339/1‑151 (MQ=255)
                     gcACCACTTCGTTGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAAc                                                                                                                           >  1:213671/1‑150 (MQ=255)
                     gcACCACTTCGTTGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACg                                                                                                                          >  2:137380/1‑151 (MQ=255)
                       accacTTCGTTGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGcc                                                                                                                        >  1:551937/1‑151 (MQ=255)
                          acTTCGTTGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGcc                                                                                                                        >  1:904127/1‑148 (MQ=255)
                          acTTCGTTGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCa                                                                                                                       <  1:221790/149‑1 (MQ=255)
                           cTTCGTTGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGa                                                                                                                     >  1:371611/1‑150 (MQ=255)
                            ttCGTTGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGaa                                                                                                                    >  2:211238/1‑150 (MQ=255)
                             tCGTTGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGaa                                                                                                                    <  1:298869/149‑1 (MQ=255)
                              cGTTGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAAc                                                                                                                  <  1:700363/150‑1 (MQ=255)
                                ttGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCg                                                                                                                >  2:83885/1‑150 (MQ=255)
                                ttGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCg                                                                                                                >  1:687497/1‑150 (MQ=255)
                                 tGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCg                                                                                                                >  1:838738/1‑149 (MQ=255)
                                        cGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGttc                                                                                                       <  2:177103/151‑1 (MQ=255)
                                         gggATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGttct                                                                                                      >  2:683909/1‑151 (MQ=255)
                                           gATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGttct                                                                                                      <  1:167225/149‑1 (MQ=255)
                                              catcatGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTg                                                                                                  >  1:891759/1‑150 (MQ=255)
                                                  atGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCaa                                                                                             <  2:60009/151‑1 (MQ=255)
                                                    gCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCaaa                                                                                            <  2:637302/150‑1 (MQ=255)
                                                     cGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGc                                                                                          >  1:357873/1‑151 (MQ=255)
                                                      gggCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCg                                                                                         <  2:594634/151‑1 (MQ=255)
                                                       ggCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCgg                                                                                        <  2:147634/151‑1 (MQ=255)
                                                           aTAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTgg                                                                                    >  1:358186/1‑151 (MQ=255)
                                                                  cttcttCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCAt                                                                                 <  2:508083/147‑1 (MQ=255)
                                                                  cttcttCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCaa                                                                              >  1:396355/1‑150 (MQ=255)
                                                                    tcttcGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTTACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCaaa                                                                             <  1:489990/149‑1 (MQ=255)
                                                                       tcGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAgcg                                                                        >  2:464912/1‑151 (MQ=255)
                                                                        cGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAgcgc                                                                       <  2:227229/151‑1 (MQ=255)
                                                                         gCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGcgcg                                                                      >  1:283634/1‑151 (MQ=255)
                                                                          cTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGCGCGt                                                                     >  2:409377/1‑151 (MQ=255)
                                                                           tGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCCGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGCGCGt                                                                     >  2:776865/1‑150 (MQ=255)
                                                                            gAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGCGCGTc                                                                    <  1:921767/150‑1 (MQ=255)
                                                                             aaaTCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGCGCGTCg                                                                   >  1:157904/1‑150 (MQ=255)
                                                                                    gcTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGCGCGTCGATGGCATc                                                           <  1:526823/151‑1 (MQ=255)
                                                                                    gcTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGCGCGTCGATGGCATc                                                           >  1:885568/1‑151 (MQ=255)
                                                                                       tCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGCGCGTCGATGGCATcgc                                                         <  1:341041/150‑1 (MQ=255)
                                                                                             ggctCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGCGCGTCGATGGCATCGCGGTAATc                                                  <  2:159641/151‑1 (MQ=255)
                                                                                             ggctCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGCGCGTCGATGGCATCGCGGTAATc                                                  <  2:378715/151‑1 (MQ=255)
                                                                                                tCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGCGCGTCGATGGCATCGCGGTAATCttt                                               >  2:448150/1‑151 (MQ=255)
                                                                                                tCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGCGCGTCGATGGCATCGCGGTAATCttt                                               >  1:560599/1‑151 (MQ=255)
                                                                                                   cTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGCGCGTCGATGGCATCGCGGTAATCTTTCt                                             <  2:122046/150‑1 (MQ=255)
                                                                                                   cTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGCGCGTCGATGGCATCGCGGTAATCTTTCt                                             <  2:571786/150‑1 (MQ=255)
                                                                                                     tCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGCGCGTCGATGGCATCGCGGTAATCTTTCTg                                            <  2:567628/149‑1 (MQ=255)
                                                                                                     tCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGCGCGTCGATGGCATCGCGGTAATCTTTCTGc                                           <  1:137378/150‑1 (MQ=255)
                                                                                                     tCTGCCAGCAGGTCTCCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGCGCGTCGATGGCATCGCGGTAATCTTTCTg                                            <  1:682995/149‑1 (MQ=255)
                                                                                                               ggTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGCGCGTCGATGGCATCGCGGTAATCTTTCTGCATCCGCTGCgg                                >  2:557818/1‑151 (MQ=255)
                                                                                                                gTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGCGCGTCGATGGCATCGCGGTAATCTTTCTGCATCCGCTGCggg                               >  2:280825/1‑151 (MQ=255)
                                                                                                                  cTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGCGCGTCGATGGCATCGCGGTAATCTTTCTGCATCCGCTGCGGGa                              >  1:147142/1‑150 (MQ=255)
                                                                                                                    tCAACGGCGGCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGCGCGTCGATGGCATCGCGGTAATCTTTCTGCATCCGCTGCGGGaa                             <  2:527976/149‑1 (MQ=255)
                                                                                                                        cggcggCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGCGCGTCGATGGCATCGCGGTAATCTTTCTGCATCCGCTGCGGGAAGCctgc                       <  2:644065/151‑1 (MQ=255)
                                                                                                                        cggcggCGTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGCGCGTCGATGGCATCGCGGTAATCTTTCTGCATCCGCTGCGGGAAGCctgc                       >  1:314276/1‑151 (MQ=255)
                                                                                                                           cggcgTCTTCTTCTTTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGCGCGTCGATGGCATCGCGGTAATCTTTCTGCATCCGCTGCGGGAAGCCtgctgc                    <  2:25989/151‑1 (MQ=255)
                                                                                                                                    cttcttTCTTCGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGCGCGTCGATGGCATCGCGGTAATCTTTCTGCATCCGCTGCGGGAAGCCTGCTGCCATGACAGc           >  1:39571/1‑151 (MQ=255)
                                                                                                                                          tcttcGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGCGCGTCGATGGCATCGCGGTAATCTTTCTGCATCCGCTGCGGGAAGCCTGCTGCCATGACAGCCTGAg      <  1:842473/150‑1 (MQ=255)
                                                                                                                                            ttcGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGCGCGTCGATGGCATCGCGGTAATCTTTCTGCATCCGCTGCGGGAAGCCTGCTGCCATGACAGCCTGAGCGt   >  1:732806/1‑151 (MQ=255)
                                                                                                                                             tcGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGCGCGTCGATGGCATCGCGGTAATCTTTCTGCATCCGCTGCGGGAAGCCTGCTGCCATGACAGCCTGAGc     <  2:891780/148‑1 (MQ=255)
                                                                                                                                             tcGGCTGACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGCGCGTCGATGGCATCGCGGTAATCTTTCTGCATCCGCTGCGGGAAGCCTGCTGCCATGACAGCCTGAGCGTg  >  2:27494/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TAATGCCTTCCTCCAGACAGCGCACCACTTCGTTGACCATCGGGATCATCATGCGGGCGATAATCTCTTCTTCGCTGAAATCGCGCTTCGGCTGGCTCACTTCTGCCAGCAGGTCTTCAACGGCGGCGTCTTCTTCTTTCTTCGGCTTACCTTTGCTGTCTTCTTTATAACGCCAGAAACCGAGGCCGTTCTTCTGACCAAAGCGGTTGGCATCAAACAGCGCGTCGATGGCATCGCGGTAATCTTTCTGCATCCGCTGCGGGAAGCCTGCTGCCATGACAGCCTGAGCGTG  >  NC_000913/4029030‑4029321

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