Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,051,262 A→C 18.5% 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.5% 414.0 / 94.6 185intergenic (‑497/‑85)yihA/yihIcell division GTP‑binding protein/activator of Der GTPase
Reads supporting (aligned to +/- strand):  ref base A (91/59);  new base C (23/11);  total (114/71)
Fisher's exact test for biased strand distribution p-value = 5.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.63e-01

Only 100 of 190 total aligned reads displayed.

GGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCT  >  NC_000913/4051163‑4051361
                                                                                                   |                                                                                                   
ggATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTaa                                                                                                     >  1:3613391/1‑100 (MQ=255)
 gATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAAt                                                                                                    >  1:3005622/1‑100 (MQ=255)
 gATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAAt                                                                                                    >  2:2018540/1‑100 (MQ=255)
 gATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGACGACAGAGTaa                                                                                                     >  1:2850754/1‑99 (MQ=255)
  aTTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATc                                                                                                   <  1:3642098/100‑1 (MQ=255)
  aTTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATAAAAAAAAAAAAAAGGCGACAGAGTAATc                                                                                                   <  2:3050661/100‑1 (MQ=255)
   ttGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCt                                                                                                  <  2:49723/100‑1 (MQ=255)
    tGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTg                                                                                                 <  1:3925080/100‑1 (MQ=255)
    tGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGAAAAAAAAACGACAGAGTAATCt                                                                                                  >  1:787601/1‑99 (MQ=255)
      ttCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTc                                                                                               >  1:1420973/1‑100 (MQ=255)
       tCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCg                                                                                              <  2:1515273/100‑1 (MQ=255)
         cGGCTAAATGAAAAACAGGGCGTGTTTCCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGcc                                                                                            <  2:1607179/100‑1 (MQ=38)
         cGGCTAAAGGAAAAACAGGGTGTGTTGGTGGCTTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGcc                                                                                            <  1:330577/100‑1 (MQ=38)
         cGGCTAAAGGAAAAACAGGGTGTGTTGGCGGTCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGcc                                                                                            <  1:71158/100‑1 (MQ=39)
         cGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGcc                                                                                            <  1:1915827/100‑1 (MQ=255)
         cGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGcc                                                                                            <  2:1523875/100‑1 (MQ=255)
         cGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGcc                                                                                            <  1:3255534/100‑1 (MQ=255)
          ggCTAAAGGAAAAATAGTGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGcc                                                                                            <  1:1965048/99‑1 (MQ=38)
          ggCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGGTTGTAAGACCCGTTAAGGGTTATGAGTCAGAAAAAAAGGCGACAGAGTAATCTGTCGc                                                                                             <  2:3807525/98‑1 (MQ=255)
          ggCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCt                                                                                           >  2:40884/1‑100 (MQ=255)
          ggCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCt                                                                                           >  1:1283644/1‑100 (MQ=255)
          ggCTAAAGGAAAAACAGGATGTGTTGGCGGCCTGCAAGGATTGTAAGACCTATTAAGAATTAAAAATCAGAAAAAAAGGCGACAGATTACTCTGTCGcc                                                                                            <  2:1817653/99‑1 (MQ=25)
              aaaGGAAAAACAGGGTGTGTTGGCGGCTTGAAAGGATTGTAAGACCCGTTTAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCttn                                                                                         <  2:2906530/98‑2 (MQ=255)
              aaaGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCttttt                                                                                       <  2:3480400/100‑1 (MQ=255)
               aaGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAAGCAACAGAGTAATCTGTCGCCtttttt                                                                                      >  1:3555381/1‑100 (MQ=255)
                aGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGCCTTTTTTc                                                                                     <  1:2675402/100‑1 (MQ=255)
                aGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTc                                                                                     >  2:1355830/1‑100 (MQ=255)
                 ggAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAATCAGGAAAAAAAGCGACAGAGTAATCTGTCGCCTTTTTTCt                                                                                    >  1:324172/1‑100 (MQ=255)
                 ggAAAAACAGGGTGTGTTGGCGGCCTGAAAGGGTTGTAAAACCCGATAAGGGGTATCGGTGAGGGAAAAAGTCGACAGAGTAATCTGTTGGCGTTTTGCt                                                                                    >  2:2276761/1‑100 (MQ=255)
                   aaaaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCttt                                                                                  >  1:2414964/1‑100 (MQ=255)
                    aaaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttg                                                                                 >  1:2181038/1‑100 (MQ=255)
                      aaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgct                                                                               >  1:68431/1‑100 (MQ=255)
                        cAGGGTTTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgc                                                                                <  1:3693851/97‑1 (MQ=255)
                           ggTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTgcttgctt                                                                          >  2:3206436/1‑100 (MQ=255)
                           ggTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTgcttgctt                                                                          >  1:2703314/1‑100 (MQ=255)
                             tgtgtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTc                                                                        <  2:1232200/100‑1 (MQ=255)
                              gtgtTTGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTc                                                                        <  2:1103281/99‑1 (MQ=255)
                                gtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAAGCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGt                                                                     >  2:3388737/1‑100 (MQ=255)
                                  tGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTa                                                                   <  1:949068/100‑1 (MQ=255)
                                   ggcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATTATTCATGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAg                                                                  <  1:265235/100‑1 (MQ=255)
                                   ggcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTa                                                                   >  2:833214/1‑99 (MQ=255)
                                   ggcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAg                                                                  >  2:2719315/1‑100 (MQ=255)
                                   ggcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTACTTTGCTTGCTTTCTGTTAg                                                                  <  2:2388563/100‑1 (MQ=255)
                                   ggcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGACGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAg                                                                  >  1:1050992/1‑100 (MQ=255)
                                    gcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGACAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGa                                                                 >  2:2166455/1‑100 (MQ=255)
                                     cggcCTGCAAGGATTGTAAGACCCGTTAAGGTTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGAt                                                                <  1:1570665/100‑1 (MQ=255)
                                     cggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAAGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGAt                                                                >  1:224404/1‑100 (MQ=255)
                                         cTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccg                                                            >  2:3005878/1‑100 (MQ=255)
                                            cAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgccg                                                         >  1:3230827/1‑100 (MQ=255)
                                               ggATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTACTTGCTTTCTGTTAGATTCCGCCGCaa                                                      >  1:2146878/1‑100 (MQ=255)
                                                  ttGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATc                                                   <  1:1118484/100‑1 (MQ=255)
                                                    gTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCta                                                 >  2:2593435/1‑100 (MQ=255)
                                                    gTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTACTTGCTTTCTGTTAGATTCCGCCGCAAATCta                                                 >  1:3454854/1‑100 (MQ=255)
                                                     tAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCtat                                                <  2:2855766/100‑1 (MQ=255)
                                                      aaGACCCGTTAAGGGTTATGAGTCAGTAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCtata                                               <  2:3277528/100‑1 (MQ=255)
                                                       aGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATAc                                              <  2:589463/100‑1 (MQ=255)
                                                         aCCCGTTAAGGGTTATGAATCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTg                                            >  2:1068170/1‑100 (MQ=255)
                                                          cccGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGa                                           <  2:3230796/100‑1 (MQ=255)
                                                            cGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCACAAATCTATAATGAAt                                         >  1:1792490/1‑100 (MQ=255)
                                                             gTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCTGCAAATCTATACTGAATa                                        >  2:3350615/1‑100 (MQ=255)
                                                                 aGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAAc                                    <  2:331920/100‑1 (MQ=255)
                                                                   ggTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACgg                                  <  2:661670/100‑1 (MQ=255)
                                                                    gTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGc                                 <  1:1838537/100‑1 (MQ=255)
                                                                    gTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTAAATAAAAc                                    >  2:3887317/1‑97 (MQ=255)
                                                                         gAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTaaaa                            <  1:2315251/100‑1 (MQ=255)
                                                                          aGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAg                           >  2:2897707/1‑100 (MQ=255)
                                                                           gTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGa                          <  1:2593459/100‑1 (MQ=255)
                                                                           gTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGa                          >  1:60870/1‑100 (MQ=255)
                                                                             cAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACg                        <  1:2117471/100‑1 (MQ=255)
                                                                             cAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACg                        <  1:3206467/100‑1 (MQ=255)
                                                                               ggAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGaa                      >  1:3216690/1‑100 (MQ=255)
                                                                               ggAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAAACGaa                      >  2:2765319/1‑100 (MQ=255)
                                                                                    aaaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCACATCTATACTGAATAAAACGGCTACAAGACGAACCAtt                 >  2:2584998/1‑100 (MQ=38)
                                                                                    aaaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCAtt                 >  2:803869/1‑100 (MQ=255)
                                                                                    aaaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAAACTATACTGAATAAAACAGCTAAAAGACGAAACAt                  >  1:1741015/1‑99 (MQ=37)
                                                                                     aaGGCGACAGATTACTCTGTCGCCTTTTTTGTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTa                >  2:3590831/1‑100 (MQ=255)
                                                                                     aaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTa                >  1:228113/1‑100 (MQ=255)
                                                                                     aaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTa                <  2:3441518/100‑1 (MQ=255)
                                                                                      aGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTAt               >  1:2914314/1‑100 (MQ=255)
                                                                                      aGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTAt               >  1:573220/1‑100 (MQ=255)
                                                                                      aGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTAt               >  1:1937739/1‑100 (MQ=255)
                                                                                      aGGCGACAGAGTAATCTGACGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACTAACCATTa                >  2:466752/1‑99 (MQ=255)
                                                                                       ggCGACATATTACTCTATCGCATAATTTCTTTGCTTGCTTTCTATTAGATTCCGCCGCAAATCTAAACTAAATAAAACGGCCAAAAGACGAAACATTAt               >  2:3096926/1‑99 (MQ=255)
                                                                                       ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGTCTAAAAGACGAActa                   >  1:3309924/1‑93 (MQ=38)
                                                                                       ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg              >  2:3882780/1‑100 (MQ=255)
                                                                                       ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAAACGAACCATTATg              >  2:1668837/1‑100 (MQ=255)
                                                                                       ggCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAGAGACGAACCATTATg              >  1:3869620/1‑100 (MQ=255)
                                                                                       ggCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg              >  1:3388701/1‑100 (MQ=255)
                                                                                       ggCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg              >  1:3795114/1‑100 (MQ=255)
                                                                                        gCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTACGCCGCAAATCTATACTAAATAAAACGGCTAAAAGACGAACCATTATGa             >  2:1032325/1‑100 (MQ=38)
                                                                                        gCGACAGATCACTCTGTAGCCTTTTATATTTGCGTGCATTCTGTTAAATTCCGGGACAAATCTATACTGAATAAAAAGGATTTAAGACGAAACATTATGa             >  2:1870399/1‑100 (MQ=255)
                                                                                        gCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTAt               >  2:1623650/1‑98 (MQ=255)
                                                                                         cGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGaa            >  1:1986110/1‑100 (MQ=255)
                                                                                            cAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAAcc         >  1:1371322/1‑100 (MQ=255)
                                                                                            cAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAAc          >  1:2478174/1‑99 (MQ=255)
                                                                                            cAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAAACGAACCATTATGAAAc          >  1:1657410/1‑99 (MQ=255)
                                                                                              gagTAATCTGTCGCCTTTTTTCTTTGCTTGTTTTCTGTTAGATTCCGCCGCAAATCTATCGTGAATAAAAAGGCTAAAAGATAAAACATTAAGAAACca        >  2:729409/1‑99 (MQ=255)
                                                                                                 tAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTACATTCCGCCGCAAATATATACTGAATAAAACGACTAAAAGACGAACCATGATGAAACcatcat    >  1:335292/1‑100 (MQ=255)
                                                                                                  aaTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcatcat    >  2:1592057/1‑99 (MQ=255)
                                                                                                   aTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCt  <  1:528000/100‑1 (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: