Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,671,759 C→A 35.7% R179R (CGG→AGG)  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,671,7590CA35.7% 220.9 / 190.5 189R179R (CGG→AGG) hcaB2,3‑dihydroxy‑2,3‑dihydrophenylpropionate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (72/49);  new base A (36/32);  total (108/81)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 192 total aligned reads displayed.

GCCTGGTATCCTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACGCCGG  >  NC_000913/2671660‑2671858
                                                                                                   |                                                                                                   
 ccTGGTATCCTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCg                                                                                                    >  2:2018722/1‑100 (MQ=255)
  cTGGTATCCTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCgg                                                                                                   >  2:1560972/1‑100 (MQ=255)
   tGGTATCCTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAggg                                                                                                  <  2:2947073/100‑1 (MQ=255)
    ggTATCCTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTCTTCGCCAACTGGCCTATGAACTGGCACCGAAAGTGCGGGt                                                                                                 >  1:3077131/1‑100 (MQ=255)
         ccTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATg                                                                                            <  1:4229463/100‑1 (MQ=255)
         ccTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATg                                                                                            >  2:3088986/1‑100 (MQ=255)
          cTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATgg                                                                                           <  2:3048332/100‑1 (MQ=255)
          cTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATgg                                                                                           >  2:2121345/1‑100 (MQ=255)
              cggcggTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTc                                                                                       >  2:4131374/1‑100 (MQ=255)
                 cggTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGc                                                                                    <  1:1037197/100‑1 (MQ=255)
                  ggTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGcc                                                                                   >  2:439069/1‑100 (MQ=255)
                  gatctcccGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGcc                                                                                   <  1:1979126/95‑1 (MQ=255)
                   gTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGccc                                                                                  >  1:1413208/1‑100 (MQ=255)
                   gTGGCCCGCTGTACACCGCCAGAAAACATGCCGCAACCGGACTTATTCGCCAACCGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGccc                                                                                  <  2:3916056/100‑1 (MQ=255)
                          gCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGt                                                                           >  2:2755462/1‑100 (MQ=255)
                          gCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGt                                                                           <  2:2933596/100‑1 (MQ=255)
                          gCTGTACACCGCCAGTAAACATGCCGCAACCGGACTGATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGt                                                                           <  2:2028445/100‑1 (MQ=255)
                           cTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTa                                                                          <  2:849015/100‑1 (MQ=255)
                            tGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTAt                                                                         >  2:3516801/1‑100 (MQ=255)
                              tACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATgg                                                                       >  1:1524997/1‑100 (MQ=255)
                              tACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATgg                                                                       <  1:2900829/100‑1 (MQ=255)
                              tACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATgg                                                                       <  2:3970804/100‑1 (MQ=255)
                              tACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATgg                                                                       >  2:3162068/1‑100 (MQ=255)
                               acacCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGc                                                                      <  2:253387/100‑1 (MQ=255)
                                cacCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGcc                                                                     >  1:2526700/1‑100 (MQ=255)
                                 acCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCa                                                                    >  2:2050662/1‑100 (MQ=255)
                                 acCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCCTACGAACTGCCACCGAACGTGCGCGTGACCGGCGTCGGCCCGTGGCGTATGGc                                                                      >  1:1016807/1‑98 (MQ=255)
                                  ccGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAg                                                                   >  2:3728609/1‑100 (MQ=255)
                                  ccGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAg                                                                   >  1:3822597/1‑100 (MQ=255)
                                   cGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAg                                                                   <  1:3937729/99‑1 (MQ=255)
                                   cGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGc                                                                  <  2:3200617/100‑1 (MQ=255)
                                   cGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGc                                                                  <  1:1827860/100‑1 (MQ=255)
                                    gCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCg                                                                 >  1:141721/1‑100 (MQ=255)
                                      cAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATCGCGACGGCCCGAGTGGTCTGGCCAGCGAc                                                               >  1:3452807/1‑100 (MQ=255)
                                       aGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGAcc                                                              >  2:1439979/1‑100 (MQ=255)
                                        gTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCt                                                             >  1:3774012/1‑100 (MQ=255)
                                         tAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTg                                                            <  2:3584061/100‑1 (MQ=255)
                                         tAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTg                                                            >  1:3915791/1‑100 (MQ=255)
                                          aaaCATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTgc                                                           >  1:2090888/1‑100 (MQ=255)
                                             cATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGcgcg                                                        >  2:3218716/1‑100 (MQ=255)
                                              aggCCGCAACCGGACTGATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCgg                                                       <  2:636583/98‑1 (MQ=255)
                                                 ccGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGccc                                                    <  2:1924814/100‑1 (MQ=255)
                                                 ccGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGccc                                                    <  2:1513646/100‑1 (MQ=255)
                                                   gCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCcac                                                  >  2:2632979/1‑100 (MQ=255)
                                                    cAACGGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCcaca                                                 <  1:2512435/100‑1 (MQ=255)
                                                      aCCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTACGCgg                                                       >  1:499446/1‑92 (MQ=255)
                                                            cTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGTGACCTGCGCGGCCCACAGGCGCTCg                                         >  2:3886851/1‑100 (MQ=255)
                                                            cTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCg                                         >  2:1249319/1‑100 (MQ=255)
                                                            cTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCg                                         <  1:923600/100‑1 (MQ=255)
                                                            cTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCg                                         >  2:3773344/1‑100 (MQ=255)
                                                             ttATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCgg                                        <  1:3925698/100‑1 (MQ=255)
                                                             ttATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCgg                                        <  2:377890/100‑1 (MQ=255)
                                                             ttATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCgg                                        <  1:1088666/100‑1 (MQ=255)
                                                             ttATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCgg                                        >  2:197634/1‑100 (MQ=255)
                                                              tATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCggg                                       >  2:132711/1‑100 (MQ=255)
                                                              tATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCggg                                       >  1:2575316/1‑100 (MQ=255)
                                                              tATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCggg                                       >  2:1814356/1‑100 (MQ=255)
                                                               aTTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGc                                      >  1:2702396/1‑100 (MQ=255)
                                                               aTTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGc                                      >  1:3976133/1‑100 (MQ=255)
                                                               aTTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGc                                      <  1:3240131/100‑1 (MQ=255)
                                                                 tCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCaa                                    >  1:3035609/1‑100 (MQ=255)
                                                                    ccAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGt                                 >  1:4165932/1‑100 (MQ=255)
                                                                    ccAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGt                                 <  2:688107/100‑1 (MQ=255)
                                                                     cAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTg                                <  1:2122738/100‑1 (MQ=255)
                                                                      aaCTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGa                               >  1:1617024/1‑100 (MQ=255)
                                                                        cTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAACGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGaaa                             >  1:1221070/1‑100 (MQ=255)
                                                                          ggCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAAcc                           >  2:2191628/1‑100 (MQ=255)
                                                                            cTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTc                         >  2:801176/1‑100 (MQ=255)
                                                                             ttATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCg                        >  1:2369633/1‑100 (MQ=255)
                                                                               aTGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGAt                      <  2:10556/100‑1 (MQ=255)
                                                                               aTGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCCGCCCGTGTTGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGAt                      <  1:984279/100‑1 (MQ=255)
                                                                                 gAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATaa                    <  1:1062761/100‑1 (MQ=255)
                                                                                   aCTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTAGGGCAAAGAGAAACCTAGATAATg                  >  2:3564602/1‑100 (MQ=255)
                                                                                    cTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGc                 <  1:1397263/100‑1 (MQ=255)
                                                                                    cTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGc                 <  2:577667/100‑1 (MQ=255)
                                                                                    cTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGc                 <  1:3807259/100‑1 (MQ=255)
                                                                                    cTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGc                 <  2:3774029/100‑1 (MQ=255)
                                                                                    cTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGc                 >  2:3466989/1‑100 (MQ=255)
                                                                                     tGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCAATAATGCa                >  1:3625427/1‑100 (MQ=255)
                                                                                     tGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCa                >  2:1040979/1‑100 (MQ=255)
                                                                                      ggCACCGAGAGTGAGGTGGAATGGCGTCGGCCCGGGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAg               <  1:171856/100‑1 (MQ=255)
                                                                                      ggCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAg               <  2:3428703/100‑1 (MQ=255)
                                                                                      ggCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAg               >  1:4102460/1‑100 (MQ=255)
                                                                                      ggCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAg               <  2:667005/100‑1 (MQ=255)
                                                                                        cACCGAAAGTGCGTGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCAGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGtc             >  1:3310668/1‑100 (MQ=255)
                                                                                        cACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGtc             >  1:2161158/1‑100 (MQ=255)
                                                                                         aCCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGtct            >  2:77811/1‑100 (MQ=255)
                                                                                          ccGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGtctc           <  2:1525001/100‑1 (MQ=255)
                                                                                          ccGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGtctc           <  2:1366999/100‑1 (MQ=255)
                                                                                           cGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTctct          >  2:1210566/1‑100 (MQ=255)
                                                                                            gAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATCATGCAGTCTCTg         >  1:1876328/1‑100 (MQ=255)
                                                                                            gAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCCGGACAAGTGAAACCTCGATAATGCAGTCTCTg         >  1:1109937/1‑100 (MQ=255)
                                                                                            gAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCAGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTg         >  2:3364610/1‑100 (MQ=255)
                                                                                              aaGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGTCAGCGACCCGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGAc       <  2:4069233/100‑1 (MQ=255)
                                                                                              aaGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGAc       >  1:602658/1‑100 (MQ=255)
                                                                                              aaGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGAc       <  2:3168255/100‑1 (MQ=255)
                                                                                                gTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGTGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACGc     >  1:3487444/1‑100 (MQ=255)
                                                                                                  gCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACGcc    <  1:822735/99‑1 (MQ=255)
                                                                                                  gCGGGTGAATGGCGTCGCCCCGTGTGGAATGGCCAGCGACCTGCTCGGCACACAGGCGCTCGGGCAAAGTGAATCCTCGATAATGCAGtc             <  1:2409207/90‑1 (MQ=255)
                                                                                                   cGGGTGAATGGCGTCGGCCCGTGTGGTATGGACAGCGACCTGCGCGGCCCACCGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACGCCgg  >  1:377417/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GCCTGGTATCCTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACGCCGG  >  NC_000913/2671660‑2671858

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