Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 368,939 A→G T369A (ACG→GCG)  mhpA → 3‑(3‑hydroxyphenyl)propionate hydroxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913368,9390AG98.3% 696.9 / 0.7 233T369A (ACG→GCG) mhpA3‑(3‑hydroxyphenyl)propionate hydroxylase
Reads supporting (aligned to +/- strand):  ref base A (0/4);  new base G (105/124);  total (105/128)
Fisher's exact test for biased strand distribution p-value = 1.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.44e-01

Only 100 of 260 total aligned reads displayed.

TCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGATGCCGCAA  >  NC_000913/368848‑369028
                                                                                           |                                                                                         
 caacaaGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGtgcg                                                                                         >  1:6262403/1‑93 (MQ=255)
 caacaaGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGtgcg                                                                                         >  1:8482866/1‑93 (MQ=255)
  agcaaGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGtgcgt                                                                                        <  2:5356202/91‑1 (MQ=255)
  aacaaGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGtgcgt                                                                                        <  1:1403129/93‑1 (MQ=255)
  aacaaGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGtgcgt                                                                                        >  1:7660061/1‑93 (MQ=255)
  aacaaGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCCAAACGCTGGCAGGGtgcgt                                                                                        >  2:5588779/1‑93 (MQ=255)
   acaaGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGtgcgtt                                                                                       <  2:1901425/93‑1 (MQ=255)
     aaGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGtgcgttac                                                                                     >  1:5364330/1‑93 (MQ=255)
      aGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGggcgttacg                                                                                    >  2:5233830/1‑93 (MQ=255)
         aCGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGa                                                                                 >  2:8627901/1‑93 (MQ=255)
          cGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGGGCGGTACGTGAc                                                                                >  1:4491465/1‑93 (MQ=255)
            tCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACgg                                                                              <  1:3961743/93‑1 (MQ=255)
               cgATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGt                                                                           >  2:9619230/1‑93 (MQ=255)
                gATCATGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGtc                                                                          <  2:5300215/93‑2 (MQ=255)
                gATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACgggggt                                                                          >  2:6474393/1‑89 (MQ=255)
                gATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGtt                                                                          <  1:2052272/93‑1 (MQ=255)
                gATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGtt                                                                          >  2:7712295/1‑93 (MQ=255)
                 aTCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGttt                                                                         >  2:4165842/1‑93 (MQ=255)
                 aTCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGttt                                                                         <  2:981669/93‑1 (MQ=255)
                 aTCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGGGCGTTAAGTGACGGCGttt                                                                         >  1:1099898/1‑93 (MQ=255)
                  tCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTc                                                                        <  2:1370887/93‑1 (MQ=255)
                  tCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTc                                                                        >  1:7041160/1‑93 (MQ=255)
                   cACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTcc                                                                       <  1:9961344/93‑1 (MQ=255)
                   cACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTcc                                                                       <  1:3443477/93‑1 (MQ=255)
                   cACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTcc                                                                       <  1:6114965/93‑1 (MQ=255)
                   cACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTcc                                                                       <  1:9323950/93‑1 (MQ=255)
                    aCGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCt                                                                      <  2:6647458/93‑1 (MQ=255)
                     cGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTg                                                                     <  1:755759/93‑1 (MQ=255)
                       ccAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGc                                                                   <  1:1693050/93‑1 (MQ=255)
                         aaaGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTg                                                                 >  2:1505463/1‑93 (MQ=255)
                         aaaGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTg                                                                 >  2:3154206/1‑93 (MQ=255)
                          aaGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGt                                                                >  1:3638035/1‑93 (MQ=255)
                          aaGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGGGCGTTACGTGACGGCGTTTCCTGGCTGt                                                                >  2:4324654/1‑93 (MQ=255)
                           aGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGtt                                                               >  2:398535/1‑93 (MQ=255)
                           aGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGtt                                                               <  2:3302134/93‑1 (MQ=255)
                             cGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGa                                                             <  2:5486696/93‑1 (MQ=255)
                             cGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGGGCGTTACGTGACGGCGTTTCCTGGCGGttta                                                             >  2:7917886/1‑91 (MQ=255)
                              gatgatTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGaa                                                            <  1:924552/93‑1 (MQ=255)
                               atgatTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAAt                                                           >  1:2735487/1‑93 (MQ=255)
                               atgatTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAAt                                                           <  2:161201/93‑1 (MQ=255)
                                 gatTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTa                                                         <  1:5899119/93‑1 (MQ=255)
                                   ttGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATc                                                       <  1:5549262/93‑1 (MQ=255)
                                   ttGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATc                                                       <  1:7970663/93‑1 (MQ=255)
                                   ttGACCTGTCCGTGACGGCGGACAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATc                                                       >  1:1000026/1‑93 (MQ=255)
                                    tGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCt                                                      <  1:6472728/93‑1 (MQ=255)
                                      aCCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTgc                                                    <  1:9720462/93‑1 (MQ=255)
                                         tGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTgccgc                                                 >  2:3436597/1‑93 (MQ=255)
                                         tGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGGGCGCTACGTGACGGCGTTTCCTGGATGTTGGATTATCcgccgc                                                 >  2:1731027/1‑93 (MQ=255)
                                           tCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCa                                               <  2:1073666/93‑1 (MQ=255)
                                           tCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCa                                               >  1:295239/1‑93 (MQ=255)
                                            ctgTGACGGCGGGCAACTTGCTGGCTCCTCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCCTTTCCTGGCTGTTGAATTATCTGCCGCCAg                                              <  1:4524490/91‑1 (MQ=255)
                                            ccGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCTAg                                              >  1:8685770/1‑93 (MQ=255)
                                            ccGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAg                                              >  2:4441830/1‑93 (MQ=255)
                                              gTGAGGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTa                                            <  2:7210084/93‑1 (MQ=255)
                                              gTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTa                                            >  1:4642625/1‑93 (MQ=255)
                                              gTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTa                                            <  1:4923822/93‑1 (MQ=255)
                                                 aCGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTaaaa                                         <  1:1259446/93‑1 (MQ=255)
                                                 aCGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTaaaa                                         <  1:6761625/93‑1 (MQ=255)
                                                 aCGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTaaaa                                         <  1:2144135/93‑1 (MQ=255)
                                                   ggcggGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACg                                       <  2:7023324/93‑1 (MQ=255)
                                                    gcggGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGc                                      <  1:1963232/93‑1 (MQ=255)
                                                      gggCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTa                                    <  2:2161615/93‑1 (MQ=255)
                                                       ggCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTAc                                   <  1:2675628/93‑1 (MQ=255)
                                                       ggCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTAc                                   <  1:8902751/93‑1 (MQ=255)
                                                       ggCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTAc                                   >  2:1152193/1‑93 (MQ=255)
                                                        gCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACt                                  <  1:7478408/93‑1 (MQ=255)
                                                         cAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACtt                                 <  2:2865673/93‑1 (MQ=255)
                                                         cAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACtt                                 >  2:9048336/1‑93 (MQ=255)
                                                          aaCGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTc                                >  2:5765982/1‑93 (MQ=255)
                                                          aaCGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTc                                >  1:6922795/1‑93 (MQ=255)
                                                           aCGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTcc                               >  1:5600251/1‑93 (MQ=255)
                                                           aCGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTcc                               <  2:974162/93‑1 (MQ=255)
                                                           aCGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTcc                               >  1:9604044/1‑93 (MQ=255)
                                                            cGTGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCt                              <  1:2237355/93‑1 (MQ=255)
                                                              tGCTGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCg                            >  2:6096054/1‑93 (MQ=255)
                                                                 tGGCTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGaaa                         >  2:935666/1‑93 (MQ=255)
                                                                    cTCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATgc                      <  2:2093782/93‑1 (MQ=255)
                                                                     tCCGCCGAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGCAATgcg                     >  2:8902093/1‑93 (MQ=255)
                                                                      ccgccgAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATgcgc                    <  1:7391545/93‑1 (MQ=255)
                                                                      ccgccgAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCCAAATgcgc                    >  1:9722799/1‑93 (MQ=255)
                                                                       cgccgAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCt                   <  1:5760837/93‑1 (MQ=255)
                                                                       cgccgAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCt                   >  2:4774480/1‑93 (MQ=255)
                                                                         ccgAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTAGAAATGCGCTTc                 >  2:5980865/1‑93 (MQ=255)
                                                                          cgAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCa                <  2:5807355/93‑1 (MQ=255)
                                                                           gAAACGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCaa               >  1:9133285/1‑93 (MQ=255)
                                                                           gAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCaa               <  2:5512343/93‑1 (MQ=255)
                                                                              aCGCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGcc            <  1:8160049/93‑1 (MQ=255)
                                                                                gCTGGCAGGGTGCGTTACGTGACGTCGTTTCCTGGCTGTTAAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGa          <  1:3299357/93‑1 (MQ=255)
                                                                                gCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGa          <  1:4325934/93‑1 (MQ=255)
                                                                                gCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGa          <  1:4202882/93‑1 (MQ=255)
                                                                                gCTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGa          <  1:6560705/93‑1 (MQ=255)
                                                                                 cTGGCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGAt         >  2:3916867/1‑93 (MQ=255)
                                                                                   ggCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGATGc       >  1:7597699/1‑93 (MQ=255)
                                                                                    gCAGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGATGcc      <  2:4739935/93‑1 (MQ=255)
                                                                                      aGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTACAACGCTACTTCCTCGAAATGCGCTTCAAGCCTATGCCGc    >  2:7002678/1‑93 (MQ=255)
                                                                                      aGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGATGCCGc    <  1:3029232/93‑1 (MQ=255)
                                                                                      aGGGTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGATGCCGc    >  2:6269011/1‑93 (MQ=255)
                                                                                       gggTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGATGCCGCa   >  1:2403512/1‑93 (MQ=255)
                                                                                        ggTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGATGCCGCaa  <  2:2519517/93‑1 (MQ=255)
                                                                                        ggTGCGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGATGCCGCaa  <  2:8747902/93‑1 (MQ=255)
                                                                                           |                                                                                         
TCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGATGCCGCAA  >  NC_000913/368848‑369028

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