Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 371,070 C→G 15.1% A265G (GCC→GGC)  mhpB → 2,3‑dihydroxyphenylpropionate 1,2‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913371,0700CG15.1% 459.2 / 67.2 165A265G (GCC→GGC) mhpB2,3‑dihydroxyphenylpropionate 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (78/62);  new base G (14/11);  total (92/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-01

Only 100 of 171 total aligned reads displayed.

CGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAAT  >  NC_000913/370971‑371169
                                                                                                   |                                                                                                   
cGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGc                                                                                                     <  1:2224197/100‑1 (MQ=255)
 gCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGcc                                                                                                    >  2:2089386/1‑100 (MQ=255)
     aaCCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGaa                                                                                                >  1:411754/1‑100 (MQ=255)
      aCCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGaaa                                                                                               >  1:465261/1‑100 (MQ=255)
       cccGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAg                                                                                              >  2:2597966/1‑100 (MQ=255)
        ccGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGt                                                                                             >  1:1579125/1‑100 (MQ=255)
        ccGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAAGGACGCATACAGGAACTGGGTGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGAAAgc                                                                                             <  2:1988009/100‑2 (MQ=255)
         cGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGAAAgct                                                                                            <  1:2789515/100‑3 (MQ=255)
         cGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTc                                                                                            >  1:3445975/1‑100 (MQ=255)
         cGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTc                                                                                            >  2:3600999/1‑100 (MQ=255)
          gATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGAAAgctc                                                                                           <  1:1417180/100‑4 (MQ=255)
          gATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGAAAgctc                                                                                           <  1:68625/100‑4 (MQ=255)
           aTTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGAAAgctct                                                                                          <  1:3390100/100‑5 (MQ=255)
            tttGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGan                                                                                         <  1:2507950/100‑2 (MQ=255)
             ttGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGAAAgctcttc                                                                                        <  1:2268103/100‑7 (MQ=38)
             ttGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGaca                                                                                        <  2:2592400/100‑1 (MQ=255)
              agggATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGacag                                                                                       >  2:1089860/2‑99 (MQ=255)
               gggATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGAAAgctcttcgt                                                                                      <  2:770144/100‑9 (MQ=37)
               gggATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGAAAgctcttcgt                                                                                      <  2:2725612/100‑9 (MQ=37)
               gggATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGAcaca                                                                                      <  1:1298567/100‑1 (MQ=255)
                     aCCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAAt                                                                                >  2:1996275/1‑100 (MQ=255)
                      ccAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATc                                                                               <  1:2521628/100‑1 (MQ=255)
                      ccAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATc                                                                               <  1:2945858/100‑1 (MQ=255)
                       cAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCa                                                                              <  2:1781080/100‑1 (MQ=255)
                       cAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCa                                                                              >  1:157537/1‑100 (MQ=255)
                         gTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCaaa                                                                            <  2:509641/100‑1 (MQ=255)
                         gTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCaaa                                                                            >  2:1324603/1‑100 (MQ=255)
                         gTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACAAATGAAATCaaa                                                                            <  2:2308505/100‑1 (MQ=255)
                          ttCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCaaaa                                                                           >  2:3318495/1‑100 (MQ=255)
                           tCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAAc                                                                          <  1:3503149/100‑1 (MQ=255)
                              tGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTg                                                                       >  2:2051588/1‑100 (MQ=255)
                              tGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTg                                                                       <  2:626994/100‑1 (MQ=255)
                               gACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTgg                                                                      <  1:937585/100‑1 (MQ=255)
                                aCTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTggg                                                                     <  1:3106809/100‑1 (MQ=255)
                                  tttGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTggg                                                                     <  2:2510685/98‑1 (MQ=255)
                                     gCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTcgcc                                                                >  2:219160/1‑100 (MQ=255)
                                      cTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTcgccg                                                               >  1:1978329/1‑100 (MQ=255)
                                          aGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCttt                                                           <  1:56183/100‑1 (MQ=255)
                                           gCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCtttt                                                          >  1:1959594/1‑100 (MQ=255)
                                            cAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCttt                                                           >  2:2206105/1‑98 (MQ=255)
                                            cAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTg                                                         <  2:1667837/100‑1 (MQ=255)
                                             aGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGc                                                        >  2:1330456/1‑100 (MQ=255)
                                              gggACGCATACAGGAACTGGATGCCGTCAGTAAGGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACATGGGTCGCCGCTCTTg                                                         >  2:3015335/1‑98 (MQ=255)
                                              gggACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGACATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGGCGCTTTTGcc                                                       >  2:3294339/1‑100 (MQ=255)
                                               ggACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCg                                                      >  1:1609077/1‑100 (MQ=255)
                                                  cGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTa                                                   >  2:1070379/1‑100 (MQ=255)
                                                    cATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTAtt                                                 >  2:2949046/1‑100 (MQ=255)
                                                    cATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTAtt                                                 >  1:1096524/1‑100 (MQ=255)
                                                    cATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTAtt                                                 >  1:2657114/1‑100 (MQ=255)
                                                      tACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTc                                               >  2:364685/1‑100 (MQ=255)
                                                      tACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTc                                               >  1:1950826/1‑100 (MQ=255)
                                                      tACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTc                                               <  1:3140445/100‑1 (MQ=255)
                                                      tACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTc                                               <  1:281357/100‑1 (MQ=255)
                                                        cAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTg                                             >  2:507961/1‑100 (MQ=255)
                                                          ggAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCg                                           >  1:2021577/1‑100 (MQ=255)
                                                          ggAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCg                                           >  1:2547296/1‑100 (MQ=255)
                                                           gAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGt                                          >  1:606500/1‑100 (MQ=255)
                                                           gAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGt                                          >  2:2363723/1‑100 (MQ=255)
                                                            aaCTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGtt                                         <  2:852862/100‑1 (MQ=255)
                                                             aCTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGttt                                        >  1:3342834/1‑100 (MQ=255)
                                                              cTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTg                                       <  1:1548944/100‑1 (MQ=255)
                                                               tGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTgg                                      >  1:85360/1‑100 (MQ=255)
                                                                ggATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGc                                     <  2:2843223/100‑1 (MQ=255)
                                                                ggATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGc                                     <  2:320580/100‑1 (MQ=255)
                                                                ggATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGc                                     >  1:1363628/1‑100 (MQ=255)
                                                                ggATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGc                                     >  1:871136/1‑100 (MQ=255)
                                                                ggATGCCGTCAGTAACGAAGAGCTTTACGCCATTGCCGGTAAGACGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGc                                     >  2:3339284/1‑100 (MQ=255)
                                                                  aTGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCaa                                   >  1:3424049/1‑100 (MQ=255)
                                                                  aTGCCGTCAGTAACGAAGAGCTTTCCGCCAGTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTgg                                      <  2:2333986/97‑1 (MQ=255)
                                                                  aTGCCGTCAGTAACGAAGAGCTCTCCGTCATTGCCGGAACGTCGACACATGAAATCACAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCa                                    <  2:282723/99‑1 (MQ=255)
                                                                        tCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCg                             >  1:2763922/1‑100 (MQ=255)
                                                                           gTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAg                          >  1:1647309/1‑100 (MQ=255)
                                                                           gTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGATTGGCAACTGGCGTAg                          >  2:1717227/1‑100 (MQ=255)
                                                                           gTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCAGCTATTTCTGCGTTTGGCAACTGGCGTa                           >  2:784897/1‑99 (MQ=255)
                                                                            tAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGc                         >  1:1675036/1‑100 (MQ=255)
                                                                            tAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGc                         >  1:217250/1‑100 (MQ=255)
                                                                               cGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGaa                      >  1:987887/1‑100 (MQ=255)
                                                                                gAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAg                     >  2:3433818/1‑100 (MQ=255)
                                                                                gAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAg                     <  1:571675/100‑1 (MQ=255)
                                                                                  tcgacTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAgg                    >  1:1730167/5‑99 (MQ=39)
                                                                                   cgacTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGc                  >  2:1645057/4‑100 (MQ=255)
                                                                                    gacTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAggg                   >  1:499916/3‑98 (MQ=255)
                                                                                    gacTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCg                 >  2:1500670/3‑100 (MQ=255)
                                                                                    gacTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCg                 >  1:3613449/3‑100 (MQ=255)
                                                                                      cTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGtt               >  1:2428539/1‑100 (MQ=255)
                                                                                       tttCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGtta              <  1:105750/100‑1 (MQ=255)
                                                                                        ttCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGttat             <  2:957041/100‑1 (MQ=255)
                                                                                          ccGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGttatta           <  1:630262/100‑1 (MQ=255)
                                                                                            gCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTAcc         >  2:3420007/1‑100 (MQ=255)
                                                                                              cATGGCCGGAAAGTCGACACATGAAGTCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGc       <  2:3040334/100‑1 (MQ=255)
                                                                                               aTTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGcc      <  1:3574489/100‑1 (MQ=255)
                                                                                                ttGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGccc     <  2:3125195/100‑1 (MQ=255)
                                                                                                ttGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGccc     <  1:609721/100‑1 (MQ=255)
                                                                                                ttGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGccc     <  1:307977/100‑1 (MQ=255)
                                                                                                ttGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGccc     >  1:3007768/1‑100 (MQ=255)
                                                                                                ttGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGcc      <  2:1131863/99‑1 (MQ=255)
                                                                                                  gCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCaa   >  2:2850432/1‑100 (MQ=255)
                                                                                                  gCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCAAAGGGCGTTATTACCGCCCaa   >  1:1022461/1‑100 (MQ=255)
                                                                                                   ccGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAAt  >  1:830786/1‑100 (MQ=255)
                                                                                                   ccGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAAt  >  1:2054730/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAAT  >  NC_000913/370971‑371169

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