Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,669,800 T→C C186R (TGT→CGT)  hcaB → 2,3‑dihydroxy‑2,3‑dihydrophenylpropionate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,669,8000TC99.7% 875.0 / ‑5.3 282C186R (TGT→CGT) hcaB2,3‑dihydroxy‑2,3‑dihydrophenylpropionate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (122/159);  total (123/159)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.29e-01

Only 100 of 303 total aligned reads displayed.

TACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACGCCGGAGAAAATTGCC  >  NC_000913/2669710‑2669889
                                                                                          |                                                                                         
 acacCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGcgtg                                                                                        <  2:5291630/93‑1 (MQ=255)
    ccGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTa                                                                                     <  1:256411/93‑1 (MQ=255)
     cGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTAt                                                                                    <  2:2431764/93‑1 (MQ=255)
      gctaGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATg                                                                                   <  1:3827026/90‑1 (MQ=255)
       ccAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATgg                                                                                  >  2:7688534/1‑93 (MQ=255)
         aGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGcc                                                                                >  1:984435/1‑93 (MQ=255)
         aGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGcc                                                                                >  1:4043299/1‑93 (MQ=255)
         aGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGcc                                                                                <  2:6792610/93‑1 (MQ=255)
          gTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCa                                                                               >  1:5526112/1‑93 (MQ=255)
          gTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCa                                                                               <  1:3016090/93‑1 (MQ=255)
            aaaCATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGc                                                                             <  2:3748431/93‑1 (MQ=255)
             aaCATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCg                                                                            >  1:2127240/1‑93 (MQ=255)
              aCATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGa                                                                           >  2:5213825/1‑93 (MQ=255)
              aCATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGa                                                                           >  2:369229/1‑93 (MQ=255)
              aCATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTACGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCgc                                                                           >  2:2754354/1‑92 (MQ=255)
                aTGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGAcc                                                                         <  2:6476885/93‑1 (MQ=255)
                aTGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGAcc                                                                         <  1:2016728/93‑1 (MQ=255)
                aTGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGAcc                                                                         >  2:7580961/1‑93 (MQ=255)
                 tGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCt                                                                        >  1:4785433/1‑93 (MQ=255)
                 tGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCt                                                                        <  1:1589664/93‑1 (MQ=255)
                 tGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGCATGGCCAGCGACCt                                                                        >  1:1359685/1‑93 (MQ=255)
                  gCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTg                                                                       <  1:798115/93‑1 (MQ=255)
                    cGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTgcg                                                                     >  2:2135126/1‑93 (MQ=255)
                     gCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGTCCTgcgc                                                                    <  2:3122654/93‑1 (MQ=255)
                     gCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTgcgc                                                                    >  1:153733/1‑93 (MQ=255)
                       aaCCGGACTTATTCGCCAACTGGCTTGTGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTCTGGCCAGCGACCTGCGCgg                                                                  >  1:2620291/1‑93 (MQ=255)
                       aaCCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCgg                                                                  >  2:3354098/1‑93 (MQ=255)
                       aaCCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCgg                                                                  >  2:2821118/1‑93 (MQ=255)
                        aCCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGc                                                                 >  2:4830532/1‑93 (MQ=255)
                          cGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGccc                                                               <  2:7216485/93‑1 (MQ=255)
                            gACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCcac                                                             <  1:5752865/93‑1 (MQ=255)
                            gACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCcac                                                             <  1:1094156/93‑1 (MQ=255)
                              cTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAg                                                           <  1:7081585/93‑1 (MQ=255)
                               ttATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAgg                                                          <  2:2247350/93‑1 (MQ=255)
                                tATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGgc                                                         >  2:5972611/1‑93 (MQ=255)
                                tATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGgc                                                         >  2:7919649/1‑93 (MQ=255)
                                tATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGgc                                                         <  2:5929084/93‑1 (MQ=255)
                                  ttCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGAATGGCCAGCGACCTGCGCGGCCCACACgcgc                                                       >  2:6209423/1‑93 (MQ=255)
                                     gCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCg                                                    <  2:5710888/93‑1 (MQ=255)
                                      ccAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCgg                                                   <  2:5287585/93‑1 (MQ=255)
                                      ccAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCgg                                                   >  2:5045403/1‑93 (MQ=255)
                                        aaCTGGCTTATGAACTGGCACCGAAAGTGCGTGTGAATGGCGTCGGCCTGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGc                                                 <  1:3640645/93‑1 (MQ=255)
                                        aaCTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGc                                                 >  1:7675143/1‑93 (MQ=255)
                                         aCTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCa                                                >  2:7626168/1‑93 (MQ=255)
                                          cTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCaa                                               >  2:6555171/1‑93 (MQ=255)
                                           tGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCaaa                                              <  1:3947995/93‑1 (MQ=255)
                                              cTTATGGACTGGCACGGAAAGTGCGGGTGAATTGCGTCGGCCCGCGGGGGATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTg                                           <  2:3975258/93‑1 (MQ=255)
                                               ttATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGa                                          >  1:4374953/1‑93 (MQ=255)
                                                tATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGaa                                         >  1:896336/1‑93 (MQ=255)
                                                 aTGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGaaa                                        <  1:7777328/93‑1 (MQ=255)
                                                  tGAACTGTCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAAc                                       <  2:642845/93‑1 (MQ=255)
                                                  tGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAAc                                       >  1:5842897/1‑93 (MQ=255)
                                                   gAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAAcc                                      >  1:4631418/1‑93 (MQ=255)
                                                    aaCTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGAGCTCGGGCAAAGTGAAACCt                                     <  1:7224202/93‑1 (MQ=255)
                                                     aCTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTc                                    <  2:3900298/93‑1 (MQ=255)
                                                     aCTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTc                                    >  1:4285187/1‑93 (MQ=255)
                                                      cTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCg                                   >  1:6630253/1‑93 (MQ=255)
                                                      cTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCg                                   >  1:6857686/1‑93 (MQ=255)
                                                        ggCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGAt                                 <  1:3274690/93‑1 (MQ=255)
                                                        ggCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAAGCGCTCGGGCAAAGTGAAAACTCGAt                                 >  2:4760855/1‑93 (MQ=255)
                                                         gCACCGAAAGTGCGGGTGAAGGGCGTCGGCCCGCGTGGTAGGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATa                                <  2:1623651/93‑1 (MQ=255)
                                                          cACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATaa                               <  2:5484294/93‑1 (MQ=255)
                                                          cACCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATaa                               <  1:3438835/93‑1 (MQ=255)
                                                           aCCGAAAGTGCGGGTGAATGGCGTCGGCCTGCTTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCTCTCGGGCAAAGTGAAACCTCGATAAt                              <  1:4972768/93‑1 (MQ=255)
                                                           aCCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAAt                              <  1:5213646/93‑1 (MQ=255)
                                                           aCCGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAAt                              >  1:5167585/1‑93 (MQ=255)
                                                            ccGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATg                             >  2:4946971/1‑93 (MQ=255)
                                                            ccGAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATg                             >  1:4418933/1‑93 (MQ=255)
                                                              gAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCa                           <  1:2491644/93‑1 (MQ=255)
                                                              gAAAGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCa                           >  1:2894616/1‑93 (MQ=255)
                                                               aaaGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAg                          <  1:3351415/93‑1 (MQ=255)
                                                               aaaGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAg                          >  1:5126095/1‑93 (MQ=255)
                                                                aaGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGt                         >  1:5448632/1‑93 (MQ=255)
                                                                aaGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGt                         <  1:2992976/93‑1 (MQ=255)
                                                                aaGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGt                         >  1:6366419/1‑93 (MQ=255)
                                                                aaGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGt                         >  2:4419182/1‑93 (MQ=255)
                                                                aaGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGt                         <  1:2442762/93‑1 (MQ=255)
                                                                 aGTGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGtc                        <  2:7655821/93‑1 (MQ=255)
                                                                   tGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGtctc                      <  1:5247105/93‑1 (MQ=255)
                                                                   tGCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGtctc                      <  2:1503798/93‑1 (MQ=255)
                                                                   ggCGGGTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGtctc                      <  1:294233/92‑1 (MQ=255)
                                                                      gggTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGa                   <  2:7370393/93‑1 (MQ=255)
                                                                      gggTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGa                   >  2:1860193/1‑93 (MQ=255)
                                                                       ggTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCACAGTGGAACCTCGATAATGCAGTCTCTGAc                  >  1:2101992/1‑93 (MQ=255)
                                                                       ggTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTgtc                  <  1:7121210/93‑3 (MQ=255)
                                                                       ggTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGAc                  <  2:4173022/93‑1 (MQ=255)
                                                                       ggTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGAc                  >  2:1085474/1‑93 (MQ=255)
                                                                        gTGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACg                 >  2:3163980/1‑93 (MQ=255)
                                                                         tGAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACGc                >  1:7878629/1‑93 (MQ=255)
                                                                          gAATGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACGcc               >  2:1339580/1‑93 (MQ=255)
                                                                            aTGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACGCCgg             <  1:1733890/93‑1 (MQ=255)
                                                                             tGGCGTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACGCCGga            <  1:3298406/93‑1 (MQ=255)
                                                                                 gTCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACGCCGGAGaaa        >  2:6954769/1‑93 (MQ=255)
                                                                                  tCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAGGTGAAACCTCGATAATGCAGTCTCTGACGCCGGAGaaaa       <  1:2754355/93‑1 (MQ=255)
                                                                                  tCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACGCCGGAGaaaa       <  1:2304443/93‑1 (MQ=255)
                                                                                  tCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACGCCGGAGaaaa       >  2:3574410/1‑93 (MQ=255)
                                                                                  tCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACGCCGGAGaaaa       <  1:6312564/93‑1 (MQ=255)
                                                                                  tCGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACGCCGGAGaaaa       <  1:5996351/93‑1 (MQ=255)
                                                                                   cGGCCCGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACGCCGGAGAAAAt      <  1:649955/93‑1 (MQ=255)
                                                                                      cccGCGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTCACGCCGGAGAAAATTgc   >  2:366464/1‑93 (MQ=255)
                                                                                          |                                                                                         
TACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACGCCGGAGAAAATTGCC  >  NC_000913/2669710‑2669889

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