Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 371,070 C→G 15.2% 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.2% 529.4 / 80.9 192A265G (GCC→GGC) mhpB2,3‑dihydroxyphenylpropionate 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (86/76);  new base G (14/15);  total (101/91)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.04e-01

Only 100 of 203 total aligned reads displayed.

CGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAAT  >  NC_000913/370971‑371169
                                                                                                   |                                                                                                   
cGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGc                                                                                                     >  2:587918/1‑100 (MQ=255)
cGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGc                                                                                                     >  1:127999/1‑100 (MQ=255)
cGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGc                                                                                                     >  1:2074171/1‑100 (MQ=255)
 gCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGcc                                                                                                    <  1:962096/100‑1 (MQ=255)
 gCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGcc                                                                                                    <  2:346217/100‑1 (MQ=255)
  ccgatctcGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCg                                                                                                   <  1:3893755/93‑1 (MQ=255)
  cTCATCCCGATTTGGGATAACCAGTTCATTACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATCCCGTCATTAACGAAGAGCTTTCCGCCATTGCCg                                                                                                   <  2:2514026/100‑1 (MQ=255)
  cTCATCCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCg                                                                                                   <  2:2072869/100‑1 (MQ=255)
  cTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCCGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCg                                                                                                   <  2:2715454/100‑1 (MQ=255)
  cTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCg                                                                                                   <  2:3380523/100‑1 (MQ=255)
  cTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCg                                                                                                   <  1:1755579/100‑1 (MQ=255)
    cAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGa                                                                                                 <  1:3471208/100‑1 (MQ=255)
     aaCCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGaa                                                                                                <  1:4082258/100‑1 (MQ=255)
      aCCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGaaa                                                                                               >  1:3919072/1‑100 (MQ=255)
       cccGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAg                                                                                              <  2:3537588/100‑1 (MQ=255)
       cccGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAg                                                                                              >  1:2015758/1‑100 (MQ=255)
       cccGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAg                                                                                              >  2:775006/1‑100 (MQ=255)
        ccGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGAAAgc                                                                                             <  1:2700956/100‑2 (MQ=255)
        ccGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGt                                                                                             >  1:258329/1‑100 (MQ=255)
        ccGATTTGCGATAACCAGTTCATGACTTTGCTGGAGCAGAGACTCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGACAATGGCGGAAAgn                                                                                             <  2:1845701/100‑2 (MQ=255)
         cGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGAAAgct                                                                                            <  1:633036/100‑3 (MQ=255)
         cGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCgagatcgg                                                                                            >  2:3893730/1‑93 (MQ=255)
         cGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTc                                                                                            <  1:429521/100‑1 (MQ=255)
           aTTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGAAAgctct                                                                                          <  2:173496/100‑5 (MQ=255)
           aTTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCg                                                                                           >  2:3421689/1‑99 (MQ=255)
            tttGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGac                                                                                         <  2:1973638/100‑1 (MQ=255)
             ttGGGATAACCAGGTCATGACTTTGCGGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGAAAgctcttc                                                                                        <  2:3741751/100‑7 (MQ=37)
               gggATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGAcaca                                                                                      >  1:903353/1‑100 (MQ=255)
               gggATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGAcaca                                                                                      >  2:4070594/1‑100 (MQ=255)
                     aCCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAAt                                                                                >  1:1550944/1‑100 (MQ=255)
                     aCCAGTTCATCACTTTGCTGCAGCATGCACGCATCCAGGAACTGTATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAATGTCTACACATg                                                                                    <  2:2192041/96‑1 (MQ=255)
                      ccAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATc                                                                               >  2:2886012/1‑100 (MQ=255)
                       cAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCa                                                                              >  1:110191/1‑100 (MQ=255)
                            cATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAAcc                                                                         >  2:1009686/1‑100 (MQ=255)
                            cATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAAcc                                                                         >  1:2288726/1‑100 (MQ=255)
                            cATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAAcc                                                                         >  1:2171666/1‑100 (MQ=255)
                              tGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTg                                                                       <  1:2536430/100‑1 (MQ=255)
                              tGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTg                                                                       <  2:2421418/100‑1 (MQ=255)
                              tGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTg                                                                       <  1:4120270/100‑1 (MQ=255)
                                 cTTTGCTGGAGCAGGGACGCATACAGGCACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGt                                                                    <  2:1426391/100‑1 (MQ=255)
                                   ttGCTGGAGCAGGGACGCATACAGGAAATGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAATCGACACATGAAATCAAAACCTGGGGcg                                                                  >  1:1467348/1‑100 (MQ=255)
                                       tGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTcgccgc                                                              >  2:1574910/1‑100 (MQ=255)
                                       tGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTcgccgc                                                              >  1:876576/1‑100 (MQ=255)
                                        ggAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCt                                                             >  1:2648095/1‑100 (MQ=255)
                                        ggAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCt                                                             >  1:732452/1‑100 (MQ=255)
                                        ggAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCt                                                             >  1:1020115/1‑100 (MQ=255)
                                        ggAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGACGGAAAGTCGACACATAAAATCACAACCAGGGTtgccgc                                                              >  1:1814221/1‑99 (MQ=255)
                                           gCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCtttt                                                          <  1:3568704/100‑1 (MQ=255)
                                            cAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTg                                                         <  2:3791411/100‑1 (MQ=255)
                                              gggACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGcc                                                       <  1:3155230/100‑1 (MQ=255)
                                                gACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGcc                                                       <  2:1524432/98‑1 (MQ=255)
                                                gACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGc                                                     >  2:3001332/1‑100 (MQ=255)
                                                gACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGc                                                     <  1:1714330/100‑1 (MQ=255)
                                                 aCGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAACTCAAAACCTGGGTCGCGGCTTTTGCCGcn                                                    <  2:1462730/100‑2 (MQ=255)
                                                  cGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTa                                                   <  1:51384/100‑1 (MQ=255)
                                                     aTACAGGAGCTGTATGCCGTCAGTAACGAAGAGCTTGCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTAtt                                                 <  2:4171973/99‑1 (MQ=255)
                                                     aTACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTAtt                                                 >  2:2464622/1‑99 (MQ=255)
                                                       aCAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCt                                              >  1:2363001/1‑100 (MQ=255)
                                                        cAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTg                                             <  2:2362996/100‑1 (MQ=255)
                                                        cAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTg                                             >  2:1304800/1‑100 (MQ=255)
                                                          ggAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCg                                           <  1:2246669/100‑1 (MQ=255)
                                                           gAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGt                                          <  1:119791/100‑1 (MQ=255)
                                                             aCTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGttt                                        <  1:3818955/100‑1 (MQ=255)
                                                                ggATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGc                                     <  1:280157/100‑1 (MQ=255)
                                                                  aTGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCag                                   >  2:3597483/1‑99 (MQ=255)
                                                                    gCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACt                                 >  1:2693930/1‑100 (MQ=255)
                                                                       gTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGc                              <  1:1249611/100‑1 (MQ=255)
                                                                         cAGTAACTAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATTAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTTCGTTTTGCAACTGGCGt                            <  2:836794/100‑1 (MQ=255)
                                                                          aGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGCAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTa                           <  2:1403399/100‑1 (MQ=255)
                                                                          aGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTa                           >  2:1103181/1‑100 (MQ=255)
                                                                          aGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTa                           >  1:2183579/1‑100 (MQ=255)
                                                                             aaCGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCg                        >  2:3648874/1‑100 (MQ=255)
                                                                               cGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGaa                      >  1:4054140/1‑100 (MQ=255)
                                                                                tgtcgacTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAg                     <  1:2878650/94‑1 (MQ=38)
                                                                                gAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAg                     <  1:4198112/100‑1 (MQ=255)
                                                                                 gtcgacTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAg                     >  2:583034/6‑99 (MQ=38)
                                                                                  agagCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAggg                   >  1:1669154/1‑100 (MQ=255)
                                                                                    gacTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCg                 >  2:300958/3‑100 (MQ=255)
                                                                                    agCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCg                 >  1:2458392/1‑100 (MQ=255)
                                                                                     acTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGt                >  2:2666447/2‑100 (MQ=255)
                                                                                      cTTTCCGGCGTTGCCGGAGAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGtt               <  2:1211915/100‑1 (MQ=255)
                                                                                      cTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGtt               >  1:1944825/1‑100 (MQ=255)
                                                                                      cTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGtt               >  2:2112971/1‑100 (MQ=255)
                                                                                      cTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGtt               >  2:102092/1‑100 (MQ=255)
                                                                                      cTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGtt               >  1:2561491/1‑100 (MQ=255)
                                                                                      cTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGtt               >  1:1688418/1‑100 (MQ=255)
                                                                                       tttCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGtta              <  2:514038/100‑1 (MQ=255)
                                                                                       tttCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGtta              <  2:3202167/100‑1 (MQ=255)
                                                                                        ttCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGttat             <  2:2140702/100‑1 (MQ=255)
                                                                                           cGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTAc          >  1:3017162/1‑100 (MQ=255)
                                                                                           cGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTAc          >  1:555989/1‑100 (MQ=255)
                                                                                            gCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTAcc         >  1:1646165/1‑100 (MQ=255)
                                                                                              cATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGc       <  2:1379232/100‑1 (MQ=255)
                                                                                              cATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGc       >  1:3630408/1‑100 (MQ=255)
                                                                                              cATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGc       <  2:473635/100‑1 (MQ=255)
                                                                                               aTTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGGCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGcc      >  2:2231211/1‑100 (MQ=255)
                                                                                                 tGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCa    >  1:1836833/1‑100 (MQ=255)
                                                                                                   ccGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAAt  >  1:3328171/1‑100 (MQ=255)
                                                                                                   ccGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAAt  >  2:4119648/1‑100 (MQ=255)
                                                                                                   ccGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAAt  >  1:710233/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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: