Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,051,262 A→C 18.1% intergenic (‑497/‑85) 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,2620AC18.1% 272.5 / 57.8 122intergenic (‑497/‑85)yihA/yihIcell division GTP‑binding protein/activator of Der GTPase
Reads supporting (aligned to +/- strand):  ref base A (56/43);  new base C (17/5);  total (73/49)
Fisher's exact test for biased strand distribution p-value = 9.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.89e-01

Only 100 of 126 total aligned reads displayed.

GGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCT  >  NC_000913/4051163‑4051361
                                                                                                   |                                                                                                   
 gATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAAt                                                                                                    <  1:1978719/100‑1 (MQ=255)
  aTTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATc                                                                                                   <  1:2063636/100‑1 (MQ=255)
  aTTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATc                                                                                                   <  2:681132/100‑1 (MQ=255)
     gTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGt                                                                                                >  1:2350258/1‑100 (MQ=255)
     gTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGt                                                                                                >  1:602702/1‑100 (MQ=255)
        ccGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGc                                                                                             <  2:578330/100‑1 (MQ=255)
        ccGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGc                                                                                             <  1:1564500/100‑1 (MQ=255)
         cGGCTAAAGGAAAAACAGGGTGTGTTGGCGGTCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGc                                                                                             <  1:594224/99‑1 (MQ=39)
         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)
               aaGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGACGACAGAGTAATCTGTCGCCtttttt                                                                                      >  1:2801811/1‑100 (MQ=255)
                aGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAAGCAACAGAGTAATCTGTCACCtttttt                                                                                      >  1:981443/1‑99 (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:1867969/100‑1 (MQ=255)
                  gAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCtt                                                                                   >  2:1516476/1‑100 (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)
                              gtgtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCt                                                                       >  1:1682019/1‑100 (MQ=255)
                               tgttgtCGGCCTGCAAGGATTGTAATACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTg                                                                      <  2:2559/100‑1 (MQ=255)
                               tgtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTg                                                                      <  1:2277518/100‑1 (MQ=255)
                                 ttGGCGGCCTGCAAGGATTGTAAGACCCGTTAATGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGtt                                                                    <  1:899237/100‑1 (MQ=255)
                                 ttGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGtt                                                                    <  2:1772206/100‑1 (MQ=255)
                                  tGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTa                                                                   >  2:338158/1‑100 (MQ=255)
                                  tGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTa                                                                   <  2:2180907/100‑1 (MQ=255)
                                  tGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTa                                                                   >  2:1148301/1‑100 (MQ=255)
                                   ggcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAg                                                                  <  1:2839754/100‑1 (MQ=255)
                                   ggcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAg                                                                  <  2:2321540/100‑1 (MQ=255)
                                   ggcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAg                                                                  >  2:857205/1‑100 (MQ=255)
                                    gcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGa                                                                 <  1:436847/100‑1 (MQ=255)
                                    gcggcCGGCAAAGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGa                                                                 <  2:3026010/100‑1 (MQ=255)
                                     cggcCTGCAAGGATTGTAAAACACGTTAAGGGTTATGAGACAGGAAAAAAGGCAACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTa                                                                   >  2:749059/1‑97 (MQ=39)
                                            cAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgccg                                                         >  1:2494205/1‑100 (MQ=255)
                                             aaGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCcgccgc                                                        >  1:1631512/1‑100 (MQ=255)
                                              aGGATTGTAAGACCTTTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCa                                                       <  1:1389982/100‑1 (MQ=255)
                                              aGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCa                                                       >  2:198453/1‑100 (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)
                                                         aCCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTg                                            <  2:269164/100‑1 (MQ=255)
                                                         aCCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTCTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTg                                            <  2:1681999/100‑1 (MQ=255)
                                                          cccGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGa                                           >  1:2937814/1‑100 (MQ=255)
                                                          cccGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATACCGCCGCAAATCTATACTGa                                           >  1:1196429/1‑100 (MQ=255)
                                                             gtgaggggTTATGAGTCAGTAAAAGAGGCGACAGAGTAATCTGTCGGGTTTTTTCTTTGCTTGCGTTGTGTTATATTCTGCCGTAAATCTATACTGAATa                                        <  2:2108792/95‑1 (MQ=255)
                                                                    gTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCGAATCTATACTGAATAAAACGGc                                 <  1:432060/100‑1 (MQ=255)
                                                                    gTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGc                                 >  2:2460773/1‑100 (MQ=255)
                                                                     ttATGAGTCAGGAAAGAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCt                                <  2:2741201/100‑1 (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)
                                                                              aagAAAAGAGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAAACGa                       >  1:1588673/3‑100 (MQ=37)
                                                                              aGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGa                       >  1:1182322/1‑100 (MQ=255)
                                                                                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)
                                                                                    aaaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCAtt                 >  2:1343027/1‑100 (MQ=255)
                                                                                     aaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAAACGAACCATTa                >  2:1340992/1‑100 (MQ=255)
                                                                                      aGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTa                >  1:330520/1‑99 (MQ=255)
                                                                                      aGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTAt               >  2:2295144/1‑100 (MQ=255)
                                                                                       ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg              >  1:2516658/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)
                                                                                       ggCGACAGATTACTCTGTCACCAATTTTCGTTGCTTGCTTTCTGATAGATTCGCGCGCCAATCTATACTGACAAAAAAGGCTGAAAGACGAACCATTATg              >  2:1328701/1‑100 (MQ=255)
                                                                                       ggCGACAGATTACTCTGTAGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCAGCAAATCTATACTGAATAAAACGGCTAAAAGACAAACCATTATg              >  2:1255748/1‑100 (MQ=38)
                                                                                       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)
                                                                                              gagTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACca        >  2:547440/1‑99 (MQ=255)
                                                                                                gTAATCTGTCGCCTTTTTTCTTTTCTTTCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAAACGAACCATTATGAAACcat       >  2:455851/1‑98 (MQ=255)
                                                                                                gTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcatca     >  1:960907/1‑100 (MQ=255)
                                                                                                gTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcatca     <  2:2185202/100‑1 (MQ=255)
                                                                                                   aTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCt  <  2:2356632/100‑1 (MQ=255)
                                                                                                   aTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCt  >  1:1942755/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCT  >  NC_000913/4051163‑4051361

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: