Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,397,557 G→T 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,5570GT96.8% 521.6 / 3.9 157G49V (GGT→GTT) mutLmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (2/3);  new base T (79/73);  total (81/76)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 157 total aligned reads displayed.

GGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCG  >  NC_000913/4397460‑4397652
                                                                                                  |                                                                                               
ggtgaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTg                                                                                               <  2:282536/100‑1 (MQ=255)
 gtgaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTgg                                                                                              <  1:3036649/100‑1 (MQ=255)
 gtgaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTgg                                                                                              <  2:1196359/100‑1 (MQ=255)
    aggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGc                                                                                           >  1:2853323/1‑100 (MQ=255)
      gtggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGa                                                                                         >  1:2756177/1‑100 (MQ=255)
       tggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGaa                                                                                        <  2:715078/100‑1 (MQ=255)
        ggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGaaa                                                                                       >  1:1581381/1‑100 (MQ=255)
         gtCGAGCCACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAAc                                                                                      <  2:2312864/100‑1 (MQ=255)
           cGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACtt                                                                                    >  1:1779287/1‑100 (MQ=255)
           cGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACtt                                                                                    >  2:3085443/1‑100 (MQ=255)
             aGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTAt                                                                                  <  2:2190751/100‑1 (MQ=255)
                gACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCCGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCg                                                                               >  1:239393/1‑100 (MQ=255)
                      cGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCg                                                                         <  1:886656/100‑1 (MQ=255)
                        tCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTg                                                                       >  2:3291789/1‑100 (MQ=255)
                        tCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTg                                                                       <  2:2397602/100‑1 (MQ=255)
                         cGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGa                                                                      <  2:2875655/100‑1 (MQ=255)
                         cGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGa                                                                      <  1:2873221/100‑1 (MQ=255)
                          ggTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGAt                                                                     >  1:421332/1‑100 (MQ=255)
                           gTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATa                                                                    >  1:888108/1‑100 (MQ=255)
                           gTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGAAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATGCGTGATa                                                                    >  1:2300771/1‑100 (MQ=255)
                             aGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAAc                                                                  <  2:2157659/100‑1 (MQ=255)
                             aGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAAc                                                                  >  1:2340742/1‑100 (MQ=255)
                              gTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATAGCGAACGCGTGGGGGCGAAACTTATCCGCATTCGTGATAACg                                                                 >  2:2801024/1‑100 (MQ=255)
                                cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGc                                                               >  2:253583/1‑100 (MQ=255)
                                cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGc                                                               <  1:2536705/100‑1 (MQ=255)
                                cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGc                                                               >  2:1724415/1‑100 (MQ=255)
                                cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGc                                                               >  2:2905219/1‑100 (MQ=255)
                                cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGc                                                               >  1:1364402/1‑100 (MQ=255)
                                  aaGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTg                                                             >  1:1856526/1‑100 (MQ=255)
                                   aGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGc                                                            <  1:510430/100‑1 (MQ=255)
                                   aGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGc                                                            >  2:21107/1‑100 (MQ=255)
                                   aGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGAAGCGTATCGATATTGATATCGAACGCGTTGGGGAGAAACTTATACACATTCGTGATAAAGGCt                                                              >  2:2710287/1‑98 (MQ=255)
                                    gAACGAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAAAGCGTTGGGGCGAAACTTATCCGCATTCGCGATAACGGCTgtg                                                           <  1:2268325/100‑3 (MQ=255)
                                     aaCTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCgg                                                          <  1:1259288/100‑1 (MQ=255)
                                     aaCTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCgg                                                          <  1:1124255/100‑1 (MQ=255)
                                       cTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTa                                                        >  2:2541346/1‑100 (MQ=255)
                                        tAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTAt                                                       <  1:811381/100‑1 (MQ=255)
                                         aGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATc                                                      <  1:1395795/100‑1 (MQ=255)
                                          gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCa                                                     <  1:83116/100‑1 (MQ=255)
                                          gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCa                                                     >  2:2455218/1‑100 (MQ=255)
                                          gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCa                                                     >  1:3046980/1‑100 (MQ=255)
                                          gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCa                                                     <  1:1663086/100‑1 (MQ=255)
                                           tGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAATCTTATCCGCATTCGTGATAACGGCTGCGGTATCaa                                                    <  1:611564/100‑1 (MQ=255)
                                            ggAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaa                                                   >  2:2223435/1‑100 (MQ=255)
                                            ggAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaa                                                   >  1:1336610/1‑100 (MQ=255)
                                            ggAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaa                                                    <  2:2242309/99‑1 (MQ=255)
                                             gAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGAAGTATAaaaa                                                  <  2:2956288/100‑1 (MQ=255)
                                               aaaCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaaaaa                                                <  1:99981/100‑1 (MQ=255)
                                                aaCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAg                                               <  2:1122898/100‑1 (MQ=255)
                                                 aCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGa                                              >  2:1362213/1‑100 (MQ=255)
                                                 aCAGCCTCGATGCAGCTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGa                                              <  1:3297860/100‑1 (MQ=255)
                                                  cAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGAt                                             >  2:2338363/1‑100 (MQ=255)
                                                    gCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGa                                           >  1:2234249/1‑100 (MQ=255)
                                                      cTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgc                                         <  2:2319964/100‑1 (MQ=255)
                                                      cTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgc                                         >  1:1637152/1‑100 (MQ=255)
                                                        gatctcAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctg                                       <  1:450934/95‑1 (MQ=255)
                                                         gATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctgg                                      >  1:2596748/1‑100 (MQ=255)
                                                         gATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctgg                                      <  1:2807986/100‑1 (MQ=255)
                                                          aTGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggc                                     >  1:349121/1‑100 (MQ=255)
                                                          aTGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggc                                     <  1:3291565/100‑1 (MQ=255)
                                                             cAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctgagatc                                  >  2:450933/1‑95 (MQ=255)
                                                              cgccccaaCGCGTTCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctg                                 >  2:764705/8‑100 (MQ=25)
                                                               ggTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctgg                                >  2:1993873/1‑100 (MQ=255)
                                                               ggTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGGATCAAAAAAGATGAgctggcgctgg                                >  2:52754/1‑100 (MQ=255)
                                                                gTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCTGTATCAAAAAAGATGAgctggcgctg                                 <  2:1598502/98‑1 (MQ=255)
                                                                ccccaaCGCGTTCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggc                               >  1:2004971/6‑100 (MQ=25)
                                                                ccccaaCGCGTTCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggc                               >  1:1060397/6‑100 (MQ=37)
                                                                  gCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgc                             <  1:2711855/100‑1 (MQ=255)
                                                                   cGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgct                            <  1:6136/100‑1 (MQ=255)
                                                                    gACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctg                           >  1:1688293/1‑100 (MQ=255)
                                                                    gACGCATATCGATATTGATATCGAACGCGTTAGGGCGAAACTTATCCGCATTCGTGATAACGGAAAAAAAAAAAAAAAAGATGAgctggcgctggcg                              <  2:888106/97‑1 (MQ=255)
                                                                    aaCGCGTTCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctg                           >  1:217751/2‑100 (MQ=37)
                                                                    aaCGCGTTCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctg                           >  2:3006399/2‑100 (MQ=37)
                                                                      cgcgTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctggc                         >  2:1496676/1‑100 (MQ=255)
                                                                      cgcgTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctggc                         <  1:2709796/100‑1 (MQ=255)
                                                                      cgcgTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctggc                         >  1:879821/1‑100 (MQ=255)
                                                                        cgTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGggctggct                        >  2:2246607/1‑99 (MQ=255)
                                                                        cgTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTc                       >  2:1584116/1‑100 (MQ=255)
                                                                        cgTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTc                       >  1:3056126/1‑100 (MQ=255)
                                                                            tCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCa                   <  2:3138510/100‑1 (MQ=255)
                                                                              gATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCAt                  <  2:2170500/99‑1 (MQ=255)
                                                                              gATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATg                 <  1:1251064/100‑1 (MQ=255)
                                                                                tatTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcc               <  1:2739227/100‑1 (MQ=255)
                                                                                 atTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcca              <  1:324809/100‑1 (MQ=255)
                                                                                  ttGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccac             >  1:3347613/1‑100 (MQ=255)
                                                                                  ttGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccac             <  2:3276601/100‑1 (MQ=255)
                                                                                  ttGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccac             >  2:375073/1‑100 (MQ=255)
                                                                                  ttGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCATGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccac             >  2:2862828/1‑100 (MQ=255)
                                                                                   tGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacc            <  1:2754630/100‑1 (MQ=255)
                                                                                    gATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAGCGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCGCGTCATGccacc            <  2:2736541/99‑1 (MQ=255)
                                                                                    gATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacca           >  1:1824783/1‑100 (MQ=255)
                                                                                       atCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTa        >  2:848983/1‑100 (MQ=255)
                                                                                       atCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTa        >  1:2457451/1‑100 (MQ=255)
                                                                                       atCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTa        >  2:1516191/1‑100 (MQ=255)
                                                                                         cGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaaa      >  2:1049297/1‑100 (MQ=255)
                                                                                         cGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaaa      >  1:2794395/1‑100 (MQ=255)
                                                                                          gAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaaaa     >  2:661407/1‑100 (MQ=255)
                                                                                          gAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaaaa     <  2:1856519/100‑1 (MQ=255)
                                                                                           aaCGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCGGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAAt    <  1:1316868/100‑1 (MQ=255)
                                                                                              cgcgTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCg  <  1:3209767/100‑1 (MQ=255)
                                                                                                  |                                                                                               
GGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCG  >  NC_000913/4397460‑4397652

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