Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,051,259 G→T 21.2% 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.2% 283.8 / 65.0 123intergenic (‑494/‑88)yihA/yihIcell division GTP‑binding protein/activator of Der GTPase
Reads supporting (aligned to +/- strand):  ref base G (51/45);  new base T (17/9);  total (69/54)
Fisher's exact test for biased strand distribution p-value = 3.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

Only 100 of 134 total aligned reads displayed.

CACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCA  >  NC_000913/4051160‑4051358
                                                                                                   |                                                                                                   
cacGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACagag                                                                                                     >  1:1519316/1‑100 (MQ=255)
cacGGATGGGTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGAAAGGATTGTAAGACCCGTTAAGGGTTAAAAAAAAAAAAAAAAGGCGACagag                                                                                                     <  1:2422894/100‑1 (MQ=255)
  cagATTGTTCCGGCTAAAGGTAAAACAGTGTTTGTTGGCGGCCTGCAAGGATTTTAAGGACCGTTAAGGGTTATGAGTCAGGAACAAAGGCGACagag                                                                                                     <  1:1597008/96‑1 (MQ=255)
  cGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTa                                                                                                   >  1:2572975/1‑100 (MQ=255)
    ggtggtagcGGCTAAAGGAAAAACAGGGTGTGATGGCGGCTTGCAAGGATTGTAAGACCCGTTAAAGGATAAAAAAAAAAAAAAAAGGCGACAGAGTa                                                                                                   <  1:1480886/90‑1 (MQ=255)
    gATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAGGCGACAGAGTAAt                                                                                                 >  1:2268143/1‑99 (MQ=255)
     aTTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATc                                                                                                <  1:1313633/100‑1 (MQ=255)
        gTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGt                                                                                             >  2:255919/1‑100 (MQ=255)
            cGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGcc                                                                                         <  2:331175/100‑1 (MQ=255)
            cGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGcc                                                                                         >  1:1069307/1‑100 (MQ=255)
            cGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGcc                                                                                         <  2:1069311/100‑1 (MQ=255)
            cGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGcc                                                                                         <  2:1172987/100‑1 (MQ=255)
             ggCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGCCt                                                                                        <  1:868585/100‑1 (MQ=255)
             gacTAAAGGAAAAACAGGGTGTGTTGGCGTCAAGCTAGGATTGTAAGACCAGATAAGGGTTATGAGTTAGAAAAAAAGGAGCCAGAGTAATCTGTCgtgt                                                                                        >  1:1647510/3‑97 (MQ=255)
              gatAAAGCAAAAACAGGGTGTGCTGGCGGCCTGCAAAGATTGTAAGCCCCGTTAAGCGTTATAAATAAAGAAAAAAGGCGACAGATTACTCTGTCGCCt                                                                                        <  1:4280/97‑1 (MQ=18)
              gCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGCCtt                                                                                       <  2:2246789/100‑1 (MQ=255)
               cTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGCCttt                                                                                      <  2:2341454/100‑1 (MQ=255)
                tAAAGGAAAAACATGGTGTGTTGGCGGGCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGcc                                                                                         <  1:593728/96‑1 (MQ=38)
                  aaGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCtttttt                                                                                   >  1:804301/1‑100 (MQ=255)
                    ggAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCt                                                                                 >  1:1058397/1‑100 (MQ=255)
                    ggAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCt                                                                                 >  2:2584430/1‑100 (MQ=255)
                      aaaaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCttt                                                                               >  2:1324960/1‑100 (MQ=255)
                       aaaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttg                                                                              <  2:353410/100‑1 (MQ=255)
                       aaaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttg                                                                              >  2:341568/1‑100 (MQ=255)
                       aaaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttg                                                                              >  2:2107687/1‑100 (MQ=255)
                        aaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAGGGCGACAGAGTAATCTGTCGCCTTTTTTCTttg                                                                              <  1:1392842/99‑1 (MQ=255)
                         aaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgct                                                                            <  2:1799945/100‑1 (MQ=255)
                          ncAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgctt                                                                           >  1:928606/2‑100 (MQ=255)
                              ggTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTgcttgctt                                                                       <  1:2527621/100‑1 (MQ=255)
                               gtgtgtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCttt                                                                      <  1:2415168/100‑1 (MQ=255)
                                   nttGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGt                                                                  >  2:142526/2‑100 (MQ=255)
                                    ttGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTTAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGt                                                                  <  1:965930/99‑1 (MQ=255)
                                     tGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGCCAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTa                                                                >  1:2089671/1‑100 (MQ=255)
                                     tGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTa                                                                >  2:2419022/1‑100 (MQ=255)
                                      ggcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAg                                                               >  1:1858973/1‑100 (MQ=255)
                                      ggcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAg                                                               <  2:1519325/100‑1 (MQ=255)
                                        cggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGAt                                                             >  1:1364558/1‑100 (MQ=255)
                                           ccTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTcc                                                          <  2:509199/100‑1 (MQ=255)
                                            cTGCAAGGATTGTGAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccg                                                         >  2:1741018/1‑100 (MQ=255)
                                            cTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccg                                                         >  1:1839998/1‑100 (MQ=255)
                                              gCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgcc                                                       >  2:1857511/1‑100 (MQ=255)
                                               cAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgccg                                                      <  2:173543/100‑1 (MQ=255)
                                                aaGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCcgccgc                                                     <  2:2479925/100‑1 (MQ=255)
                                                 aGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCcgccgc                                                     >  1:1895630/1‑99 (MQ=255)
                                                 aGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCa                                                    <  2:1992876/100‑1 (MQ=255)
                                                 aGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCa                                                    >  1:285845/1‑100 (MQ=255)
                                                 aGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCa                                                    <  1:2314310/100‑1 (MQ=255)
                                                 aGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCa                                                    >  2:619889/1‑100 (MQ=255)
                                                   gATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCaaa                                                  >  1:1116845/1‑100 (MQ=255)
                                                   gATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCaaa                                                  <  2:21074/100‑1 (MQ=255)
                                                    aTTGTAAGACCCGTTAAGGGTTATTAGTCAGGATAAAAGGCGACTTAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGAt                                                             <  1:629024/88‑1 (MQ=255)
                                                       gTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCta                                              >  2:752509/1‑100 (MQ=255)
                                                       gTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCATTTTTTCTTTGCTTGCTGTCTGTTAGATTCCGCCGCAAATCta                                              >  1:1013403/1‑100 (MQ=255)
                                                          aGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATAc                                           <  2:2187464/100‑1 (MQ=255)
                                                          aGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTACGCCGCAAATCTATAc                                           >  1:1750749/1‑100 (MQ=255)
                                                           gACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACt                                          <  1:703361/100‑1 (MQ=255)
                                                             cccGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGa                                        >  1:1548497/1‑100 (MQ=255)
                                                              ccGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGaa                                       <  2:1116849/100‑1 (MQ=255)
                                                                    aGGGTTATGAGTCAGGTAAAAAGGCGACAGTGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAAc                                 <  1:1239769/100‑1 (MQ=255)
                                                                    aGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAAc                                 <  2:456485/100‑1 (MQ=255)
                                                                       gTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGc                              <  2:2070038/100‑1 (MQ=255)
                                                                       gTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGc                              <  2:2089685/100‑1 (MQ=255)
                                                                        ttATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCt                             <  1:237429/100‑1 (MQ=255)
                                                                         tATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTa                            >  2:620752/1‑100 (MQ=255)
                                                                         tATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTa                            <  2:1546541/100‑1 (MQ=255)
                                                                         tATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTa                            <  2:1587108/100‑1 (MQ=255)
                                                                           tGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTaaa                          >  1:699611/1‑100 (MQ=255)
                                                                            gAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTaaa                          >  1:1018996/1‑99 (MQ=255)
                                                                              gTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGa                       >  1:1675570/1‑100 (MQ=255)
                                                                                cAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACg                     <  2:786163/100‑1 (MQ=255)
                                                                                cAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTCAGCTTCCGCAGCATATCTATACTGAATAAAACGGCTAAAACACg                     >  1:2144717/1‑100 (MQ=38)
                                                                                 aagAAAAAAGGCGACAGATTACTCTGTCGCCTTTTTTCATTGCTTGCTGTCTGCTAGATTACGCCACTAATCTATACTGAATAAAACGGCTAAGAGACGa                    >  1:1018602/3‑100 (MQ=21)
                                                                                   gacaaaaGGCGGCCGAGTGATCTGAAGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCAGCCGCAAATCTATACTGAGTAAAACGGCTAAAAGACGaa                   <  1:1751775/96‑1 (MQ=255)
                                                                                     aaaaaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCa                >  2:514543/1‑100 (MQ=255)
                                                                                     aaaaaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCa                >  2:1035345/1‑100 (MQ=255)
                                                                                     aaaaaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCa                >  1:1838752/1‑100 (MQ=255)
                                                                                        ctggCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTa             <  2:680803/98‑1 (MQ=255)
                                                                                        aaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTCGATTCCGCCGCAACTCTATACTGAATAACACGGCTAAAAGAGGAACCATTa             >  1:8384/1‑100 (MQ=37)
                                                                                        aaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTa             >  2:185594/1‑100 (MQ=255)
                                                                                        aaGGCGACAGAGTAATCTGTCACCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTa             <  2:1230561/100‑1 (MQ=255)
                                                                                         aGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTAt            >  1:1921772/1‑100 (MQ=255)
                                                                                          ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGCTTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGAAGAACCATTAt            >  2:710426/1‑99 (MQ=38)
                                                                                          ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTCAAAGACGAACCATTAt            >  1:524433/1‑99 (MQ=255)
                                                                                          ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg           >  2:1532117/1‑100 (MQ=255)
                                                                                          ggCGACAGAGTAATCTGTCGCCTTTTTTTTTTTTTTTCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTaa            >  1:680800/1‑98 (MQ=255)
                                                                                          ggCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg           <  2:606652/100‑1 (MQ=255)
                                                                                          ggCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg           >  1:1969919/1‑100 (MQ=255)
                                                                                           gCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg           >  1:1662758/1‑99 (MQ=255)
                                                                                            cGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGGCGAACCATTATGaa         >  2:2306165/1‑100 (MQ=255)
                                                                                            cGACAGAGTAATCTGTCGCCTTTTGTCTACGCGTGCATTCTGTTAGATTCCGGCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGaa         <  1:1498332/100‑1 (MQ=255)
                                                                                             gACAGAGTAATCTGTCGCCTTTTTTTTTTTTTTTCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGaaa        >  1:1266197/1‑100 (MQ=255)
                                                                                             gACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGaaa        <  2:1320207/100‑1 (MQ=255)
                                                                                              aCAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAAc       <  2:699614/100‑1 (MQ=255)
                                                                                                agagTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACca     <  2:1609087/100‑1 (MQ=255)
                                                                                                 gagTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcat    >  1:2153491/1‑100 (MQ=255)
                                                                                                 gagTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcat    >  2:2538350/1‑100 (MQ=255)
                                                                                                 gagTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcat    >  2:1032399/1‑100 (MQ=255)
                                                                                                 gagTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcat    >  2:2353756/1‑100 (MQ=255)
                                                                                                 gagTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcat    >  2:1813938/1‑100 (MQ=255)
                                                                                                   gTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcatca  <  2:164370/100‑1 (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: