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,5570GT98.4% 429.5 / ‑4.0 128G49V (GGT→GTT) mutLmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (75/51);  total (76/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 = 1.00e+00

Only 100 of 128 total aligned reads displayed.

GGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGC  >  NC_000913/4397460‑4397653
                                                                                                  |                                                                                                
ggtgaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTg                                                                                                >  2:623430/1‑100 (MQ=255)
   gaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTgggg                                                                                             <  1:2553499/100‑1 (MQ=255)
    aggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGc                                                                                            >  1:2876083/1‑100 (MQ=255)
     ggtggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGAGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGc                                                                                            >  1:1018228/1‑99 (MQ=255)
     ggtgggCGAGCGACCTGCGTCGGGAGTCAAAGAATTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑AGGCGGTGGGGcn                                                                                           <  1:1531539/100‑2 (MQ=255)
        ggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGaaa                                                                                        >  2:1878720/1‑100 (MQ=255)
        ggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGaaa                                                                                        >  2:1321381/1‑100 (MQ=255)
         gtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAAc                                                                                       <  1:1361836/100‑1 (MQ=255)
         gtCGAGCGACCCGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAAc                                                                                       <  1:1420948/100‑1 (MQ=255)
          tCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACt                                                                                      <  2:12216/100‑1 (MQ=255)
            gAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTa                                                                                    >  1:1913540/1‑100 (MQ=255)
             aGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTAt                                                                                   >  2:1470279/1‑100 (MQ=255)
             aGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTAt                                                                                   >  2:2216294/1‑100 (MQ=255)
             aGCGACCTGCGTCGGTAGTCAAAGAACTAGTGAAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTAt                                                                                   >  2:2752653/1‑100 (MQ=255)
               cGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATcc                                                                                 <  2:1284777/100‑1 (MQ=255)
                gACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCg                                                                                <  2:2100376/100‑1 (MQ=255)
                gACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCg                                                                                <  2:2445538/100‑1 (MQ=255)
                gACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCg                                                                                <  1:2300568/100‑1 (MQ=255)
                gACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCg                                                                                >  1:2630878/1‑100 (MQ=255)
                  ccTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCa                                                                              >  2:2736923/1‑100 (MQ=255)
                   cTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCAt                                                                             >  1:2160220/1‑100 (MQ=255)
                   cTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCAt                                                                             <  2:881327/100‑1 (MQ=255)
                     gCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGGTGGGGCGAAACTTATCCGCATTc                                                                           >  1:1356802/1‑100 (MQ=255)
                      cGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCg                                                                          <  2:1646312/100‑1 (MQ=255)
                       ttccgatctCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGt                                                                         <  1:781356/92‑1 (MQ=255)
                        tCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTg                                                                        >  1:1685608/1‑100 (MQ=255)
                        tCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTg                                                                        <  1:676242/100‑1 (MQ=255)
                        tCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTg                                                                        <  2:268736/100‑1 (MQ=255)
                         cGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGa                                                                       >  1:2389883/1‑100 (MQ=255)
                          ggTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGAt                                                                      <  1:491758/100‑1 (MQ=255)
                             aGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAAc                                                                   >  1:1622955/1‑100 (MQ=255)
                              gTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACg                                                                  <  2:827724/100‑1 (MQ=255)
                               tCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACgg                                                                 <  1:2338316/100‑1 (MQ=255)
                               tCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACgg                                                                 >  1:2073999/1‑100 (MQ=255)
                               tCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACg                                                                  >  1:2068484/1‑99 (MQ=255)
                                cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGc                                                                >  1:588972/1‑100 (MQ=255)
                                  aaGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTg                                                              >  1:716911/1‑100 (MQ=255)
                                   aGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGc                                                             >  1:2092312/1‑100 (MQ=255)
                                    gAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCg                                                            <  1:73513/100‑1 (MQ=255)
                                    gAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCg                                                            >  1:1572759/1‑100 (MQ=255)
                                     aaCTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCgg                                                           <  2:2074027/100‑1 (MQ=255)
                                       cTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTa                                                         >  2:1525987/1‑100 (MQ=255)
                                       cTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTa                                                         <  2:577876/100‑1 (MQ=255)
                                       cTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTa                                                         <  1:298533/100‑1 (MQ=255)
                                         aGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATc                                                       <  1:2752617/100‑1 (MQ=255)
                                         aGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATc                                                       >  1:2622104/1‑100 (MQ=255)
                                          gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCa                                                      <  1:966037/100‑1 (MQ=255)
                                          gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATCGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCa                                                      >  1:1591110/1‑100 (MQ=255)
                                            ggAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaa                                                    >  1:2786623/1‑100 (MQ=255)
                                             gAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaaa                                                   <  1:2735171/100‑1 (MQ=255)
                                              agaaCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaaaa                                                  <  2:1499241/98‑1 (MQ=255)
                                              aaaaCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaaaa                                                  <  2:2371598/100‑1 (MQ=255)
                                               aaaCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaaaaa                                                 <  2:1498179/100‑1 (MQ=255)
                                                 aCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGa                                               <  2:2845279/100‑1 (MQ=255)
                                                      cTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgc                                          >  2:562378/1‑100 (MQ=255)
                                                      cTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgc                                          >  1:1716381/1‑100 (MQ=255)
                                                        cGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctg                                        <  1:47134/100‑1 (MQ=255)
                                                        cGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgct                                         <  1:641793/99‑1 (MQ=255)
                                                          aTGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggc                                      <  2:1546450/100‑1 (MQ=255)
                                                           tGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcg                                     >  1:2483375/1‑100 (MQ=255)
                                                            gCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgc                                    >  2:2839984/1‑100 (MQ=255)
                                                             cAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgct                                   >  2:2522454/1‑100 (MQ=255)
                                                             cAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgct                                   >  1:2831318/1‑100 (MQ=255)
                                                               ggTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctgg                                 >  2:91532/1‑100 (MQ=255)
                                                               ggTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctgg                                 >  2:1223695/1‑100 (MQ=255)
                                                               ggTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctgg                                 <  1:2211132/100‑1 (MQ=255)
                                                                ccccaaCGCGT‑TCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggc                                >  2:285707/6‑100 (MQ=25)
                                                                ccccaaCGCGT‑TCGATATCAATATCGATACGCGTCGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggc                                >  2:1084618/6‑100 (MQ=25)
                                                                  ccaaCGCGT‑TCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgc                              >  2:389281/4‑100 (MQ=37)
                                                                   caaCGCGT‑TCGATATCAATATCGATACGAGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgct                             >  1:103453/3‑100 (MQ=255)
                                                                   cGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgct                             >  2:607502/1‑100 (MQ=255)
                                                                   cGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgct                             <  2:2731651/100‑1 (MQ=255)
                                                                    gACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctg                            <  2:1913567/100‑1 (MQ=255)
                                                                    aaCGCGT‑TCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctg                            >  2:1257306/2‑100 (MQ=37)
                                                                      cgcgTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctggc                          <  1:115718/100‑1 (MQ=255)
                                                                       gcgTATGGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCAAAAAAAAAAAAAAAGATGAGCTGGCGCTGGCGCTGGCt                         <  1:1099115/100‑1 (MQ=255)
                                                                       gcgTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCt                         <  2:63616/100‑1 (MQ=255)
                                                                        cgTATCGATATCGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGAATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTc                        >  1:1526637/1‑100 (MQ=255)
                                                                         gTATCGATATTGATATCGA‑ACGCGTTGGGGCGAGACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCg                       >  1:257097/1‑100 (MQ=255)
                                                                         gTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCg                       >  1:1181668/1‑100 (MQ=255)
                                                                         gTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCg                       <  2:2348047/100‑1 (MQ=255)
                                                                         gTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCg                       <  1:2403327/100‑1 (MQ=255)
                                                                         cgttCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCt                         >  1:1616419/4‑99 (MQ=38)
                                                                          tATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGt                      >  2:1905304/1‑100 (MQ=255)
                                                                          tATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGt                      <  1:2756643/100‑1 (MQ=255)
                                                                            tCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCa                    >  2:1006551/1‑100 (MQ=255)
                                                                             cGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCAt                   >  1:2012361/1‑100 (MQ=255)
                                                                                tatTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcc                >  2:1795458/1‑100 (MQ=255)
                                                                                  tcaatatCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccac              >  1:1233539/4‑100 (MQ=255)
                                                                                   tGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacc             <  1:1027585/100‑1 (MQ=255)
                                                                                   tGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacc             >  2:1318670/1‑100 (MQ=255)
                                                                                   tGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacc             <  2:411160/100‑1 (MQ=255)
                                                                                   tGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacc             <  1:1239578/100‑1 (MQ=255)
                                                                                   ccgatctgA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacc             <  2:1505992/93‑1 (MQ=255)
                                                                                    gATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacca            <  2:1146990/100‑1 (MQ=255)
                                                                                    gATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacca            >  1:157252/1‑100 (MQ=255)
                                                                                    gATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacca            >  1:503566/1‑100 (MQ=255)
                                                                                      tatCGA‑ACGCGTTGGGGCGAAACCTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCGGGCGCTGGCTCGTAATGacacca            >  1:2734776/1‑98 (MQ=255)
                                                                                        tCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaa        >  1:625680/1‑100 (MQ=255)
                                                                                               gcgTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGc  <  2:1018239/100‑1 (MQ=255)
                                                                                                  |                                                                                                
GGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGC  >  NC_000913/4397460‑4397653

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: ‑