Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 371,064 C→G 16.4% A263G (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,0640CG16.4% 529.9 / 87.6 200A263G (GCC→GGC) mhpB2,3‑dihydroxyphenylpropionate 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (95/70);  new base G (19/14);  total (116/84)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

Only 100 of 206 total aligned reads displayed.

TGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCG  >  NC_000913/370965‑371163
                                                                                                   |                                                                                                   
tGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGc                                                                                                     <  1:1988694/100‑1 (MQ=255)
tGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGc                                                                                                     >  2:205696/1‑100 (MQ=255)
    tCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCAtt                                                                                                 >  2:2426944/1‑100 (MQ=255)
     ccGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTg                                                                                                >  1:3908430/1‑100 (MQ=255)
     ccGCTCAACCCGATTTGGGATAACCAGTGCATGACTTTGCTGGAGCAGGGACGCAGACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTg                                                                                                <  2:1627124/100‑1 (MQ=255)
      cGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGc                                                                                               >  2:2264981/1‑100 (MQ=255)
          cAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCgg                                                                                            >  1:642307/1‑99 (MQ=255)
          cAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCg                                                                                             <  2:1562191/98‑1 (MQ=255)
            aCCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGCACTGGATGCCGTCAGTAACGAAGAGCTTTCCGACATTGCCGGaaa                                                                                         >  2:176257/1‑100 (MQ=255)
            aCCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGaaa                                                                                         >  1:692516/1‑100 (MQ=255)
            aCCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGaaa                                                                                         >  1:2127398/1‑100 (MQ=255)
            aCCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGaaa                                                                                         >  2:1655005/1‑100 (MQ=255)
            aCCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGaaa                                                                                         >  1:2606174/1‑100 (MQ=255)
             cccGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAg                                                                                        <  2:1928085/100‑1 (MQ=255)
              ccGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGt                                                                                       >  2:3217415/1‑100 (MQ=255)
               cGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGAAAgct                                                                                      <  2:1175774/100‑3 (MQ=255)
               cGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGAAAgct                                                                                      <  2:1284247/100‑3 (MQ=255)
               cGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGAAAgct                                                                                      <  2:810373/100‑3 (MQ=255)
               cGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTc                                                                                      <  1:133847/100‑1 (MQ=255)
               cGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGCAACGAAGAGCTTTCCGGCAATGGCGGAAAgct                                                                                      <  1:1292954/100‑3 (MQ=255)
                gATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGTACGCATACAGGAACTGGATTCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGAAAgctc                                                                                     <  1:2391119/100‑4 (MQ=255)
                gATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCg                                                                                     <  2:2611127/100‑1 (MQ=255)
                  tttGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGAAAgctctt                                                                                   <  2:1715824/100‑6 (MQ=38)
                  tttGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGAAAgctctt                                                                                   <  2:329097/100‑6 (MQ=38)
                  tttGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGac                                                                                   >  1:2934203/1‑100 (MQ=255)
                   ttGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGAAAgctcttc                                                                                  <  2:3971518/100‑7 (MQ=38)
                   ttGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGAAAgctctt                                                                                   <  1:2060126/99‑6 (MQ=38)
                   ttGGGATAACCAGTTCATGACTTTGCTGGAGCAAGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGaca                                                                                  <  1:1772383/100‑1 (MQ=255)
                     gggATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGAcaca                                                                                <  1:1165860/100‑1 (MQ=255)
                      ggATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACAt                                                                               >  2:350014/1‑100 (MQ=255)
                      ggATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAACAGCTTTCCGCCATTGCCGGAAAGTCGACACAt                                                                               >  1:1583250/1‑100 (MQ=255)
                       gATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCTGAAAgctcttcgt                                                                                <  1:3655445/98‑9 (MQ=37)
                       gATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGAcaca                                                                                >  1:1425573/1‑98 (MQ=255)
                       gATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATg                                                                              <  1:4044876/100‑1 (MQ=255)
                        aTAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGa                                                                             <  2:33142/100‑1 (MQ=255)
                           aCCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAAt                                                                          >  2:1134810/1‑100 (MQ=255)
                            ccAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATc                                                                         >  2:362810/1‑100 (MQ=255)
                             cAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATc                                                                         >  1:3022297/1‑99 (MQ=255)
                              aGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCaa                                                                       >  1:229507/1‑100 (MQ=255)
                              aGTTCATGACTTTGCTGGAGCAGGGACGAATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCTGCCATTGCCGGAAAGTCGACACATGAAAt                                                                          <  1:1454190/97‑1 (MQ=255)
                                 tCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAAc                                                                    <  1:1042636/100‑1 (MQ=255)
                                  cATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAAcc                                                                   >  2:2492926/1‑100 (MQ=255)
                                    tGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTg                                                                 >  2:3732959/1‑100 (MQ=255)
                                        tttGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTc                                                             <  2:3554545/100‑1 (MQ=255)
                                          tGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTcgc                                                           >  2:1248294/1‑100 (MQ=255)
                                           gCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACTCATGAAATCCAAACCTGGGTcgcc                                                          <  1:940990/100‑1 (MQ=255)
                                           gCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTcgcc                                                          >  1:3173262/1‑100 (MQ=255)
                                             tGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTcgccgc                                                        >  1:4149074/1‑100 (MQ=255)
                                             tGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTcgccgc                                                        >  2:3466096/1‑100 (MQ=255)
                                                 gCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCtttt                                                    >  1:3849548/1‑100 (MQ=255)
                                                   aGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGc                                                  <  2:4361621/100‑1 (MQ=255)
                                                     ggACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGcc                                                 <  1:2614738/99‑1 (MQ=255)
                                                      gACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGc                                               <  2:13406/100‑1 (MQ=255)
                                                          cATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTAtt                                           >  1:720202/1‑100 (MQ=255)
                                                           aTACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTAttt                                          >  2:1055815/1‑100 (MQ=255)
                                                             aCAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCt                                        >  2:4334534/1‑100 (MQ=255)
                                                              cAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTg                                       <  1:808465/100‑1 (MQ=255)
                                                               aGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGc                                      <  2:2581467/100‑1 (MQ=255)
                                                               aGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGc                                      >  2:547245/1‑100 (MQ=255)
                                                                ggAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCg                                     <  1:2430697/100‑1 (MQ=255)
                                                                 gAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGt                                    <  2:418370/100‑1 (MQ=255)
                                                                    cTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTg                                 <  2:2447162/100‑1 (MQ=255)
                                                                    cTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTg                                 <  1:3309241/100‑1 (MQ=255)
                                                                       gATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAATACCTGGGGCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCa                              <  2:851488/100‑1 (MQ=255)
                                                                        aTGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCaa                             >  2:2749498/1‑100 (MQ=255)
                                                                         tGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAAc                            >  1:785017/1‑100 (MQ=255)
                                                                         tGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAAc                            <  1:4291153/100‑1 (MQ=255)
                                                                           ccGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTg                          >  2:1225934/1‑100 (MQ=255)
                                                                           ccGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTg                          >  1:3983776/1‑100 (MQ=255)
                                                                             gTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGc                        <  2:686792/100‑1 (MQ=255)
                                                                              tCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCg                       <  1:7536/100‑1 (MQ=255)
                                                                               cAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGt                      >  1:1404368/1‑100 (MQ=255)
                                                                                aGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTa                     <  2:412207/100‑1 (MQ=255)
                                                                                  tAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGc                   <  1:1105795/100‑1 (MQ=255)
                                                                                   aaCGAAGAGCTTGCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCg                  <  1:3487274/100‑1 (MQ=255)
                                                                                    aCGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGTCGCTTTTGCCGCTATTTCTGCGTTTGGCTACTGGCGTAGcca                 <  2:2980923/100‑3 (MQ=255)
                                                                                     cGAAGAGCTTTCCGCCAGTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGaa                <  1:1544332/100‑1 (MQ=255)
                                                                                      tgtcgacTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAg               >  2:3244367/7‑100 (MQ=38)
                                                                                       gtcgacTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAgg              >  2:640194/6‑100 (MQ=38)
                                                                                          gacTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCg           >  2:837562/3‑100 (MQ=255)
                                                                                          gacTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCg           >  2:1794160/3‑100 (MQ=255)
                                                                                          gacTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCg           >  1:904845/3‑100 (MQ=255)
                                                                                          gacTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCg           >  1:2231773/3‑100 (MQ=255)
                                                                                          gacTTTCCGGCAATGGCGGAAAGTCGAAACATGAAATCAAAACCTGGGTCGCCGCTATTGCCGCTAATTCTGAGTTTGGCAACTGGCGTAg                    >  1:1816965/3‑91 (MQ=255)
                                                                                          agCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCg           <  2:1464667/100‑1 (MQ=255)
                                                                                           gCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGt          <  1:50697/100‑1 (MQ=255)
                                                                                           acTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGt          >  2:1860333/2‑100 (MQ=255)
                                                                                           acTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGt          >  2:3242357/2‑100 (MQ=255)
                                                                                           aatttCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGt          >  2:4080411/3‑100 (MQ=255)
                                                                                            cTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAgaccg           >  1:2801480/1‑94 (MQ=255)
                                                                                            cTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGt          <  2:3060519/99‑1 (MQ=255)
                                                                                              ttCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGttat       >  2:3988087/1‑100 (MQ=255)
                                                                                                ccGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGttatta     >  1:1675769/1‑100 (MQ=255)
                                                                                                ccGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCCACTGGCGTAGCGAAGGGCGttatta     >  1:3746255/1‑100 (MQ=255)
                                                                                                ccGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGttatta     >  1:2748786/1‑100 (MQ=255)
                                                                                                ccGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGttatta     >  2:1114393/1‑100 (MQ=255)
                                                                                                 cGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTAc    <  2:2870818/100‑1 (MQ=255)
                                                                                                 cGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTAc    >  1:3510971/1‑100 (MQ=255)
                                                                                                 cGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTAc    <  1:1843722/100‑1 (MQ=255)
                                                                                                   ccATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTAcc   <  1:2781716/99‑1 (MQ=255)
                                                                                                   |                                                                                                   
TGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCG  >  NC_000913/370965‑371163

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: