Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,311,831 C→T 15.9% intergenic (+23/+17) tonB → / ← yciA membrane spanning protein in TonB‑ExbB‑ExbD transport complex/acyl‑CoA esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,311,8310CT15.9% 313.3 / 42.1 131intergenic (+23/+17)tonB/yciAmembrane spanning protein in TonB‑ExbB‑ExbD transport complex/acyl‑CoA esterase
Reads supporting (aligned to +/- strand):  ref base C (51/58);  new base T (10/11);  total (62/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 136 total aligned reads displayed.

ATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTT  >  NC_000913/1311732‑1311930
                                                                                                   |                                                                                                   
aTGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGAc                                                                                                     <  2:3160491/100‑1 (MQ=255)
 tGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGAcc                                                                                                    >  1:248577/1‑100 (MQ=255)
  gAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCg                                                                                                   >  1:3503041/1‑100 (MQ=255)
  gAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCg                                                                                                   <  1:2505554/100‑1 (MQ=255)
   aGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCgg                                                                                                  <  1:1622913/100‑1 (MQ=255)
   aGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCgg                                                                                                  >  1:3223386/1‑100 (MQ=255)
   aGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCgg                                                                                                  <  1:733347/100‑1 (MQ=255)
    gCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGa                                                                                                 <  1:3709693/100‑1 (MQ=255)
    gCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGa                                                                                                 >  1:1610210/1‑100 (MQ=255)
    gCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGa                                                                                                 >  1:600426/1‑100 (MQ=255)
    gCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGa                                                                                                 <  2:2396330/100‑1 (MQ=255)
     ccGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAg                                                                                                <  1:4139729/100‑1 (MQ=255)
      cGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAgg                                                                                               >  1:793134/1‑100 (MQ=255)
       gggTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGc                                                                                              >  1:2523035/1‑100 (MQ=255)
         gTAAACCAGGCAGTGGGATTGTGCAGAATATCCTGTTTAAAATTAACGGCACCTCCGAAATTCAGTAATCAGAAAGTCAAAAGCCTCCGACCGGAGGc                                                                                              <  2:395740/98‑1 (MQ=255)
           aaGCCAGTCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCtttt                                                                                          <  2:2064309/100‑1 (MQ=39)
           aaGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCtttt                                                                                          <  2:2447328/100‑1 (MQ=255)
           aaGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCtttt                                                                                          <  2:3657614/100‑1 (MQ=255)
            agggaGGGAGTGGGATTGTGGTGGATGTCCTGTTTAAAATTAACGGCACCACCGAGATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTg                                                                                         <  2:3436339/96‑1 (MQ=255)
            aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCTAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTg                                                                                         <  2:2321395/100‑1 (MQ=39)
            aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTg                                                                                         <  1:817981/100‑1 (MQ=255)
            aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTg                                                                                         <  1:2603262/100‑1 (MQ=255)
            aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTg                                                                                         <  1:2019738/100‑1 (MQ=255)
            aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTg                                                                                         <  1:2984079/100‑1 (MQ=255)
            aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTg                                                                                         <  2:3436587/100‑1 (MQ=255)
            aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTg                                                                                         <  2:3962342/100‑1 (MQ=255)
            aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTg                                                                                         >  1:1569524/1‑100 (MQ=255)
             gCCAGGCAGTGGGATTGTGGTGAATATCCTTTTTAAAATTAACGGCACCACCGAGATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTg                                                                                         <  2:5035/99‑1 (MQ=255)
             gCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGa                                                                                        >  2:3407940/1‑100 (MQ=255)
             gCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGa                                                                                        >  1:322835/1‑100 (MQ=255)
             gCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGa                                                                                        >  2:236835/1‑100 (MQ=255)
              ccAGGCAGTGGGATTGTGGTGAATAGCCTGTTTAAAATTAACGCCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGAc                                                                                       <  2:1223342/100‑1 (MQ=255)
                aGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTGACtt                                                                                     <  1:505009/100‑2 (MQ=255)
                 ggCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTAt                                                                                    >  1:4107834/1‑100 (MQ=255)
                 ggCAGGGGGATTGTTGTGAATATCCTGTGTAAAGTGAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGAc                                                                                       <  2:1222388/97‑1 (MQ=255)
                  gCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTGACtttc                                                                                   <  1:58614/100‑4 (MQ=38)
                   cAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTa                                                                                  <  2:1263420/100‑1 (MQ=255)
                      tGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTc                                                                               <  1:3830636/100‑1 (MQ=255)
                             gtggtgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGt                                                                        <  2:1611825/100‑1 (MQ=255)
                              tgttgAATATCTTGTTTAAAATTACAGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTa                                                                       <  2:1688707/97‑1 (MQ=255)
                              tggtgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTa                                                                       >  1:2122849/1‑100 (MQ=255)
                                gtgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAg                                                                     <  2:2902393/100‑1 (MQ=255)
                                gtgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAg                                                                     >  1:1926149/1‑100 (MQ=255)
                                 tgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAgg                                                                    <  2:3222300/100‑1 (MQ=255)
                                 tgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAgg                                                                    >  1:1082008/1‑100 (MQ=255)
                                       tCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGa                                                               <  1:95408/99‑1 (MQ=255)
                                       tCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAg                                                              <  2:2653853/100‑1 (MQ=255)
                                         cTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGt                                                            <  2:1673611/100‑1 (MQ=255)
                                           gTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGttt                                                          >  1:3386681/1‑100 (MQ=255)
                                             atAAAATTAACGGCACCACCGAAATTCAGTAATCAGAAAGTCAAAAGCCGCCGACCGGAGGCGTTTGACTATTACTCAACAGGTGAGGCGCGAGGttt                                                          >  2:1393332/2‑98 (MQ=255)
                                              tAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTTAAAAGCCTCCGACCGGAGGCTTTTGACTATTCATCAACAGGTAAGGCGCGAGGattt                                                         >  2:3751632/1‑98 (MQ=255)
                                              tAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTcc                                                       >  2:2884554/1‑100 (MQ=255)
                                              tAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTcc                                                       >  2:3567589/1‑100 (MQ=255)
                                              tAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTc                                                        >  2:352131/1‑99 (MQ=255)
                                                  aTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCa                                                   >  2:2971690/1‑100 (MQ=255)
                                                     aaCGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGa                                                >  1:3367042/1‑100 (MQ=255)
                                                      aCGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGAt                                               <  2:1480362/100‑1 (MQ=255)
                                                        ggCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCa                                             <  1:509659/100‑1 (MQ=255)
                                                        ggCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCa                                             <  1:1502775/100‑1 (MQ=255)
                                                             caccGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACcgc                                        >  2:454959/1‑100 (MQ=255)
                                                              accGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACcgcg                                       <  1:858613/100‑1 (MQ=255)
                                                                cGAAATTCAGTAAGCAGAAAGTCAAATGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGAc                                     <  2:2223251/100‑1 (MQ=255)
                                                                cGAAATTCAGTAAGCAGAAAGTCAAACGCCTCCGACCGGAGGCTTTTGACTATTACTCGACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACcgcg                                       <  1:1874168/98‑1 (MQ=255)
                                                                 gAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACa                                    >  2:4001204/1‑100 (MQ=255)
                                                                 gAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACa                                    >  1:2040193/1‑100 (MQ=255)
                                                                 gAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACa                                    <  2:486910/100‑1 (MQ=255)
                                                                      tCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACtt                               >  1:2734456/1‑100 (MQ=255)
                                                                       cAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTa                              <  1:2451723/100‑1 (MQ=255)
                                                                       cAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTa                              <  2:3503024/100‑1 (MQ=255)
                                                                        aGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTCTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTaa                             >  1:2964832/1‑100 (MQ=255)
                                                                         gttgagtaataGTCAAAAGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTaaa                            >  1:3577764/11‑100 (MQ=25)
                                                                         gTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTaaa                            >  1:84350/1‑100 (MQ=255)
                                                                         gTAAGCAGAAAATCAAAAGCCTCCGACCGGAGGCTTTCGACTATTACTCAACAGGTAAGGCGCGAGGGTTTCCATCAGGATCAACCGCGAAATACTaaaa                            >  2:1719573/1‑100 (MQ=255)
                                                                          tAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTTACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAaat                           <  2:1052863/100‑1 (MQ=255)
                                                                          tAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAaat                           <  1:3206767/100‑1 (MQ=255)
                                                                            aGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAaataa                         >  2:2479719/1‑100 (MQ=255)
                                                                              cAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGCGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATg                       >  1:102808/1‑100 (MQ=255)
                                                                               aGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGc                      <  2:3527956/100‑1 (MQ=255)
                                                                                gAAAGTCAAAAGCCTCCGACCGGAGGGTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCt                     >  1:2550114/1‑100 (MQ=255)
                                                                                 aaaGTCAAAAGCCTCCGACCGGAGTCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCtt                    <  2:621442/100‑1 (MQ=255)
                                                                                 aaaGTCAAAAGCCTCCGACCGGAGGCTTTTGACGATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCtn                    <  2:2517016/100‑2 (MQ=255)
                                                                                   aGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCt                  <  2:2523028/100‑1 (MQ=255)
                                                                                    ggcAAAAGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACCGCTAAGGCGCGAGGTTTTCGTTCAGGGTCAACCGCGACATACTTaactta                         >  2:181132/3‑88 (MQ=25)
                                                                                    gTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTg                 >  1:466268/1‑100 (MQ=255)
                                                                                    gTCAAAAGCCTACGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTGTTCCTTCAGGATCAACCGCGACATACTTAAATAATCCTTCt                  >  2:3822654/1‑99 (MQ=38)
                                                                                     tCAAAAGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGt                >  1:1333172/1‑100 (MQ=255)
                                                                                      cAAAAGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTc               >  1:1952645/1‑100 (MQ=255)
                                                                                       aaaaGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCg              <  2:1569521/100‑1 (MQ=255)
                                                                                       aaaaGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCg              <  1:4001231/100‑1 (MQ=255)
                                                                                        aaaGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGc             <  2:272946/100‑1 (MQ=255)
                                                                                          aGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCtt           <  1:1042831/100‑1 (MQ=255)
                                                                                           gCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCttt          >  1:2627674/1‑100 (MQ=255)
                                                                                           gCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCttt          <  2:1766704/100‑1 (MQ=255)
                                                                                            ccTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTta         >  1:3378931/1‑100 (MQ=255)
                                                                                             cTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTtat        <  1:3567606/100‑1 (MQ=255)
                                                                                              tCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTtata       <  2:3551981/100‑1 (MQ=255)
                                                                                                 gACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCg    >  2:1337708/1‑100 (MQ=255)
                                                                                                  aCCGGAGGCTTTTGACTCTTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGt   >  2:608172/1‑100 (MQ=255)
                                                                                                   ccGGATGCTTTTGACTAGTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGtt  <  1:1947799/100‑1 (MQ=255)
                                                                                                   ccGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGtt  >  2:3102945/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
ATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTT  >  NC_000913/1311732‑1311930

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: