Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,339,802 A→G 19.2% A240A (GCT→GCC adiA ← arginine decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,339,8020AG19.2% 428.0 / 95.1 210A240A (GCT→GCCadiAarginine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (91/78);  new base G (24/16);  total (116/94)
Fisher's exact test for biased strand distribution p-value = 5.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

Only 100 of 212 total aligned reads displayed.

ACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCA  >  NC_000913/4339704‑4339901
                                                                                                  |                                                                                                   
aCCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAg                                                                                                    <  1:3818183/100‑1 (MQ=255)
aCCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAg                                                                                                    <  2:1482124/100‑1 (MQ=255)
 ccATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGc                                                                                                   >  1:3595829/1‑100 (MQ=255)
 ccATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGc                                                                                                   >  1:374062/1‑100 (MQ=255)
 ccATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGc                                                                                                   <  2:159723/100‑1 (MQ=255)
  cATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGcc                                                                                                  >  1:1915544/1‑100 (MQ=255)
   atatAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCt                                                                                                 >  2:1859916/1‑100 (MQ=255)
    tataGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTg                                                                                                <  1:2034938/100‑1 (MQ=255)
    tataGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGACTg                                                                                                >  2:2850370/1‑100 (MQ=255)
    tataGACCGGGTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACTACCACGACATCGTTATCGGTCATGCAAGCCTg                                                                                                >  2:562403/1‑100 (MQ=255)
     ataGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGc                                                                                               <  2:775902/100‑1 (MQ=255)
       aGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCat                                                                                             <  2:825815/100‑1 (MQ=255)
        gACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCatg                                                                                            <  1:1116529/100‑1 (MQ=255)
          ccGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCGGTTACGGTGAACGACCAGGACATCGTTATCGGTCATGCAGGCTTGCatgac                                                                                          <  2:1718521/100‑2 (MQ=255)
          ccGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAGGCTTGCatgac                                                                                          <  1:1026732/100‑2 (MQ=255)
          ccGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCatgat                                                                                          <  2:51976/100‑1 (MQ=255)
          ccGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCatgat                                                                                          >  1:2605477/1‑100 (MQ=255)
          ccGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCatgat                                                                                          <  1:2677093/100‑1 (MQ=255)
          ccGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGAGGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCatgat                                                                                          <  2:2102107/100‑1 (MQ=255)
            ggTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAGGCTTGCatgaccg                                                                                        <  1:2495781/100‑4 (MQ=255)
            ggTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATgg                                                                                        <  1:3465632/100‑1 (MQ=255)
             gTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAGGCTTGCatgaccga                                                                                       <  1:2190522/100‑5 (MQ=255)
             gTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAGGCTTGCatgaccga                                                                                       <  1:2229926/100‑5 (MQ=255)
             gTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAGGCTTGCatgaccga                                                                                       <  2:769304/100‑5 (MQ=255)
             gTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGt                                                                                       >  2:645677/1‑100 (MQ=255)
             gTTGCGCGCCTGTCAGCATCAAACCTTGTTCGATGGCTGTATGGCAGTTACGGTCAACGTCCACGACATCGTTATCGGTCATGCAGGCTTGCatgaccg                                                                                        <  2:236119/99‑4 (MQ=255)
               ttCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGc                                                                                     >  1:300802/1‑100 (MQ=255)
                 cgcgCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCgg                                                                                   <  2:940201/100‑1 (MQ=255)
                 cgcgCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCgg                                                                                   >  1:1843835/1‑100 (MQ=255)
                  gcgcCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGt                                                                                  <  1:2954121/100‑1 (MQ=255)
                   cgcCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGtt                                                                                 >  1:3282995/1‑100 (MQ=255)
                   cgcCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGtt                                                                                 >  1:425392/1‑100 (MQ=255)
                   cgcCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACAACGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGtt                                                                                 >  2:2567214/1‑100 (MQ=255)
                      cTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTaga                                                                              <  1:142074/100‑1 (MQ=255)
                       tGTCAGCATCAAACCTTGTTCGATGGATTTATCGCAGTTACGGTCAACGACCACGACATCGTTACGGGTCATGCAAGCCTGCATGATCGTGAGGTTagag                                                                             >  2:1559190/1‑100 (MQ=255)
                              aTCAAACCTTGTTCGATGGATTTATGTCAGTTACGGTTAACGCCCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGaa                                                                       <  2:2305731/99‑1 (MQ=255)
                                  aaCCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTAcc                                                                  <  1:3117651/100‑1 (MQ=255)
                                    ccTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACcga                                                                >  1:118068/1‑100 (MQ=255)
                                    ccTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACcga                                                                >  2:2251643/1‑100 (MQ=255)
                                      ttGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACcgacg                                                              <  2:1915533/100‑1 (MQ=255)
                                      ttGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACcgacg                                                              <  2:2288797/100‑1 (MQ=255)
                                      ttGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACcgac                                                               >  2:864546/1‑99 (MQ=255)
                                        gTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCgacgac                                                            <  2:349314/100‑1 (MQ=255)
                                          tCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCATGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTa                                                          >  2:2451246/1‑100 (MQ=255)
                                           cGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTAc                                                         >  2:2909186/1‑100 (MQ=255)
                                            naTGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCAAGATGGCGCGGTTAGGGCCGGAAGTACCGACGACTAcc                                                        >  2:2579056/2‑100 (MQ=255)
                                            gATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTAcc                                                        >  2:900701/1‑100 (MQ=255)
                                            gATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTAcc                                                        <  1:1799814/100‑1 (MQ=255)
                                            gATGGATATATGGCAGGTACGGTCAACGACCCCGACCTCGTTAGAGGTCCTGCAAGCCTGCATGATTGTGCGGCTAGAGCCGGAAGTACctacgcccac                                                         >  2:696986/1‑94 (MQ=255)
                                             aTGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCg                                                       >  1:283158/1‑100 (MQ=255)
                                                gATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGAcc                                                    <  1:3097806/100‑1 (MQ=255)
                                                gATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGAcc                                                    >  1:3534715/1‑100 (MQ=255)
                                                gATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGAcc                                                    >  1:2771281/1‑100 (MQ=255)
                                                     aTGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGTATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAg                                               >  1:3353927/1‑100 (MQ=255)
                                                     aTGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAg                                               <  1:141555/100‑1 (MQ=255)
                                                     aTGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAg                                               >  1:87719/1‑100 (MQ=255)
                                                       ggCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCg                                             >  1:50244/1‑100 (MQ=255)
                                                          aGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATc                                          <  2:304339/100‑1 (MQ=255)
                                                          aGTTACGGTCAACGACCACGACATCGTTATCGGTCAGGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATc                                          <  2:199979/100‑1 (MQ=255)
                                                          aGTTACGGTCAACGACCACGACAGCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATc                                          <  1:2924152/100‑1 (MQ=255)
                                                          aGTGACGGTCAACGACCACTACATCGTTATCGGTCATGCAATCCTGCATGATGGTGCGGTTAGAGCCGTAAGTACCGACGACTACCGACCAGGAGCGATc                                          <  2:2211184/100‑1 (MQ=255)
                                                           gTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCg                                         >  2:1001469/1‑100 (MQ=255)
                                                             tACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCgtc                                       <  1:2580251/100‑3 (MQ=255)
                                                             tACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGc                                       >  2:3130053/1‑100 (MQ=255)
                                                              aCGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGTCTACCGACCAGGAGCGATCGGc                                       <  2:1380049/99‑1 (MQ=255)
                                                              aCGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGAGATCGGCa                                      >  2:2172543/1‑100 (MQ=255)
                                                                 gTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGCCCAGGAGCGATCGGCACCa                                   >  2:3717713/1‑100 (MQ=255)
                                                                  tCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCa                                   <  2:3208969/99‑1 (MQ=255)
                                                                   cAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGAAGGTGCGGTTAGAGCCGGAAGTCCCGACCACTCACGACCAGGAGCGATCGGCACCaaa                                 >  2:1242894/1‑100 (MQ=255)
                                                                     aCGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATa                               <  2:2495158/100‑1 (MQ=255)
                                                                      cGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATAc                              >  2:128395/1‑100 (MQ=255)
                                                                       gACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATAcg                             <  2:1028544/100‑1 (MQ=255)
                                                                        aCCACGAGATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATAcgc                            >  1:2004655/1‑100 (MQ=255)
                                                                         ccACGACATCGTTATCGGTCATGCAGGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATAcgc                            >  2:1970142/1‑99 (MQ=255)
                                                                            cGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGc                        <  2:1112618/100‑1 (MQ=255)
                                                                             gACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCg                       <  1:3634684/100‑1 (MQ=255)
                                                                               cATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGc                     >  1:1203217/1‑100 (MQ=255)
                                                                                 aaccgcaccatCATGCAGGCTTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCat                   >  1:242872/11‑100 (MQ=38)
                                                                                 aaccgcaccatCATGCAGGCTTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCat                   >  1:1868535/11‑100 (MQ=38)
                                                                                   nttATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCatat                 >  2:2725465/2‑100 (MQ=255)
                                                                                   gTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCatat                 <  2:3367428/100‑1 (MQ=255)
                                                                                     tATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATAttt               >  1:3429053/1‑100 (MQ=255)
                                                                                       accatCATGCAGGCTTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATAttttt             >  1:955290/5‑100 (MQ=255)
                                                                                        ccatCATGCAGGCTTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttc            >  1:136777/4‑100 (MQ=255)
                                                                                        ccatCATGCAGGCTTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttc            >  1:1487798/4‑100 (MQ=255)
                                                                                        ccatCATGCAGGCTTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttc            >  2:159069/4‑100 (MQ=255)
                                                                                        cGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCTACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttc            >  2:2434264/1‑100 (MQ=255)
                                                                                         ggTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcg           >  1:2654050/1‑100 (MQ=255)
                                                                                         ggTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcg           <  1:209388/100‑1 (MQ=255)
                                                                                         ggTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcg           <  1:2220233/100‑1 (MQ=255)
                                                                                         ggTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcg           <  1:672510/100‑1 (MQ=255)
                                                                                         ggTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcg           <  1:556789/100‑1 (MQ=255)
                                                                                          atCATGCAGGCTTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCACATACGCGTCCGGCCTATTtttcg           >  2:2825975/2‑99 (MQ=255)
                                                                                          atCATGCAGGCTTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcgc          >  1:83830/2‑100 (MQ=255)
                                                                                          atCATGCAGGCTTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcgc          >  1:2291188/2‑100 (MQ=255)
                                                                                            cATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcgctt        >  1:1744929/1‑100 (MQ=255)
                                                                                             aTGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcgcttt       >  1:1493876/1‑100 (MQ=255)
                                                                                               gCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcgctttcg     >  1:759572/1‑100 (MQ=255)
                                                                                                cAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcgctttcgc    >  2:355651/1‑100 (MQ=255)
                                                                                                 aaGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGcc   <  2:1502682/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
ACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCA  >  NC_000913/4339704‑4339901

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: