Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,051,262 A→C 21.2% 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,2620AC21.2% 443.2 / 139.0 222intergenic (‑497/‑85)yihA/yihIcell division GTP‑binding protein/activator of Der GTPase
Reads supporting (aligned to +/- strand):  ref base A (91/83);  new base C (27/20);  total (118/104)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-03

Only 100 of 224 total aligned reads displayed.

GATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCT  >  NC_000913/4051164‑4051361
                                                                                                  |                                                                                                   
gATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAAt                                                                                                    <  2:4203488/100‑1 (MQ=255)
  ttGTTCTGGCTAAAGAAAAAACAGGGTGTCTTGGCGGCATATAAGGATTGTAAGACCTGTTAAGGTTTAAAAAAAAAAAAAAAAGGCGACAGAGTAATCt                                                                                                  <  1:3868236/100‑1 (MQ=255)
  ttGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCt                                                                                                  <  2:158654/100‑1 (MQ=255)
    gTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGt                                                                                                <  1:1132273/100‑1 (MQ=255)
    gTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGt                                                                                                >  2:3643483/1‑100 (MQ=255)
    gTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGt                                                                                                >  2:3807527/1‑100 (MQ=255)
       ccGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGc                                                                                             <  2:134343/100‑1 (MQ=255)
       ccGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGc                                                                                             <  1:2838606/100‑1 (MQ=255)
       ccGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTc                                                                                               <  1:963920/98‑1 (MQ=255)
       ccGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGc                                                                                             <  2:2166058/100‑1 (MQ=255)
        cGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGcc                                                                                            <  2:2549547/100‑1 (MQ=255)
        cGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGcc                                                                                            <  2:916423/100‑1 (MQ=255)
        cGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCg                                                                                              <  1:3984985/98‑1 (MQ=255)
         ggCTAGAGGCAAAACAGGGGGTGTTGGCCGCCTGCAAGGATAGTAAGACCCGCAGAGGTGTATGATTCGGGAAAAAAGGCGACAGAGTAATCTGt                                                                                                >  1:3992130/1‑95 (MQ=255)
         ggCTAAAGTAAAAACAGGGTTTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGACTCAGGAAAAAAGGCGACAGATTACTCTGTCGcc                                                                                            <  1:2756641/99‑1 (MQ=38)
          gCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGCCtt                                                                                          <  2:165731/100‑1 (MQ=255)
          gCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCtt                                                                                          >  1:1345782/1‑100 (MQ=255)
          gCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCtt                                                                                          <  2:1134704/100‑1 (MQ=255)
           cTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGCCttt                                                                                         <  2:1726637/100‑1 (MQ=255)
            taaattaaaaaTATGTTGTGTTGGCGGTCTGCAAGGATTGTAAGATCCGTTAAGGGTTATGAGTTAGGAAAAAAGGCGACAGATTACTCTGTCGCCt                                                                                           <  1:3379394/91‑1 (MQ=25)
            tAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCtttt                                                                                        >  1:1461737/1‑100 (MQ=255)
              aaGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCtttttt                                                                                      <  2:900122/100‑1 (MQ=255)
               aGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGCCTTTTTTc                                                                                     <  2:514435/100‑1 (MQ=255)
               aGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCtttttt                                                                                      >  2:4144202/1‑99 (MQ=255)
               aGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTc                                                                                     >  1:3129547/1‑100 (MQ=255)
               aGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTc                                                                                     >  2:4313788/1‑100 (MQ=255)
                ggAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTACGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGCCTTTTTTcc                                                                                    <  2:1239693/100‑2 (MQ=37)
                ggAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCt                                                                                    >  2:3908267/1‑100 (MQ=255)
                 gAAAAAGAGGGTGTGTTGGGGGCCTGGAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGCCTTTTTTcct                                                                                   <  1:1278502/100‑3 (MQ=37)
                  aaaaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGCCTTTTTTcctg                                                                                  <  2:1547698/100‑4 (MQ=37)
                   aaaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttg                                                                                 >  2:661506/1‑100 (MQ=255)
                   aaaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttg                                                                                 >  1:1042266/1‑100 (MQ=255)
                    aaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTGATGAGTCAGTAAAAAAGGCGACAGGGTGATCTGTCGCCTGTTGTCTttgc                                                                                >  2:3507492/1‑100 (MQ=255)
                           gtgtgtTGGCGGCCTGCAAGGATTGTAGGACCCGTTAAGGTTTATGGGTTAGTAAATAAGGCCACAGAGTAATCTGTCGCTTTTTTTCTTTGCTTGCttt                                                                         <  1:582078/100‑1 (MQ=255)
                            tgtgtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTc                                                                        >  2:2396884/1‑100 (MQ=255)
                            tgtgtTGGCGGCCTGCAAGGACTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTc                                                                        <  2:120793/100‑1 (MQ=255)
                            tgtgtTGGCGCACTTCAAGGACTCTAGGCACCGTCAATTGTTATCAGTCAGGAAAAAAGGCGACAGAGTTTTTTGTCGCCTTTCTTGtt                                                                                   <  1:3782540/89‑1 (MQ=255)
                             gtgtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTATTTGCTTGCTTTCt                                                                       >  1:2147216/1‑100 (MQ=255)
                               gtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGt                                                                     <  1:366502/100‑1 (MQ=255)
                                ttGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGtt                                                                    <  2:2345468/100‑1 (MQ=255)
                                 tGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTa                                                                   <  2:1461742/100‑1 (MQ=255)
                                  ggcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAg                                                                  <  2:3195429/100‑1 (MQ=255)
                                  ggcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAg                                                                  >  2:4103252/1‑100 (MQ=255)
                                   gcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTGATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGa                                                                 >  1:1992847/1‑100 (MQ=255)
                                    cggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGAt                                                                >  2:446900/1‑100 (MQ=255)
                                       ccTGGAAGGATTGTAAGACCCGTTAAGGTTTATGAGTCAGTAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTa                                                                   <  1:4062587/94‑1 (MQ=255)
                                           cAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgccg                                                         >  1:3153941/1‑100 (MQ=255)
                                           cAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgccg                                                         <  1:1667039/100‑1 (MQ=255)
                                           cAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgcc                                                          >  2:116501/1‑99 (MQ=255)
                                            aaGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCcgccgc                                                        <  1:3488361/100‑1 (MQ=255)
                                             aGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCcgccgc                                                        >  1:4137214/1‑99 (MQ=255)
                                             aGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCa                                                       <  2:583366/100‑1 (MQ=255)
                                             aGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCa                                                       >  1:1527124/1‑100 (MQ=255)
                                               gATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACATAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCaaa                                                     >  1:2441142/1‑100 (MQ=255)
                                               gATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCaaa                                                     <  1:4276915/100‑1 (MQ=255)
                                                 ttGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATc                                                   <  2:304379/100‑1 (MQ=255)
                                                  tGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCATTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCt                                                  <  1:3507474/100‑1 (MQ=255)
                                                    tAAGACCCGTTCAGGGTTATGAGTCAGGAAAAACGGCCACACCGTAATCTGCCGCCCCCTTTCTTTGCTCGCTTTCTGTTAGAATCCGCCGCAAATCtat                                                <  1:1296138/100‑1 (MQ=255)
                                                       gACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACt                                             >  2:934296/1‑100 (MQ=255)
                                                        aCCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTg                                            <  2:2369230/100‑1 (MQ=255)
                                                         cccGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGa                                           >  1:541368/1‑100 (MQ=255)
                                                          cctttAAGGGTTATGAGTTAGGAAAAAAGGCGACAGAGTACTCTGTCGCCTTTTTTCTTTGGTTGATTTCTGTTAGATTGCGTCGCAAATCTATAc                                              <  1:3452671/93‑1 (MQ=25)
                                                           cGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAAt                                         <  1:2607362/100‑1 (MQ=255)
                                                           cGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAAt                                         <  1:2383632/100‑1 (MQ=255)
                                                            gTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATa                                        >  2:2800658/1‑100 (MQ=255)
                                                             ttAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATaa                                       <  2:746588/100‑1 (MQ=255)
                                                               aaGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGACTaaaa                                     >  1:91351/1‑100 (MQ=255)
                                                                       tGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCAGCAAATCTATACTGAATAAAACGGCTaaa                             >  1:1775714/1‑100 (MQ=255)
                                                                        gAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTaaaa                            >  2:1746823/1‑100 (MQ=255)
                                                                          gTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGa                          >  2:1854154/1‑100 (MQ=255)
                                                                          gTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGa                          <  1:831793/100‑1 (MQ=255)
                                                                          gTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGa                          <  2:3535771/100‑1 (MQ=255)
                                                                            cAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACg                        >  1:1654029/1‑100 (MQ=255)
                                                                             aGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGa                       <  2:861663/100‑1 (MQ=255)
                                                                                 aaaaaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCa                   >  1:4222009/1‑100 (MQ=255)
                                                                                 aaaaaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCa                   >  2:3605727/1‑100 (MQ=255)
                                                                                 aaaaaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCa                   >  2:201399/1‑100 (MQ=255)
                                                                                   aaaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCAtt                 >  2:697241/1‑100 (MQ=255)
                                                                                   aaaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTGTCTGTTAGATTCCGCCGCTGATCTATACTGAATATAACGGCTAAAAGACGAAc                     >  1:2963678/1‑96 (MQ=37)
                                                                                    aaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTa                >  2:3520355/1‑100 (MQ=255)
                                                                                    aaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTa                <  2:1993075/100‑1 (MQ=255)
                                                                                     aGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGCATCAAACGGCTAAAAGACGAACCATTAt               >  1:4032832/1‑100 (MQ=38)
                                                                                     aGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTAt               <  2:2828562/100‑1 (MQ=255)
                                                                                      ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATGAAACGGCTAAAAGACGAACCATTAt               >  2:3842629/1‑99 (MQ=38)
                                                                                      ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg              >  2:3369016/1‑100 (MQ=255)
                                                                                      ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg              >  2:2147046/1‑100 (MQ=255)
                                                                                      ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg              >  2:2736076/1‑100 (MQ=255)
                                                                                      ggCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg              >  2:2588326/1‑100 (MQ=255)
                                                                                       gCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACAGATAAAAGACGAACAATTATGa             >  1:2188045/1‑100 (MQ=38)
                                                                                       gCGACAGAGTAATCTGTCGCCTTTTTTCTTTTCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGa             >  1:3264078/1‑100 (MQ=255)
                                                                                        cGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGaa            <  1:229218/100‑1 (MQ=255)
                                                                                         gACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGaaa           <  2:4334013/100‑1 (MQ=255)
                                                                                         gACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGaaa           <  2:13261/100‑1 (MQ=255)
                                                                                          aCAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAAc          <  1:1826820/100‑1 (MQ=255)
                                                                                           cAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAAcc         >  1:191554/1‑100 (MQ=255)
                                                                                           cAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAAcc         >  2:3121902/1‑100 (MQ=255)
                                                                                             gagTAATCTGTCGCCTTTTTTTTTTGTTTTCTTTCTGTTAGATTTCGCCGCAAATTTATTTTGAATAATATGGCTAAAAGTCGAACCTTTTTGTAACcat       >  1:788561/1‑100 (MQ=255)
                                                                                              agTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcatc      <  2:4230404/100‑1 (MQ=255)
                                                                                               gTAATCTGTCGCCTTTTTTTTTTTTTTTCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcatca     >  1:2458781/1‑100 (MQ=255)
                                                                                               gTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcatca     >  2:3724396/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
GATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCT  >  NC_000913/4051164‑4051361

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: