Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,397,557 G→T 100% G49V (GGT→GTT)  mutL → methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,397,5570GT99.0% 337.3 / ‑4.1 106G49V (GGT→GTT) mutLmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (52/52);  total (54/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

Only 100 of 106 total aligned reads displayed.

GGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCG  >  NC_000913/4397460‑4397652
                                                                                                  |                                                                                               
ggtgaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTg                                                                                               <  1:809124/100‑1 (MQ=255)
ggtgaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTg                                                                                               <  2:2298096/100‑1 (MQ=255)
   gaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTgggg                                                                                            <  1:2263387/100‑1 (MQ=255)
   gaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTgggg                                                                                            <  1:796555/100‑1 (MQ=255)
     ggtggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTgggg                                                                                            <  1:1535565/98‑1 (MQ=255)
     ggtggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTgggg                                                                                            <  1:130497/98‑1 (MQ=255)
     ggtggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCg                                                                                          >  2:963683/1‑100 (MQ=255)
     ggtggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCg                                                                                          >  2:36302/1‑100 (MQ=255)
      gtggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGa                                                                                         >  2:133717/1‑100 (MQ=255)
         gtCGGGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGGGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAAc                                                                                      <  1:710459/100‑1 (MQ=255)
          tCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACt                                                                                     >  1:1527424/1‑100 (MQ=255)
           cGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACtt                                                                                    >  1:447096/1‑100 (MQ=255)
           cGAGCGACCTGCGCCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACtt                                                                                    >  2:1739875/1‑100 (MQ=255)
            gAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTa                                                                                   <  1:802959/100‑1 (MQ=255)
             aGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTAt                                                                                  >  1:373123/1‑100 (MQ=255)
             aGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTAt                                                                                  <  1:1059783/100‑1 (MQ=255)
                gACCTTCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCg                                                                               <  1:782583/100‑1 (MQ=255)
                gACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCg                                                                               >  2:452462/1‑100 (MQ=255)
                  ccTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGGGTTGGGGCGAAACTTATCCGCa                                                                             >  2:487652/1‑100 (MQ=255)
                   cTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCAt                                                                            >  2:2193269/1‑100 (MQ=255)
                       ctctatAGCCAAAAAACTAGTGGAAAACAGCCTCGATGCAGGTGCGCCACATAGCGCCATAGATCTCGAACGCGATGGGGCGAACCTTATCCGCGtt                                                                           >  1:2366394/6‑97 (MQ=255)
                        tCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTg                                                                       <  1:781054/100‑1 (MQ=255)
                        tCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTg                                                                       <  2:2075797/100‑1 (MQ=255)
                          ggTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGAt                                                                     >  1:1809975/1‑100 (MQ=255)
                          ggTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTCGGGGCGAAACTTATCCGCATTCGTGAt                                                                     >  2:2513849/1‑100 (MQ=255)
                            tAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATaa                                                                   <  2:1478515/100‑1 (MQ=255)
                             aGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGAGATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAAc                                                                  >  1:1814416/1‑100 (MQ=255)
                              gTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACg                                                                 >  2:1222/1‑100 (MQ=255)
                              gTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACg                                                                 >  1:1995021/1‑100 (MQ=255)
                              gTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACg                                                                 <  2:1435255/100‑1 (MQ=255)
                                cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACgg                                                                >  1:1032759/1‑99 (MQ=255)
                                cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGc                                                               >  2:1887229/1‑100 (MQ=255)
                                cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTAGCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACggn                                                               <  1:2256258/100‑2 (MQ=255)
                                  aaGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTg                                                             <  2:2114434/100‑1 (MQ=255)
                                   aGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGc                                                            >  2:1263089/1‑100 (MQ=255)
                                     aaCTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCTATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCgg                                                          <  2:1334260/100‑1 (MQ=255)
                                     aaCTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCgg                                                          <  1:1342119/100‑1 (MQ=255)
                                          gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATc                                                      >  2:534840/1‑99 (MQ=255)
                                          gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCa                                                     <  2:1776258/100‑1 (MQ=255)
                                          gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCa                                                     >  1:1899601/1‑100 (MQ=255)
                                          gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGAGATTGGTATCGAACGCGTGGGGGCGCAAATTATCAGCACTCGTGATAACGGATGCGGTATc                                                      >  1:727109/1‑99 (MQ=255)
                                           tGGACAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCCTTCGTGATAACGGCTGCGGTATCaa                                                    >  2:2499261/1‑100 (MQ=255)
                                            ggAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaa                                                   >  1:779357/1‑100 (MQ=255)
                                            ggAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaa                                                   >  2:1271638/1‑100 (MQ=255)
                                            ggAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAAATTATCCGCATTAGTGATAAAGGATGCGGTATCa                                                     >  2:476527/1‑98 (MQ=255)
                                             gAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaaa                                                  >  2:286963/1‑100 (MQ=255)
                                                 aCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGa                                              <  2:431653/100‑1 (MQ=255)
                                                   aGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATg                                            <  2:1841154/100‑1 (MQ=255)
                                                   aGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATg                                            <  1:1003228/100‑1 (MQ=255)
                                                   aGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATg                                            <  2:2029156/100‑1 (MQ=255)
                                                   aGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATg                                            >  1:257907/1‑100 (MQ=255)
                                                     ccTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAACAAGATGAg                                          >  2:2331421/1‑100 (MQ=255)
                                                      cTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgc                                         >  2:1041753/1‑100 (MQ=255)
                                                      cTCGATGCAGGTGCGACGCGTATCGATATTGATAACGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAg                                          >  2:2315154/1‑99 (MQ=255)
                                                       tCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgct                                        <  1:1112379/100‑1 (MQ=255)
                                                        agATGCAGGTGCGACGCGTTTGTATAGTGATATCGAACGCGTTGGGGGGAAACTTAACCGCATTCGTGATAGCGGCTGCGGTATCAAAAAAGATg                                            <  1:199568/94‑1 (MQ=255)
                                                         gATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctgg                                      <  2:1307512/100‑1 (MQ=255)
                                                           tGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcg                                    <  1:1220304/100‑1 (MQ=255)
                                                             cAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgct                                  >  1:2517055/1‑100 (MQ=255)
                                                              aGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctg                                 <  1:1264182/100‑1 (MQ=255)
                                                               ggTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctgg                                >  1:1560634/1‑100 (MQ=255)
                                                               ggTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctgg                                <  1:1855949/100‑1 (MQ=255)
                                                               ggTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctgg                                <  2:1888122/100‑1 (MQ=255)
                                                                gTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggc                               <  2:569208/100‑1 (MQ=255)
                                                                gTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggc                               <  1:384880/100‑1 (MQ=255)
                                                                ccccaaCGCGTTCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggc                               >  2:2291356/6‑100 (MQ=25)
                                                                ccccaaCGCGTTCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggc                               >  1:718657/6‑100 (MQ=25)
                                                                ccccaaCGCGTTCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggc                               >  2:1652401/6‑100 (MQ=25)
                                                                 tGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcg                              <  1:2054742/100‑1 (MQ=255)
                                                                 tGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcg                              <  2:177819/100‑1 (MQ=255)
                                                                   cGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgct                            <  1:1263085/100‑1 (MQ=255)
                                                                    gACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctn                           <  1:1047308/100‑2 (MQ=255)
                                                                    gACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctg                           >  1:1035826/1‑100 (MQ=255)
                                                                        gcgttCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCt                        >  1:1961647/5‑100 (MQ=37)
                                                                        cgTATCGATATTGATCTCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTc                       <  2:2550557/100‑1 (MQ=255)
                                                                        cgTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTc                       <  2:1911630/100‑1 (MQ=255)
                                                                        cgTATCGATATTGATATCAAAACCGTTGGGGCGAAAATTATCCGCATTCGTGATAACAGCTCCGGTATCAAAAAAGATAAgctggcgcgggcgc                             >  1:2437491/1‑94 (MQ=255)
                                                                         gTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCg                      <  2:843438/100‑1 (MQ=255)
                                                                         gTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCg                      <  2:1350249/100‑1 (MQ=255)
                                                                         cgttCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTc                       >  2:317315/4‑100 (MQ=38)
                                                                          tATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCg                      >  1:25095/1‑99 (MQ=255)
                                                                          tATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGt                     >  2:1518507/1‑100 (MQ=255)
                                                                            tCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCa                   <  1:36301/100‑1 (MQ=255)
                                                                            tCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCa                   <  2:1546234/100‑1 (MQ=255)
                                                                             cGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGCATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCAt                  <  2:2278054/100‑1 (MQ=255)
                                                                                tatTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcc               <  2:2377852/100‑1 (MQ=255)
                                                                                  ttGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccac             >  2:891879/1‑100 (MQ=255)
                                                                                  ttGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccac             >  1:71342/1‑100 (MQ=255)
                                                                                  ttGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccac             >  2:1564876/1‑100 (MQ=255)
                                                                                   tGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacc            <  2:1017836/100‑1 (MQ=255)
                                                                                   tGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacc            <  2:866766/100‑1 (MQ=255)
                                                                                   tGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacc            >  1:828808/1‑100 (MQ=255)
                                                                                   tGATATCGAACGCGTTGGGGCGAAACTTATCCGCATCCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacc            <  1:954810/100‑1 (MQ=255)
                                                                                     atatCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAg          >  1:1276184/1‑100 (MQ=255)
                                                                                      tatCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAg          <  2:1469496/99‑1 (MQ=255)
                                                                                      tatCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGt         >  1:2015033/1‑100 (MQ=255)
                                                                                      tatCGAACGCGTTGGAGCGAAACTTATCCGTATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGCCTCTGGCGCTGGCTCGTCATGCCACCAg          <  1:1171801/99‑1 (MQ=255)
                                                                                        tCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaa       >  1:2233522/1‑100 (MQ=255)
                                                                                         cGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTGGTAATGCCACCAGTaaa      >  1:1230411/1‑100 (MQ=255)
                                                                                              cgcgTTGGGGCGAAACTTATCCGCATTCGTGATAACTGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCg  <  2:1814426/100‑1 (MQ=255)
                                                                                                  |                                                                                               
GGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCG  >  NC_000913/4397460‑4397652

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: