Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,051,262 A→C 20.0% 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,2620AC20.0% 270.0 / 72.8 130intergenic (‑497/‑85)yihA/yihIcell division GTP‑binding protein/activator of Der GTPase
Reads supporting (aligned to +/- strand):  ref base A (56/48);  new base C (17/9);  total (73/57)
Fisher's exact test for biased strand distribution p-value = 3.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

Only 100 of 130 total aligned reads displayed.

GATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCAT  >  NC_000913/4051164‑4051359
                                                                                                  |                                                                                                 
gATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAGGCGACAGAGTAAt                                                                                                  >  1:2268143/1‑99 (MQ=255)
 aTTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATc                                                                                                 <  1:1313633/100‑1 (MQ=255)
  ttGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCt                                                                                                >  1:1663744/1‑100 (MQ=255)
    gTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGt                                                                                              >  2:255919/1‑100 (MQ=255)
    gTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGt                                                                                              >  1:1486056/1‑100 (MQ=255)
    attCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAATAATCTGt                                                                                              >  1:1019393/2‑100 (MQ=255)
        cGGGTAAAGGAAAAACAGGGTGTGTTGGCGGCGTGAAAGGATTGTAAGCCCCGTTAAGGGTTAAAAAAAAGAAAAAAAGGCGACAGATTACTCTGTCGcc                                                                                          <  1:738394/100‑1 (MQ=19)
        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)
         gacTAAAGGAAAAACAGGGTGTGTTGGCGTCAAGCTAGGATTGTAAGACCAGATAAGGGTTATGAGTTAGAAAAAAAGGAGCCAGAGTAATCTGTCgtgt                                                                                         >  1:1647510/3‑97 (MQ=255)
          gatAAAGCAAAAACAGGGTGTGCTGGCGGCCTGCAAAGATTGTAAGCCCCGTTAAGCGTTATAAATAAAGAAAAAAGGCGACAGATTACTCTGTCGCCt                                                                                         <  1:4280/97‑1 (MQ=18)
          gCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGCCtt                                                                                        <  2:2246789/100‑1 (MQ=255)
            tAAAGGAAAAACAGGGTTTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGTGTCTGGAAAAAAGGCGACAGAGTAACCTGACGCCttt                                                                                       <  1:1353867/99‑1 (MQ=255)
              aaGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCtttttt                                                                                    >  1:804301/1‑100 (MQ=255)
                  aaaaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCttt                                                                                >  2:1324960/1‑100 (MQ=255)
                   aaaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttg                                                                               >  2:2107687/1‑100 (MQ=255)
                   aaaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttg                                                                               >  2:341568/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)
                              tgtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTg                                                                    >  2:2495839/1‑100 (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:1247952/100‑1 (MQ=255)
                                 tGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTa                                                                 >  2:2419022/1‑100 (MQ=255)
                                  ggcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAg                                                                >  1:1858973/1‑100 (MQ=255)
                                       ccTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTcc                                                           <  2:509199/100‑1 (MQ=255)
                                        cTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccg                                                          >  1:1839998/1‑100 (MQ=255)
                                          gCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgcc                                                        >  1:2570292/1‑100 (MQ=255)
                                          gCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgcc                                                        >  2:1857511/1‑100 (MQ=255)
                                           cAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgccg                                                       <  2:173543/100‑1 (MQ=255)
                                            aaGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCcgccgc                                                      >  2:2347113/1‑100 (MQ=255)
                                             aGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCcgccgc                                                      >  1:1895630/1‑99 (MQ=255)
                                             aGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCa                                                     <  2:1992876/100‑1 (MQ=255)
                                             aGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCa                                                     >  2:619889/1‑100 (MQ=255)
                                               gATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCaaa                                                   <  2:21074/100‑1 (MQ=255)
                                               gATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCaaa                                                   >  1:1116845/1‑100 (MQ=255)
                                                aTTGTAAGACCCGTTAAGGGTTATTAGTCAGGATAAAAGGCGACTTAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGAt                                                              <  1:629024/88‑1 (MQ=255)
                                                   gTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCta                                               >  2:1734332/1‑100 (MQ=255)
                                                   gTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCta                                               >  2:752509/1‑100 (MQ=255)
                                                   gTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCATTTTTTCTTTGCTTGCTGTCTGTTAGATTCCGCCGCAAATCta                                               >  1:1013403/1‑100 (MQ=255)
                                                      aGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATAc                                            <  1:203836/100‑1 (MQ=255)
                                                      aGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTACGCCGCAAATCTATAc                                            >  1:1750749/1‑100 (MQ=255)
                                                       gACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACt                                           <  1:703361/100‑1 (MQ=255)
                                                         cccGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGa                                         >  2:1902258/1‑100 (MQ=255)
                                                         cccGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGa                                         >  1:1548497/1‑100 (MQ=255)
                                                          ccGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGaa                                        <  2:1116849/100‑1 (MQ=255)
                                                                aGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAAc                                  <  2:2099372/100‑1 (MQ=255)
                                                                aGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAAc                                  <  2:456485/100‑1 (MQ=255)
                                                                aGGGCTATGAGTCAGGAAAAAAGGCGCCAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCAGTTAGATTCCGCCGCAACTCTATACTGAATAAAAc                                  <  1:514176/100‑1 (MQ=39)
                                                                 gggTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACg                                 <  2:1486065/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)
                                                                        gAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTaaa                           >  1:1018996/1‑99 (MQ=255)
                                                                         aGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAg                         <  2:445870/100‑1 (MQ=255)
                                                                          gTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGa                        >  1:224520/1‑100 (MQ=255)
                                                                            cAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACg                      <  2:786163/100‑1 (MQ=255)
                                                                            aaagAAAAAAGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTTTACTGAATAAAACGGCTAAAAGACg                      >  1:1831764/4‑100 (MQ=37)
                                                                             aagAAAAAAGGCGACAGATTACTCTGTCGCCTTTTTTCATTGCTTGCTGTCTGCTAGATTACGCCACTAATCTATACTGAATAAAACGGCTAAGAGACGa                     >  1:1018602/3‑100 (MQ=21)
                                                                             aGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACg                      >  2:1807960/1‑99 (MQ=255)
                                                                               gacaaaaGGCGGCCGAGTGATCTGAAGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCAGCCGCAAATCTATACTGAGTAAAACGGCTAAAAGACGaa                    <  1:1751775/96‑1 (MQ=255)
                                                                               gAAAAAAGGCGACAGAGTAATCTGTCGCCGTTTTTCTATGCTTGCTTGCTGTTAGATTACGCTGTTCATTTATGCTGGACAAAACGATTTAAAGGCGaaa                   >  1:1877867/1‑99 (MQ=16)
                                                                                 aaaaaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCa                 >  2:1035345/1‑100 (MQ=255)
                                                                                  aaaaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCAt                <  1:662555/100‑1 (MQ=255)
                                                                                   aaaCGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATATATACTGAATAAAACGGCTAAAAGACGTACCAtt               >  2:296908/1‑100 (MQ=38)
                                                                                    ctggCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTa              <  2:680803/98‑1 (MQ=255)
                                                                                    aaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTCGATTCCGCCGCAACTCTATACTGAATAACACGGCTAAAAGAGGAACCATTa              >  1:8384/1‑100 (MQ=37)
                                                                                    aaGGCGACAGAGTAATCTGTCACCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTa              <  2:1230561/100‑1 (MQ=255)
                                                                                     aGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTAt             >  2:531765/1‑100 (MQ=255)
                                                                                      ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGCTTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGAAGAACCATTAt             >  2:710426/1‑99 (MQ=38)
                                                                                      ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTCAAAGACGAACCATTAt             >  1:524433/1‑99 (MQ=255)
                                                                                      ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg            >  2:157887/1‑100 (MQ=255)
                                                                                      ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg            >  2:1532117/1‑100 (MQ=255)
                                                                                      ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCATACTGTTAGATTCCGCCGAGAATCTATACTGGATAAAACGGCCACCACACGAACCAtt               >  1:1193533/1‑97 (MQ=21)
                                                                                      ggCGACAGAGTAATCTGTCGCCTTTTTTTTTTTTTTTCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTaa             >  1:680800/1‑98 (MQ=255)
                                                                                      ggCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg            >  2:2239921/1‑100 (MQ=255)
                                                                                      ggCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg            >  1:1969919/1‑100 (MQ=255)
                                                                                      ggCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg            <  2:606652/100‑1 (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         <  1:899964/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:1813938/1‑100 (MQ=255)
                                                                                             gagTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcat     >  2:2538350/1‑100 (MQ=255)
                                                                                             gagTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcat     >  2:2353756/1‑100 (MQ=255)
                                                                                             gagTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcat     >  2:1032399/1‑100 (MQ=255)
                                                                                               gTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcatca   <  2:164370/100‑1 (MQ=255)
                                                                                                tAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcatcat  <  2:2538718/100‑1 (MQ=255)
                                                                                                tAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcatcat  <  2:762523/100‑1 (MQ=255)
                                                                                                  |                                                                                                 
GATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCAT  >  NC_000913/4051164‑4051359

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: