Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 368,911 G→A V359V (GTG→GTA)  mhpA → 3‑(3‑hydroxyphenyl)propionate hydroxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913368,9110GA83.9% 686.7 / 95.5 292V359V (GTG→GTA) mhpA3‑(3‑hydroxyphenyl)propionate hydroxylase
Reads supporting (aligned to +/- strand):  ref base G (21/26);  new base A (140/105);  total (161/131)
Fisher's exact test for biased strand distribution p-value = 1.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-05

Only 100 of 292 total aligned reads displayed.

AGCCCGCGATGCGCTGCTCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCT  >  NC_000913/368821‑369000
                                                                                          |                                                                                         
   ccGCGATGCGCTGCTCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGc                                                                                      >  1:1939952/1‑93 (MQ=255)
   ccGCGATGCGCTGCTCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGc                                                                                      <  2:1829380/93‑1 (MQ=255)
    cgcgATGCGCTGCTCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAAAGTACTggat                                                                                     >  1:6572792/1‑91 (MQ=255)
       gATGCGCTGCTCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTccg                                                                                  <  1:6436052/93‑1 (MQ=255)
        aTGCGCTGCTCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTccgc                                                                                 >  2:6389829/1‑93 (MQ=255)
        aTGCGCTGCTCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTccgc                                                                                 <  1:228852/93‑1 (MQ=255)
          gcgcTGCTCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTccgccg                                                                               >  1:5257426/1‑93 (MQ=255)
           cgcTGCTCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGa                                                                              >  2:1024398/1‑93 (MQ=255)
            gctgctCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGaa                                                                             >  1:4268883/1‑93 (MQ=255)
            gctgctCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGaa                                                                             >  1:3096194/1‑93 (MQ=255)
             ctgctCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGaaa                                                                            <  2:5651857/93‑1 (MQ=255)
               gctCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACg                                                                          >  1:3701979/1‑93 (MQ=255)
                ctCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGc                                                                         >  2:1929850/1‑93 (MQ=255)
                 tCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCt                                                                        >  1:6918938/1‑93 (MQ=255)
                 tCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCt                                                                        <  2:3299905/93‑1 (MQ=255)
                 tCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCt                                                                        >  1:361419/1‑93 (MQ=255)
                  cGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTg                                                                       <  2:1776202/93‑1 (MQ=255)
                  cGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTg                                                                       >  2:3114187/1‑93 (MQ=255)
                     tACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCa                                                                    >  1:4909274/1‑93 (MQ=255)
                      aCCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAg                                                                   >  1:3450572/1‑93 (MQ=255)
                       ccTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAgg                                                                  <  2:6495368/93‑1 (MQ=255)
                       ccTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAgg                                                                  >  2:319826/1‑93 (MQ=255)
                        cTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAggg                                                                 <  2:671969/93‑1 (MQ=255)
                         tATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAgggg                                                                >  1:2645752/1‑92 (MQ=255)
                         tATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAgggg                                                                >  1:4204535/1‑92 (MQ=255)
                          accaacaaGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGta                                                               >  2:3635893/3‑93 (MQ=255)
                          aTCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGta                                                               <  2:3902014/93‑1 (MQ=255)
                           tCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGtac                                                              <  2:2879039/93‑1 (MQ=255)
                           tCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGtac                                                              <  1:5516406/93‑1 (MQ=255)
                            caacaaGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGtacg                                                             >  2:5828350/1‑93 (MQ=255)
                              acaaGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGtacgtt                                                           <  2:1097752/93‑1 (MQ=255)
                              acaaGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGgacgtt                                                           >  1:1724990/1‑93 (MQ=255)
                               caaGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGtacgtta                                                          >  2:3014500/1‑93 (MQ=255)
                               caaGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGtacgtta                                                          <  1:4757142/93‑1 (MQ=255)
                               caaGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGtacgtta                                                          <  2:6110129/93‑1 (MQ=255)
                                 aGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGtacgttacg                                                        >  1:347561/1‑93 (MQ=255)
                                 aGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGtacgttacg                                                        >  2:3984670/1‑93 (MQ=255)
                                 aGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGaacgttacg                                                        >  2:521955/1‑93 (MQ=255)
                                  gAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGtacgttacgt                                                       <  1:1333845/93‑1 (MQ=255)
                                  gAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGgacgttacgt                                                       >  1:3918604/1‑93 (MQ=255)
                                  gAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCCCCGAAACGCTGGCAGGGtacgttacgt                                                       <  2:6418198/93‑1 (MQ=255)
                                   aaCGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGGACGTTACGTg                                                      >  2:7396345/1‑93 (MQ=255)
                                    aCGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGa                                                     <  1:2962176/93‑1 (MQ=255)
                                    aCGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGGACGTTACGTGa                                                     >  1:2703421/1‑93 (MQ=255)
                                     cGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGAc                                                    <  2:4801901/93‑1 (MQ=255)
                                     cGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGAc                                                    <  2:5337393/93‑1 (MQ=255)
                                       tCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACgg                                                  >  2:4161484/1‑93 (MQ=255)
                                         gcgATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCg                                                <  2:3414542/93‑1 (MQ=255)
                                          cgATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGt                                               <  1:211626/93‑1 (MQ=255)
                                           gATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGGACGTTACGTGACGGCGtt                                              >  1:6182282/1‑93 (MQ=255)
                                            aTCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGttt                                             >  1:6284699/1‑93 (MQ=255)
                                            aTCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGttt                                             >  2:4909932/1‑93 (MQ=255)
                                            aTCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGttt                                             >  2:7007312/1‑93 (MQ=255)
                                             tCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTc                                            <  2:6605924/93‑1 (MQ=255)
                                              cACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTcc                                           <  2:2803/93‑1 (MQ=255)
                                               aCGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCt                                          <  2:63083/93‑1 (MQ=255)
                                               aCGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCt                                          >  2:3701624/1‑93 (MQ=255)
                                                cGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTg                                         <  2:551726/93‑1 (MQ=255)
                                                  ccAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGc                                       >  2:2105154/1‑93 (MQ=255)
                                                  ccAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGGACGTTACGTGACGGCGTTTCCTGGc                                       >  1:7066330/1‑93 (MQ=255)
                                                   cAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCt                                      <  2:4507226/93‑1 (MQ=255)
                                                   cAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCt                                      >  1:1040860/1‑93 (MQ=255)
                                                    aaaGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTg                                     >  2:627703/1‑93 (MQ=255)
                                                    aaaGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGGACGTTACGTGACGGCGTTTCCTGGCTg                                     >  2:2720442/1‑93 (MQ=255)
                                                    aaaGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGGACGTTACGTGACGGCGTTTCCTGGCTg                                     >  1:25342/1‑93 (MQ=255)
                                                     aaGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGGACGTTACGTGACGGCGTTTCCTGGCTGt                                    >  1:2849492/1‑93 (MQ=255)
                                                         gatgatTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGaa                                <  1:2208312/93‑1 (MQ=255)
                                                         gatgatTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGaa                                <  1:5795248/93‑1 (MQ=255)
                                                         gatgatTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGaa                                >  2:6827715/1‑93 (MQ=255)
                                                         gatgatTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGGACGTTACGTGACGGCGTTTCCTGGCTGTTGaa                                >  1:1633696/1‑93 (MQ=255)
                                                         gatgatTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGGACGTTACGTGACGGCGTTTCCTGGCTGTTGaa                                >  1:4094390/1‑93 (MQ=255)
                                                          atgatTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAAt                               >  2:1656065/1‑93 (MQ=255)
                                                           tgatTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAAtt                              >  2:1597/1‑93 (MQ=255)
                                                            gatTGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTa                             <  2:5393604/93‑1 (MQ=255)
                                                              ttGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTAtt                           >  1:823611/1‑92 (MQ=255)
                                                               tGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCt                          >  2:7233925/1‑93 (MQ=255)
                                                               tGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCt                          <  1:2340223/93‑1 (MQ=255)
                                                               tGACCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCt                          <  1:534146/93‑1 (MQ=255)
                                                                 aCCTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTgc                        >  2:836505/1‑93 (MQ=255)
                                                                  ccTGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTgcc                       <  1:5092769/93‑1 (MQ=255)
                                                                    tGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTgccgc                     >  1:1996906/1‑93 (MQ=255)
                                                                    tGTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTgccgc                     >  2:4822833/1‑93 (MQ=255)
                                                                     gTCCGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGGACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTgccgcc                    >  1:5631258/1‑93 (MQ=255)
                                                                       ccGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAg                  >  2:2593561/1‑93 (MQ=255)
                                                                       ccGTGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAg                  <  2:5939045/93‑1 (MQ=255)
                                                                          tGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTaa               <  2:4354852/93‑1 (MQ=255)
                                                                          tGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTaa               <  2:2395358/93‑1 (MQ=255)
                                                                          tGACGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTaa               >  2:2834217/1‑93 (MQ=255)
                                                                            aCGGCGGGCAACGTACTTGCTCCGCCGAAACGCTGGCAGGGGACGTTACGTGACGGCGTTTTCTTGCTGTTGAATTATCTGCCGCCAGTaaaa             >  1:3334513/1‑93 (MQ=255)
                                                                            aCGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTaaaa             <  2:163285/93‑1 (MQ=255)
                                                                            aCGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTaaaa             <  2:6564057/93‑1 (MQ=255)
                                                                            aCGGCGGGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTaaaa             <  1:4611382/93‑1 (MQ=255)
                                                                               gcggGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGc          <  2:142394/93‑1 (MQ=255)
                                                                                cggGCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCt         >  1:973056/1‑93 (MQ=255)
                                                                                 gggCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTa        <  2:7093937/93‑1 (MQ=255)
                                                                                  ggCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTAc       >  1:447504/1‑93 (MQ=255)
                                                                                   gCAACGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACt      >  2:1563168/1‑93 (MQ=255)
                                                                                     aaCGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTc    >  1:2786674/1‑93 (MQ=255)
                                                                                     aaCGTACTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTc    >  1:1139024/1‑93 (MQ=255)
                                                                                      aCGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTcc   <  2:2178516/93‑1 (MQ=255)
                                                                                          |                                                                                         
AGCCCGCGATGCGCTGCTCGATACCTATCAACAAGAACGTCGCGATCACGCCAAAGCGATGATTGACCTGTCCGTGACGGCGGGCAACGTGCTGGCTCCGCCGAAACGCTGGCAGGGTACGTTACGTGACGGCGTTTCCTGGCTGTTGAATTATCTGCCGCCAGTAAAACGCTACTTCCT  >  NC_000913/368821‑369000

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