Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 368,947 T→C R371R (CGT→CGC mhpA → 3‑(3‑hydroxyphenyl)propionate hydroxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913368,9470TC83.8% 566.5 / 99.9 235R371R (CGT→CGCmhpA3‑(3‑hydroxyphenyl)propionate hydroxylase
Reads supporting (aligned to +/- strand):  ref base T (20/18);  new base C (83/114);  total (103/132)
Fisher's exact test for biased strand distribution p-value = 2.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.91e-01

Only 100 of 276 total aligned reads displayed.

AACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGATGCCGCAATATTACGG  >  NC_000913/368856‑369036
                                                                                           |                                                                                         
  cGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGAc                                                                                        >  1:3861253/1‑93 (MQ=255)
  cGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGAc                                                                                        <  2:5212769/93‑1 (MQ=255)
    tCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGGACGTTACGTGACgg                                                                                      >  2:2329692/1‑93 (MQ=255)
      gcgATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCg                                                                                    <  2:2527091/93‑1 (MQ=255)
      gcgATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGGTACGCGACGGCg                                                                                    >  2:4715957/1‑93 (MQ=255)
      gcgATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGGACGTTACGCGACGGCg                                                                                    >  2:5103571/1‑93 (MQ=255)
        gATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGtt                                                                                  <  1:4040056/93‑1 (MQ=255)
         aTCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGttt                                                                                 >  1:2089814/1‑93 (MQ=255)
         aTCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGttt                                                                                 >  1:1500807/1‑93 (MQ=255)
         aTCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGttt                                                                                 >  1:3225811/1‑93 (MQ=255)
          tCACGCGAAGGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTTCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTc                                                                                <  1:3300118/93‑1 (MQ=255)
          tCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTc                                                                                >  2:6053646/1‑93 (MQ=255)
          tCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTc                                                                                <  1:2607986/93‑1 (MQ=255)
          tCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTc                                                                                <  2:2473301/93‑1 (MQ=255)
          tCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTc                                                                                >  1:835498/1‑93 (MQ=255)
             cGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTg                                                                             <  2:1522390/93‑1 (MQ=255)
             cGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTg                                                                             <  2:5161541/93‑1 (MQ=255)
             cGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTg                                                                             <  2:2759674/93‑1 (MQ=255)
              gCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTgg                                                                            >  2:433235/1‑93 (MQ=255)
                 aaaGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGGTCCGCCGAAACGCTGGCAGGGTACGGTACGCGACGGCGTTTTCTGGcgg                                                                         >  1:288548/1‑91 (MQ=255)
                 aaaGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTg                                                                         >  1:3633200/1‑93 (MQ=255)
                 aaaGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGGACGTTACGTGACGGCGTTTCCTGGCTg                                                                         >  1:5933836/1‑93 (MQ=255)
                  aaGCGATGATTGACCTGTCCGTGACGGCGGGAAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGt                                                                        <  2:562624/93‑1 (MQ=255)
                   aGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGtt                                                                       <  1:2930271/93‑1 (MQ=255)
                    gCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTg                                                                      <  1:2559410/93‑1 (MQ=255)
                     cGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGa                                                                     >  1:3709296/1‑93 (MQ=255)
                       atgatTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAAt                                                                   >  2:6473947/1‑93 (MQ=255)
                       atgatTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGTCGTTTCCTGGCTGTTGAAt                                                                   <  1:4984092/93‑1 (MQ=255)
                       atgatTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAAt                                                                   <  1:1161882/93‑1 (MQ=255)
                           ttGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATc                                                               <  2:593243/93‑1 (MQ=255)
                           ttGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGCTACGCGACGGCGTTTCCTGGCTGTTGAATTATc                                                               >  2:2662150/1‑93 (MQ=255)
                            tGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCt                                                              <  2:285393/93‑1 (MQ=255)
                             gACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTg                                                             <  2:1087498/93‑1 (MQ=255)
                              aCCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTgc                                                            <  1:2405149/93‑1 (MQ=255)
                               ccTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTgcc                                                           <  2:4747211/93‑1 (MQ=255)
                                cTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTgccg                                                          >  2:13953/1‑93 (MQ=255)
                                 tGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTgccgc                                                         >  2:6266717/1‑93 (MQ=255)
                                 tGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGATTCCTGGCTGTTGAATTATCTgccgc                                                         >  1:5355425/1‑93 (MQ=255)
                                   tCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCa                                                       <  2:3489312/93‑1 (MQ=255)
                                   tCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCa                                                       >  2:3593718/1‑93 (MQ=255)
                                    ccGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAg                                                      <  2:3961945/93‑1 (MQ=255)
                                     cGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGATTCCTGGCTGTTGAATTATCTGCCGCCAga                                                     >  1:442279/1‑92 (MQ=255)
                                        gACGGCGGGCACCGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTaaa                                                  <  2:5497796/93‑1 (MQ=255)
                                         aCGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTaaaa                                                 <  1:4804919/93‑1 (MQ=255)
                                          cggcggGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAAc                                                <  2:3522990/93‑1 (MQ=255)
                                            gcggGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGGACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCCGTAAAACGc                                              >  2:3471321/1‑93 (MQ=255)
                                             ggggCAAGGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCt                                             <  1:5989738/92‑1 (MQ=255)
                                             cggGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCt                                             <  2:4235978/93‑1 (MQ=255)
                                             cggGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCt                                             <  1:894716/93‑1 (MQ=255)
                                               ggCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTAc                                           >  1:1704852/1‑93 (MQ=255)
                                               ggCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTAc                                           <  2:1621008/93‑1 (MQ=255)
                                               ggCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTAc                                           >  2:1739022/1‑93 (MQ=255)
                                                 cAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACtt                                         <  1:4920304/93‑1 (MQ=255)
                                                   aCGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTcc                                       <  1:1006590/93‑1 (MQ=255)
                                                   aCGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTcc                                       <  2:6422656/93‑1 (MQ=255)
                                                   aCGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTcc                                       >  1:413407/1‑93 (MQ=255)
                                                   aCGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTcc                                       <  2:85485/93‑1 (MQ=255)
                                                   aCGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTcc                                       <  1:948507/93‑1 (MQ=255)
                                                     gTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTc                                     <  2:1279039/93‑1 (MQ=255)
                                                        cTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGaa                                  >  2:3326075/1‑93 (MQ=255)
                                                        cTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGaa                                  <  2:5282957/93‑1 (MQ=255)
                                                         tGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGaaa                                 <  2:1727550/93‑1 (MQ=255)
                                                         tGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGaaa                                 <  2:4524215/93‑1 (MQ=255)
                                                          ggCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAAt                                >  1:515957/1‑93 (MQ=255)
                                                          ggCTCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAAt                                >  1:6486760/1‑93 (MQ=255)
                                                             tCCGCCGAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATgcg                             >  2:5771284/1‑93 (MQ=255)
                                                              ccgccgAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATgcgc                            <  1:5671159/93‑1 (MQ=255)
                                                              ccgccgAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATgcgc                            <  1:1443002/93‑1 (MQ=255)
                                                              ccgccgAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATgcgc                            <  2:3744894/93‑1 (MQ=255)
                                                               cgccgAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCt                           <  1:4659999/93‑1 (MQ=255)
                                                                   gAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCaa                       >  2:6037445/1‑93 (MQ=255)
                                                                   gAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCaa                       >  1:5795658/1‑93 (MQ=255)
                                                                   gAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCat                       >  1:6269474/1‑92 (MQ=255)
                                                                   gAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCaa                       >  1:4733751/1‑93 (MQ=255)
                                                                   gAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCaa                       <  2:5689503/93‑1 (MQ=255)
                                                                   gAAACGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCaa                       <  2:719574/93‑1 (MQ=255)
                                                                     aaCGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGc                     >  2:1434206/1‑93 (MQ=255)
                                                                     aaCGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGc                     >  1:1656584/1‑93 (MQ=255)
                                                                      aCGCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGGTTCAAGcc                    >  1:1791497/1‑93 (MQ=255)
                                                                       cGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCg                   <  1:5598157/93‑1 (MQ=255)
                                                                       cGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCg                   <  2:6635976/93‑1 (MQ=255)
                                                                        gCTGGCAGGTTACGTTACGCGACGGCGTTTCCTGTCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGa                  <  2:2943875/93‑1 (MQ=255)
                                                                        gCTGGCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGa                  <  1:3951156/93‑1 (MQ=255)
                                                                          tGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGATg                >  1:2517689/1‑93 (MQ=255)
                                                                           tgCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGATGc               >  1:6382081/2‑93 (MQ=255)
                                                                           ggCAGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGATGc               <  2:4427192/93‑1 (MQ=255)
                                                                              aGGGTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGATGCCGc            >  1:4519456/1‑93 (MQ=255)
                                                                               gggTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGATGCCGCa           >  1:3798250/1‑93 (MQ=255)
                                                                               gggTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGATGCCGCa           >  1:4379187/1‑93 (MQ=255)
                                                                                ggTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGATGCCGCaa          <  1:3420651/93‑1 (MQ=255)
                                                                                ggTACGTTACGCGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGATGCCGCaa          <  2:6502051/93‑1 (MQ=255)
                                                                                  tacgttacgtGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGATGCCGCAata        <  2:4645251/93‑1 (MQ=255)
                                                                                    cgttacgcGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGATGCCGCAATAtt      >  2:2264359/1‑93 (MQ=255)
                                                                                     gttacgcGACGGGGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGATGCCGCAATATTa     <  1:4300054/93‑1 (MQ=255)
                                                                                       tacgcGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGATGCCGCAATATTAcg   <  2:2073735/93‑1 (MQ=255)
                                                                                       tacgcGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGATGCCGCAATATTAcg   <  2:3388491/93‑1 (MQ=255)
                                                                                       tacgcGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGATGCCGCAATATTAcg   >  2:1844928/1‑93 (MQ=255)
                                                                                        acgtGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGATGCCGCAATATTAcgg  <  1:3095829/93‑1 (MQ=255)
                                                                                        acgcGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGATGCCGCAATATTAcgg  <  2:5002497/93‑1 (MQ=255)
                                                                                        acgcGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGATGCCGCAATATTAcgg  >  1:3896214/1‑93 (MQ=255)
                                                                                           |                                                                                         
AACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCTCGAAATGCGCTTCAAGCCGATGCCGCAATATTACGG  >  NC_000913/368856‑369036

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