Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,051,259 G→T 21.4% intergenic (‑494/‑88) yihA ← / → yihI cell division GTP‑binding protein/activator of Der GTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,051,2590GT21.4% 482.9 / 130.0 216intergenic (‑494/‑88)yihA/yihIcell division GTP‑binding protein/activator of Der GTPase
Reads supporting (aligned to +/- strand):  ref base G (85/84);  new base T (27/19);  total (113/103)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-02

Only 100 of 227 total aligned reads displayed.

CACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCA  >  NC_000913/4051160‑4051358
                                                                                                   |                                                                                                   
cacGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACagag                                                                                                     >  1:1206752/1‑100 (MQ=255)
cacGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACagag                                                                                                     >  2:450826/1‑100 (MQ=255)
 acGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGt                                                                                                    <  1:3171393/100‑1 (MQ=255)
  cGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTa                                                                                                   >  1:3235313/1‑100 (MQ=255)
    gATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAAt                                                                                                 >  2:1001399/1‑100 (MQ=255)
      ttGTTCTGGCTAAAGAAAAAACAGGGTGTCTTGGCGGCATATAAGGATTGTAAGACCTGTTAAGGTTTAAAAAAAAAAAAAAAAGGCGACAGAGTAATCt                                                                                               <  1:3868236/100‑1 (MQ=255)
        gTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGt                                                                                             >  2:3807527/1‑100 (MQ=255)
          tCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCg                                                                                           <  2:3664033/100‑1 (MQ=255)
           ccGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGc                                                                                          <  2:134343/100‑1 (MQ=255)
           ccGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGc                                                                                          <  1:2838606/100‑1 (MQ=255)
           ccGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGc                                                                                          <  2:2166058/100‑1 (MQ=255)
            cGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGcc                                                                                         <  2:2549547/100‑1 (MQ=255)
            cGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGcc                                                                                         <  2:916423/100‑1 (MQ=255)
            cGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCg                                                                                           <  1:3984985/98‑1 (MQ=255)
             ggCTAGAGGCAAAACAGGGGGTGTTGGCCGCCTGCAAGGATAGTAAGACCCGCAGAGGTGTATGATTCGGGAAAAAAGGCGACAGAGTAATCTGt                                                                                             >  1:3992130/1‑95 (MQ=255)
              gCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGCCtt                                                                                       <  2:165731/100‑1 (MQ=255)
              gCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCtt                                                                                       <  2:1134704/100‑1 (MQ=255)
               cTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGCCttt                                                                                      <  2:1726637/100‑1 (MQ=255)
                taaattaaaaaTATGTTGTGTTGGCGGTCTGCAAGGATTGTAAGATCCGTTAAGGGTTATGAGTTAGGAAAAAAGGCGACAGATTACTCTGTCGCCt                                                                                        <  1:3379394/91‑1 (MQ=25)
                tAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCtttt                                                                                     >  1:1461737/1‑100 (MQ=255)
                  aaGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCtttttt                                                                                   <  2:900122/100‑1 (MQ=255)
                   aGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGCCTTTTTTc                                                                                  <  2:514435/100‑1 (MQ=255)
                   aGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCtttttt                                                                                   >  2:4144202/1‑99 (MQ=255)
                   aGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTc                                                                                  >  1:3129547/1‑100 (MQ=255)
                   aGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTc                                                                                  >  2:4313788/1‑100 (MQ=255)
                     gAAAAAGAGGGTGTGTTGGGGGCCTGGAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGCCTTTTTTcct                                                                                <  1:1278502/100‑3 (MQ=37)
                       aaaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttg                                                                              >  2:661506/1‑100 (MQ=255)
                       aaaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttg                                                                              >  1:1042266/1‑100 (MQ=255)
                        aaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTGATGAGTCAGTAAAAAAGGCGACAGGGTGATCTGTCGCCTGTTGTCTttgc                                                                             >  2:3507492/1‑100 (MQ=255)
                         aaCAGGGTGTATTGGCGGCCTGCAAGGATCGTAAGACCGGTTAGGGGTTATGCGTCAGGCAGAAAGACGACGGAGTAATCTGTCGCCTTTTTTCTttgct                                                                            >  1:2787324/1‑100 (MQ=255)
                           cAGGTTGTGTTGGCGTCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgcttg                                                                          <  2:363416/100‑1 (MQ=255)
                                tgtgtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTc                                                                     >  2:2396884/1‑100 (MQ=255)
                                tgtgtTGGCGGCCTGCAAGGACTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTc                                                                     <  2:120793/100‑1 (MQ=255)
                                tgtgtTGGCGCACTTCAAGGACTCTAGGCACCGTCAATTGTTATCAGTCAGGAAAAAAGGCGACAGAGTTTTTTGTCGCCTTTCTTGtt                                                                                <  1:3782540/89‑1 (MQ=255)
                                 gtgtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTATTTGCTTGCTTTCt                                                                    >  1:2147216/1‑100 (MQ=255)
                                  tgtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTg                                                                   <  2:1585642/100‑1 (MQ=255)
                                   gtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCt                                                                    >  2:1917179/1‑98 (MQ=255)
                                   gtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGt                                                                  <  1:366502/100‑1 (MQ=255)
                                    ttGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGtt                                                                 <  2:2345468/100‑1 (MQ=255)
                                      ggcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAg                                                               <  2:3195429/100‑1 (MQ=255)
                                      ggcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAg                                                               >  2:4103252/1‑100 (MQ=255)
                                       gcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTGATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGa                                                              >  1:1992847/1‑100 (MQ=255)
                                           ccTGGAAGGATTGTAAGACCCGTTAAGGTTTATGAGTCAGTAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTa                                                                <  1:4062587/94‑1 (MQ=255)
                                               cAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgcc                                                       >  2:116501/1‑99 (MQ=255)
                                                aaGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCcgccgc                                                     <  1:3488361/100‑1 (MQ=255)
                                                 aGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCcgccgc                                                     >  1:4137214/1‑99 (MQ=255)
                                                 aGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCa                                                    >  1:1527124/1‑100 (MQ=255)
                                                  ggATTTTAAGACCCGTTAAGGGTTATGAGTCATGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCTCCGCaa                                                   <  1:2404014/100‑1 (MQ=255)
                                                   gATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACATAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCaaa                                                  >  1:2441142/1‑100 (MQ=255)
                                                   gATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCaaa                                                  <  1:4276915/100‑1 (MQ=255)
                                                   gATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCaaa                                                  <  2:4211577/100‑1 (MQ=255)
                                                     ttGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATc                                                <  2:304379/100‑1 (MQ=255)
                                                      tGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCATTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCt                                               <  1:3507474/100‑1 (MQ=255)
                                                        tAAGACCCGTTCAGGGTTATGAGTCAGGAAAAACGGCCACACCGTAATCTGCCGCCCCCTTTCTTTGCTCGCTTTCTGTTAGAATCCGCCGCAAATCtat                                             <  1:1296138/100‑1 (MQ=255)
                                                         aaGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCtata                                            >  2:820821/1‑100 (MQ=255)
                                                         aaGACCCGTTAAGGGTTATGAATCAGGAAAAAAGGCGACAGAGTAATCTGTCTCCTTTTTTATTTGCGTGCTTTCTGTTAGATTCCGCCGCAAAtatata                                            >  1:1662505/1‑100 (MQ=38)
                                                           gACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACt                                          >  2:934296/1‑100 (MQ=255)
                                                            aCCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTg                                         <  2:2369230/100‑1 (MQ=255)
                                                             cccGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGa                                        >  1:541368/1‑100 (MQ=255)
                                                             cccGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGa                                        <  2:1246853/100‑1 (MQ=255)
                                                              cctttAAGGGTTATGAGTTAGGAAAAAAGGCGACAGAGTACTCTGTCGCCTTTTTTCTTTGGTTGATTTCTGTTAGATTGCGTCGCAAATCTATAc                                           <  1:3452671/93‑1 (MQ=25)
                                                               cGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAAt                                      <  1:2607362/100‑1 (MQ=255)
                                                                gTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATa                                     >  2:1600097/1‑100 (MQ=255)
                                                                gTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATa                                     >  2:2800658/1‑100 (MQ=255)
                                                                 ttAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATaa                                    <  2:746588/100‑1 (MQ=255)
                                                                   aaGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGACTaaaa                                  >  1:91351/1‑100 (MQ=255)
                                                                      ggTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACgg                               <  2:3889348/100‑1 (MQ=255)
                                                                           tGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTaaa                          >  1:1405453/1‑100 (MQ=255)
                                                                            gAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTaaaa                         >  2:1746823/1‑100 (MQ=255)
                                                                              gTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGa                       <  1:831793/100‑1 (MQ=255)
                                                                              gTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGa                       >  2:1854154/1‑100 (MQ=255)
                                                                              gTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGa                       <  2:3535771/100‑1 (MQ=255)
                                                                                cAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCTAATCTATACTTAATAAAACGGCTAAAAGACg                     >  1:1840463/1‑100 (MQ=255)
                                                                                  ggAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGaa                   >  2:3989136/1‑100 (MQ=255)
                                                                                     aaaaaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCa                >  1:4222009/1‑100 (MQ=255)
                                                                                     aaaaaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCa                >  2:201399/1‑100 (MQ=255)
                                                                                       aaaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCAtt              >  2:697241/1‑100 (MQ=255)
                                                                                       aaaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTGTCTGTTAGATTCCGCCGCTGATCTATACTGAATATAACGGCTAAAAGACGAAc                  >  1:2963678/1‑96 (MQ=37)
                                                                                        aaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTa             >  2:3520355/1‑100 (MQ=255)
                                                                                         aGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCTCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCAt               >  1:1730343/1‑97 (MQ=38)
                                                                                         aGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGCATCAAACGGCTAAAAGACGAACCATTAt            >  1:4032832/1‑100 (MQ=38)
                                                                                         aGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTAt            >  2:4371669/1‑100 (MQ=255)
                                                                                         aGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTAt            >  2:1543183/1‑100 (MQ=255)
                                                                                         aGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTAt            <  2:2828562/100‑1 (MQ=255)
                                                                                          ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg           >  2:2147046/1‑100 (MQ=255)
                                                                                          ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg           >  2:2736076/1‑100 (MQ=255)
                                                                                          ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg           >  2:3369016/1‑100 (MQ=255)
                                                                                          ggCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg           <  2:541370/100‑1 (MQ=255)
                                                                                          ggCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg           >  2:2588326/1‑100 (MQ=255)
                                                                                           gCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACAGATAAAAGACGAACAATTATGa          >  1:2188045/1‑100 (MQ=38)
                                                                                            cGACAGAGTAATCTGTCGCCTTTTTTTTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGaa         >  1:700082/1‑100 (MQ=255)
                                                                                            cGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGaa         <  1:229218/100‑1 (MQ=255)
                                                                                             gACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGaaa        <  2:13261/100‑1 (MQ=255)
                                                                                              aCAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAAc       <  1:321469/100‑1 (MQ=255)
                                                                                              aCAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAAc       <  1:1826326/100‑1 (MQ=255)
                                                                                               cAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAAcc      >  2:3121902/1‑100 (MQ=255)
                                                                                                 gagTAATCTGTCGCCTTTTTTTTTTGTTTTCTTTCTGTTAGATTTCGCCGCAAATTTATTTTGAATAATATGGCTAAAAGTCGAACCTTTTTGTAACcat    >  1:788561/1‑100 (MQ=255)
                                                                                                   gTAATCTGTCGCCTTTTTTTTTTTTTTTCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcatca  >  1:2458781/1‑100 (MQ=255)
                                                                                                   gTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcatca  >  2:4249066/1‑100 (MQ=255)
                                                                                                   gTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcatca  >  2:3724396/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCA  >  NC_000913/4051160‑4051358

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: