Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,051,259 G→T 20.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,2590GT20.4% 264.1 / 57.9 117intergenic (‑494/‑88)yihA/yihIcell division GTP‑binding protein/activator of Der GTPase
Reads supporting (aligned to +/- strand):  ref base G (51/41);  new base T (17/7);  total (68/49)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.07e-02

Only 100 of 125 total aligned reads displayed.

ACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCA  >  NC_000913/4051161‑4051358
                                                                                                  |                                                                                                   
atggATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAAAAAGAAAAAAAGGCGACAGAGt                                                                                                    <  2:2552014/98‑1 (MQ=255)
  ggATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTaa                                                                                                  <  1:918750/100‑1 (MQ=255)
   gATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAAt                                                                                                 <  1:1978719/100‑1 (MQ=255)
    aTTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATc                                                                                                <  2:681132/100‑1 (MQ=255)
       gTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGt                                                                                             >  2:317380/1‑100 (MQ=255)
       gTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGt                                                                                             >  1:602702/1‑100 (MQ=255)
          ccGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGc                                                                                          <  1:1564500/100‑1 (MQ=255)
           cGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGcc                                                                                         <  2:923137/100‑1 (MQ=255)
            ggCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGcc                                                                                         <  1:2431336/99‑1 (MQ=255)
             gatAAAGAAAAAATAGGGTGTGTTGGCGGCCTGCAAGTATTGTAAGACCCGTTAAGGTTTAAAAAAAAAAAAAAAAGGCGACAGATTACTCTGTCGCCt                                                                                        <  2:125565/97‑1 (MQ=16)
             gCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAATGGTTCTGAATCAGGAAAAAAGGCGACAGATTACTCTGTCGcc                                                                                         <  2:1607152/98‑1 (MQ=38)
              cTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCttt                                                                                      >  2:1665733/1‑100 (MQ=255)
               tAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAAACGACAGAGTAATCTGTCGCCtttt                                                                                     >  1:1289267/1‑100 (MQ=255)
                 aaGGAAAAACAGGGTGTGTTGGGGGCCTGCAAGGATTGTAAGACCCGTTAAGCGTTATGAGTCAGGATAAAAGGCGACAGATTACTCTGTCGCCtttttt                                                                                   <  2:2755136/100‑1 (MQ=38)
                  aGGAACAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTc                                                                                  <  2:2282357/100‑1 (MQ=255)
                  aGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAAGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTc                                                                                  >  1:801544/1‑100 (MQ=255)
                   ggAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCt                                                                                 >  2:1286342/1‑100 (MQ=255)
                   ggAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCt                                                                                 >  1:1225154/1‑100 (MQ=255)
                   ggAAAAACAGGGTGTGTTGGCGGCCTGCAAGGAATGTAAGACCCGTTAAGGGTTATGCGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCt                                                                                 >  2:889075/1‑100 (MQ=255)
                   ggAAAAACAGGGTGTGTTGGCGGCATGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTGAGGAAAAAAGGAGAAAGAGTAATCTGTCGCCTTTTTTCt                                                                                 >  2:491055/1‑100 (MQ=255)
                    gAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCtt                                                                                >  2:1516476/1‑100 (MQ=255)
                    gAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCtt                                                                                <  2:1867969/100‑1 (MQ=255)
                    gAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAAGCGACAGAGTAATCTGTCGCCTTTTTTCtt                                                                                >  1:1580468/1‑100 (MQ=255)
                     aaaaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCttt                                                                               >  1:540338/1‑100 (MQ=255)
                      aaaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttg                                                                              >  1:584654/1‑100 (MQ=255)
                       aaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCATTTTTTCTttgc                                                                             >  2:1147156/1‑100 (MQ=255)
                        aagaGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgct                                                                            >  2:2153620/4‑100 (MQ=255)
                           aGGGTGTGTTGGCGGCTTGCAAGTATTGTAAGATCCGTTAAGGGTTATGTGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgcttg                                                                          <  1:1602616/99‑1 (MQ=38)
                           aGGGTGTGTTGGCGGCCTGCAAGTATTGTAAGATCCGTTAAGGTTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgcttgc                                                                         <  1:1696228/100‑1 (MQ=255)
                            gggTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgcttgc                                                                         >  2:2947262/1‑99 (MQ=255)
                            gggTGTGTGGGAGGCCTGCAAGGATTGTAAGAGCCGTTATGGGTTATGAGTTAGGAAAAAAGGCGACAGTGTAATCTGTCGCCGTTTTTCTttgct                                                                            <  2:2615595/96‑1 (MQ=255)
                             ggTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGGAAAAAGGCGACAGATTATTCTGTCGCCTTTTTTcctgactcat                                                                         <  2:2319573/98‑10 (MQ=25)
                               tgtgtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTc                                                                     <  2:1709421/100‑1 (MQ=255)
                                gtgtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCt                                                                    >  1:1682019/1‑100 (MQ=255)
                                 tgttgtCGGCCTGCAAGGATTGTAATACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTg                                                                   <  2:2559/100‑1 (MQ=255)
                                 tgtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTg                                                                   <  1:2277518/100‑1 (MQ=255)
                                   ttGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGtt                                                                 <  2:1772206/100‑1 (MQ=255)
                                    tGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTa                                                                <  2:2180907/100‑1 (MQ=255)
                                    tGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTa                                                                >  2:338158/1‑100 (MQ=255)
                                    tGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTa                                                                >  2:1148301/1‑100 (MQ=255)
                                     ggcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAg                                                               <  2:2321540/100‑1 (MQ=255)
                                     ggcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAg                                                               >  2:857205/1‑100 (MQ=255)
                                     ggcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAg                                                               <  1:2839754/100‑1 (MQ=255)
                                      gcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGa                                                              <  1:436847/100‑1 (MQ=255)
                                      gcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGa                                                              >  1:2323402/1‑100 (MQ=255)
                                             gCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgcc                                                       >  1:2979595/1‑100 (MQ=255)
                                              cAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgccg                                                      >  1:2494205/1‑100 (MQ=255)
                                               aaGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCcgccgc                                                     >  1:1631512/1‑100 (MQ=255)
                                                aGGATTGTAAGACCTTTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCa                                                    <  1:1389982/100‑1 (MQ=255)
                                                aGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGATTGCTTTCTGTTAGATTCCGCCGCa                                                    >  1:1714569/1‑100 (MQ=255)
                                                   aTTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTACGCCGCaaa                                                  >  2:2766516/1‑99 (MQ=255)
                                                    ttGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATc                                                <  1:1865684/100‑1 (MQ=255)
                                                      gTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCta                                              >  2:1337334/1‑100 (MQ=255)
                                                      gTAAGACCAGTTAATGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTAGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCta                                              >  2:356362/1‑100 (MQ=255)
                                                       tatgtcccTTTATGGGTTATTGCTCATGAAAAAAGGCGACATACTAATCTGTCGCTTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATc                                                <  2:2596573/92‑1 (MQ=255)
                                                           aCCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTg                                         <  2:269164/100‑1 (MQ=255)
                                                           aCCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTg                                         <  2:1409697/100‑1 (MQ=255)
                                                            cccGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGa                                        >  1:2937814/1‑100 (MQ=255)
                                                            cccGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGa                                        >  1:2704990/1‑100 (MQ=255)
                                                            cccGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATACCGCCGCAAATCTATACTGa                                        >  1:1196429/1‑100 (MQ=255)
                                                               gtgaggggTTATGAGTCAGTAAAAGAGGCGACAGAGTAATCTGTCGGGTTTTTTCTTTGCTTGCGTTGTGTTATATTCTGCCGTAAATCTATACTGAATa                                     <  2:2108792/95‑1 (MQ=255)
                                                                   aGGGTTATGAGTCAGGAAAAAAAGCGACAGAGAAATCTGACGACTTTTTTCTTTGCTTGCTATATGTTAGATTAAGCCGCAAATCTATATTGAATaaaa                                  >  2:810084/1‑99 (MQ=255)
                                                                      gTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCGAATCTATACTGAATAAAACGGc                              <  1:432060/100‑1 (MQ=255)
                                                                      gTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGc                              >  2:2460773/1‑100 (MQ=255)
                                                                           gAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTaaaa                         <  1:83896/100‑1 (MQ=255)
                                                                            aGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAg                        >  1:2319646/1‑100 (MQ=255)
                                                                             gTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGa                       >  1:2765103/1‑100 (MQ=255)
                                                                             gTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGa                       <  1:1886969/100‑1 (MQ=255)
                                                                              tCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGAc                      <  1:913312/100‑1 (MQ=255)
                                                                               cAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACg                     >  1:2215209/1‑100 (MQ=255)
                                                                               cAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACg                     >  2:482114/1‑100 (MQ=255)
                                                                                aagAAAAGAGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAAACGa                    >  1:1588673/3‑100 (MQ=37)
                                                                                aGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGa                    >  1:1182322/1‑100 (MQ=255)
                                                                                aGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGa                    >  2:1140507/1‑100 (MQ=255)
                                                                                 ggAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACgta                   <  1:2046048/100‑3 (MQ=255)
                                                                                  gAAAAAAGGCGACAGATTACTCTGTTGCCTTTTTTATTTGCGTGCGTTCTGTTAGATACCGCCGGGAAGCTATACTGAATAAAACGGCTAAAAAACGa                    >  2:2003111/1‑98 (MQ=21)
                                                                                  gAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGTTTGTTTTTTGTTAGATTTTGCCGAAAATTTATACTGAATAAAACGGCTAAAAGACGAAc                  >  1:2528807/1‑100 (MQ=37)
                                                                                  gAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAAc                  <  2:624233/100‑1 (MQ=255)
                                                                                    aaaaaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCGGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAAcc                 >  2:2700924/1‑99 (MQ=255)
                                                                                    aaaaaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCa                >  1:1319943/1‑100 (MQ=255)
                                                                                    aaaaaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAAACGAACCa                >  1:304452/1‑100 (MQ=38)
                                                                                     aaaaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCAt               >  1:2364828/1‑100 (MQ=255)
                                                                                      naaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAAcc                 >  2:616910/2‑97 (MQ=255)
                                                                                       aaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAAACGAACCATTa             >  2:1340992/1‑100 (MQ=255)
                                                                                        aGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTAt            >  2:2295144/1‑100 (MQ=255)
                                                                                         ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg           >  2:414728/1‑100 (MQ=255)
                                                                                         ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTAAATAAAACGGCTAAAAGACGAACCATTa             >  1:1057508/1‑98 (MQ=38)
                                                                                         ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCAGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg           >  1:802057/1‑100 (MQ=39)
                                                                                         ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTACGCCGCAAATCTATACTGAATAAAAAGGCTAAAAGACGAACCATAATg           >  2:1952149/1‑100 (MQ=38)
                                                                                         ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGAGTCCGCCGCAAATCTATACTGACTAAAACGGCTAAAAGACGAACAACTATg           >  2:578811/1‑100 (MQ=37)
                                                                                         ggCGACAGATTACTCTGTCACCAATTTTCGTTGCTTGCTTTCTGATAGATTCGCGCGCCAATCTATACTGACAAAAAAGGCTGAAAGACGAACCATTATg           >  2:1328701/1‑100 (MQ=255)
                                                                                         ggCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg           <  1:1771444/100‑1 (MQ=255)
                                                                                           cGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATAATGAATAAAACGGCTAAAAGACGAACCATTATGaa         >  1:2662200/1‑100 (MQ=255)
                                                                                            gACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGaaa        >  2:2489994/1‑100 (MQ=255)
                                                                                              cAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAGCCATTATGAAAcc      <  1:927727/100‑1 (MQ=255)
                                                                                               agagTAATCTGTCGCCTTTTTTATTTGATTGCTTTCTGTTAGATTACGCCGCAAATCTATACTGAATAAAACGGATAAAAGACGAACCATTATGAAACca     >  2:2767810/1‑100 (MQ=255)
                                                                                                gagTAATCTGTCGCCTTTTTTTTTTTCTTTCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACca     >  2:566602/1‑99 (MQ=255)
                                                                                                gagTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcat    <  1:2439430/100‑1 (MQ=255)
                                                                                                gagTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACca     >  2:547440/1‑99 (MQ=255)
                                                                                                  gTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcatca  >  1:960907/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
ACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCA  >  NC_000913/4051161‑4051358

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