Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 371,064 C→G 16.0% 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.0% 418.7 / 69.0 158A263G (GCC→GGC) mhpB2,3‑dihydroxyphenylpropionate 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (77/54);  new base G (14/11);  total (92/66)
Fisher's exact test for biased strand distribution p-value = 8.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

Only 100 of 162 total aligned reads displayed.

TGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCG  >  NC_000913/370965‑371163
                                                                                                   |                                                                                                   
tGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGc                                                                                                     >  1:2590236/1‑100 (MQ=255)
 gCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGcc                                                                                                    >  1:742920/1‑100 (MQ=255)
    tCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCAtt                                                                                                 >  1:1049380/1‑100 (MQ=255)
     ccGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTg                                                                                                <  1:3067384/100‑1 (MQ=255)
      cGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGc                                                                                               <  1:2224197/100‑1 (MQ=255)
      cGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGc                                                                                               >  1:620191/1‑100 (MQ=255)
       gCTCAACCCGCTTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGcc                                                                                              >  1:468977/1‑100 (MQ=255)
       gCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGcc                                                                                              >  2:2089386/1‑100 (MQ=255)
          cAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTGAGTAACGAAGAGCTTTCCGCCATTGCCgg                                                                                            >  2:25065/1‑99 (MQ=255)
          cAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGa                                                                                           >  2:585948/1‑100 (MQ=255)
            aCCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGaaa                                                                                         >  1:465261/1‑100 (MQ=255)
             cccGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAg                                                                                        >  2:1750684/1‑100 (MQ=255)
              ccGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAAGGACGCATACAGGAACTGGGTGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGAAAgc                                                                                       <  2:1988009/100‑2 (MQ=255)
               cGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGAAAgct                                                                                      <  1:2789515/100‑3 (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)
                  tttGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGac                                                                                   <  1:1460443/100‑1 (MQ=255)
                  tttGGGATAACCAGGTCATGACTTTGCTGGAGCAGGGACGCATACAGGAGCTGGATGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGAAAgctct                                                                                    <  2:43145/99‑5 (MQ=38)
                   ttGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGAAAgctcttc                                                                                  <  1:2268103/100‑7 (MQ=38)
                   ttGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGaca                                                                                  <  2:2592400/100‑1 (MQ=255)
                   ttGGGATAACCAGTTAATGACTTTGCTGTAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTTTCGAGGAGCTTTCCGCCTTTGCCg                                                                                             <  2:1130364/89‑1 (MQ=255)
                    tGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGAAAgctcttcg                                                                                 <  1:119811/100‑8 (MQ=38)
                    tGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGacac                                                                                 >  1:1781088/1‑100 (MQ=255)
                    tGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGacac                                                                                 >  1:157740/1‑100 (MQ=255)
                    tGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTACGCCATTGCCGGAAAGTCGacac                                                                                 >  1:609153/1‑100 (MQ=255)
                     gggATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGAAAgctcttcgt                                                                                <  2:2725612/100‑9 (MQ=37)
                     gggATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGAAAgctcttcgt                                                                                <  2:770144/100‑9 (MQ=37)
                       gATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATg                                                                              <  2:2685181/100‑1 (MQ=255)
                       gATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATg                                                                              <  1:585950/100‑1 (MQ=255)
                       gATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATg                                                                              >  1:1623456/1‑100 (MQ=255)
                            ccAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATc                                                                         <  1:2945858/100‑1 (MQ=255)
                            ccAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATc                                                                         >  1:1972263/1‑100 (MQ=255)
                            ccAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATc                                                                         <  1:2521628/100‑1 (MQ=255)
                             cAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCa                                                                        >  1:157537/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)
                                  cATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAAcc                                                                   <  2:3199962/100‑1 (MQ=255)
                                    tGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTg                                                                 <  2:626994/100‑1 (MQ=255)
                                      aCTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTggg                                                               <  1:3106809/100‑1 (MQ=255)
                                            cTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTcgccg                                                         >  1:2733688/1‑100 (MQ=255)
                                            cTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTcgccg                                                         >  1:1978329/1‑100 (MQ=255)
                                             tGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTcgccgc                                                        >  2:1231882/1‑100 (MQ=255)
                                              ggAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCt                                                       >  2:1661519/1‑100 (MQ=255)
                                                aGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCttt                                                     <  1:56183/100‑1 (MQ=255)
                                                   aGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGc                                                  >  2:1330456/1‑100 (MQ=255)
                                                    gggACGCATACAGGAACTGGATGCCGTCAGTAAGGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACATGGGTCGCCGCTCTTg                                                   >  2:3015335/1‑98 (MQ=255)
                                                    gggACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAACGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGcc                                                 <  2:2043459/100‑1 (MQ=255)
                                                    gggACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGcc                                                 >  2:3614526/1‑100 (MQ=255)
                                                    gggACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGACATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGGCGCTTTTGcc                                                 >  2:3294339/1‑100 (MQ=255)
                                                        cGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTa                                             >  2:1070379/1‑100 (MQ=255)
                                                          cATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTAtt                                           >  1:1096524/1‑100 (MQ=255)
                                                          cATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTAtt                                           >  2:2526433/1‑100 (MQ=255)
                                                          cATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTAtt                                           >  2:2949046/1‑100 (MQ=255)
                                                           aTACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTAttt                                          >  1:3040357/1‑100 (MQ=255)
                                                            tACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTc                                         <  1:3140445/100‑1 (MQ=255)
                                                            tACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTc                                         <  1:281357/100‑1 (MQ=255)
                                                            tACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTc                                         >  1:1950826/1‑100 (MQ=255)
                                                              cAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTg                                       >  2:507961/1‑100 (MQ=255)
                                                               aGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGc                                      <  1:3529530/100‑1 (MQ=255)
                                                                ggAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCg                                     >  1:2021577/1‑100 (MQ=255)
                                                                 gAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGt                                    >  2:2363723/1‑100 (MQ=255)
                                                                 gAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGt                                    >  1:1715239/1‑100 (MQ=255)
                                                                 gAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGt                                    >  1:963189/1‑100 (MQ=255)
                                                                  aaCTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGtt                                   <  2:852862/100‑1 (MQ=255)
                                                                  aaCTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGtt                                   <  1:838476/100‑1 (MQ=255)
                                                                     tGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGttt                                  <  2:282895/98‑1 (MQ=255)
                                                                     tGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTgg                                <  1:2472457/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)
                                                                        aTGCCGTCAGTAACGAAGAGCTTTCCGCCAGTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTgg                                <  2:2333986/97‑1 (MQ=255)
                                                                           ccGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTg                          >  2:2083381/1‑100 (MQ=255)
                                                                              tCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCg                       >  1:2763922/1‑100 (MQ=255)
                                                                                 gTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAg                    >  1:1647309/1‑100 (MQ=255)
                                                                                 gTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGATTGGCAACTGGCGTAg                    >  2:1717227/1‑100 (MQ=255)
                                                                                 gTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTCTTGCAGCTATTTCTGCGTTTGGCAACTGGCGTAg                    >  2:150252/1‑100 (MQ=255)
                                                                                  tAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGc                   >  1:3417231/1‑100 (MQ=255)
                                                                                      gAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAg               >  2:3433818/1‑100 (MQ=255)
                                                                                       gtcgacTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAgg              >  1:3635071/6‑100 (MQ=38)
                                                                                          gacTTTCCGGCAATGGCGGAAAGTCGACACATGGAATCAAAACCTGGGTCGCCGCTTGTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGc            >  2:2226491/3‑99 (MQ=255)
                                                                                          gacTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCg           >  2:1500670/3‑100 (MQ=255)
                                                                                          gacTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCg           >  1:890556/3‑100 (MQ=255)
                                                                                           gCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCg           >  2:840293/1‑99 (MQ=255)
                                                                                           acTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCg           >  2:1008446/2‑99 (MQ=255)
                                                                                            cTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGtt         <  2:1978318/100‑1 (MQ=255)
                                                                                            cTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGtt         >  1:2202903/1‑100 (MQ=255)
                                                                                             tttCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCg                  >  2:2887899/1‑90 (MQ=255)
                                                                                             tttCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGtta        <  1:105750/100‑1 (MQ=255)
                                                                                             tttCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGtta        <  1:3204130/100‑1 (MQ=255)
                                                                                              ttCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGttat       <  2:957041/100‑1 (MQ=255)
                                                                                               tCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGttatt      >  1:1581804/1‑100 (MQ=255)
                                                                                               tCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGttatt      <  1:3433843/100‑1 (MQ=255)
                                                                                                ccGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGttatta     >  2:2201809/1‑100 (MQ=255)
                                                                                                 cGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTAc    <  1:2526449/100‑1 (MQ=255)
                                                                                                  gCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTAcc   >  2:3420007/1‑100 (MQ=255)
                                                                                                   ccATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCg  <  1:540357/100‑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: