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,5570GT99.1% 367.8 / NA 109G49V (GGT→GTT) mutLmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (60/48);  total (61/48)
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 109 total aligned reads displayed.

GGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCG  >  NC_000913/4397460‑4397652
                                                                                                  |                                                                                               
 gtgaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTgt                                                                                              <  2:1242691/100‑2 (MQ=255)
  tgaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTggg                                                                                             <  2:2552037/100‑1 (MQ=255)
   gaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTgggg                                                                                            >  1:1209531/1‑100 (MQ=255)
   gaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAAATAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTgggg                                                                                            <  1:2293948/100‑1 (MQ=255)
     ggtggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCg                                                                                          >  2:1972377/1‑100 (MQ=255)
     ggtggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCg                                                                                          >  2:272027/1‑100 (MQ=255)
        ggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGaaa                                                                                       >  2:2719978/1‑100 (MQ=255)
        ggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGaaa                                                                                       >  2:2531202/1‑100 (MQ=255)
          tCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACt                                                                                     <  2:538525/100‑1 (MQ=255)
          tCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACt                                                                                     <  2:1209517/100‑1 (MQ=255)
           cGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACtt                                                                                    >  2:2790594/1‑100 (MQ=255)
           cGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACtt                                                                                    >  2:1094597/1‑100 (MQ=255)
                gACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCg                                                                               <  2:1599106/100‑1 (MQ=255)
                 aCCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGc                                                                              >  2:2182901/1‑100 (MQ=255)
                  ccTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCa                                                                             <  2:2995391/100‑1 (MQ=255)
                   cTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCAt                                                                            >  2:1236051/1‑100 (MQ=255)
                     gCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTc                                                                          >  2:1247561/1‑100 (MQ=255)
                      cGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCg                                                                         <  1:2065067/100‑1 (MQ=255)
                      cGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCg                                                                         <  2:3184445/100‑1 (MQ=255)
                       gTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGt                                                                        <  2:1651779/100‑1 (MQ=255)
                        tCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTg                                                                       <  2:1766467/100‑1 (MQ=255)
                        tCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTg                                                                       <  1:244541/100‑1 (MQ=255)
                        tCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTg                                                                       <  2:1830193/100‑1 (MQ=255)
                         cGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGa                                                                      <  1:3124696/100‑1 (MQ=255)
                             aGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAAc                                                                  >  2:2083762/1‑100 (MQ=255)
                              gTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACg                                                                 >  2:495887/1‑100 (MQ=255)
                              gTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACg                                                                 >  1:295745/1‑100 (MQ=255)
                              gTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACg                                                                 >  2:380106/1‑100 (MQ=255)
                                cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGc                                                               >  2:1301493/1‑100 (MQ=255)
                                     aaCTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCgg                                                          <  2:2369008/100‑1 (MQ=255)
                                     aaCTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCgg                                                          <  2:2583118/100‑1 (MQ=255)
                                      aCTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGt                                                         >  1:1943648/1‑100 (MQ=255)
                                       cTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTAATAACGGCTGCGGTa                                                        >  2:2236183/1‑100 (MQ=255)
                                        tAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTAt                                                       >  2:1167398/1‑100 (MQ=255)
                                         aGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATc                                                      <  2:1296631/100‑1 (MQ=255)
                                         aGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATc                                                      <  2:1201893/100‑1 (MQ=255)
                                         aGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTAt                                                       >  2:2021757/1‑99 (MQ=255)
                                          ntGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCa                                                     >  2:2840079/2‑100 (MQ=255)
                                          gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCa                                                     <  2:2372654/100‑1 (MQ=255)
                                          gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCa                                                     <  2:1043936/100‑1 (MQ=255)
                                           tGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaa                                                    >  1:2937946/1‑100 (MQ=255)
                                            ggAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaa                                                   <  2:2744718/100‑1 (MQ=255)
                                              aaaaCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGCAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaaaa                                                 >  1:2854878/1‑100 (MQ=255)
                                               aaaCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaaaaa                                                <  1:1236065/100‑1 (MQ=255)
                                                aaCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAg                                               >  1:1142973/1‑100 (MQ=255)
                                                aaCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAg                                               <  1:2922088/100‑1 (MQ=255)
                                                 aCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGa                                              <  2:1370780/100‑1 (MQ=255)
                                                 aCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGa                                              <  1:1972400/100‑1 (MQ=255)
                                                  cAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGAt                                             >  2:2866723/1‑100 (MQ=255)
                                                  cAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGAt                                             <  1:2996507/100‑1 (MQ=255)
                                                   aGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATg                                            <  1:1975310/100‑1 (MQ=255)
                                                   aGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATg                                            <  1:2034725/100‑1 (MQ=255)
                                                   aGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATg                                            <  2:3036213/100‑1 (MQ=255)
                                                    gCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATACCGGTTGCGGTATCAAAAAAGATGa                                           >  1:445476/1‑100 (MQ=255)
                                                     ccTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAg                                          >  1:75665/1‑100 (MQ=255)
                                                     ccTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAg                                          >  2:248010/1‑100 (MQ=255)
                                                      cTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgc                                         >  1:114527/1‑100 (MQ=255)
                                                       tCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgct                                        >  1:879047/1‑100 (MQ=255)
                                                        ngATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctg                                       >  2:3145475/2‑100 (MQ=255)
                                                        cGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctg                                       <  1:617349/100‑1 (MQ=255)
                                                          aTGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggc                                     >  1:2748493/1‑100 (MQ=255)
                                                           tGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcg                                    >  2:2760835/1‑100 (MQ=255)
                                                             cAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgct                                  >  1:2518652/1‑100 (MQ=255)
                                                             cAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgct                                  >  2:2650313/1‑100 (MQ=255)
                                                             cAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTACCAAAAAAGATGAgctggcgct                                  >  1:2443554/1‑100 (MQ=255)
                                                               ggTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctgg                                >  2:564018/1‑100 (MQ=255)
                                                               ggTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctgg                                <  1:1179462/100‑1 (MQ=255)
                                                               gccccaaCGCGTTCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctgg                                >  1:397533/7‑100 (MQ=25)
                                                                 tGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcg                              >  2:1151125/1‑100 (MQ=255)
                                                                 cccaaCGCGTTCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcg                              >  1:1964784/5‑100 (MQ=25)
                                                                   cGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgct                            >  1:361784/1‑100 (MQ=255)
                                                                    gaagcgTATGGATATTGATATCGAACGCGTTGAGGCGAAACTTATCCGCGTTCGTGATTACGGAAAAAAAAAAAAAAAAGATGAgctggcgctggcgctg                           <  2:437797/97‑1 (MQ=255)
                                                                    gACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctg                           <  1:2178887/100‑1 (MQ=255)
                                                                      cgcgTATTGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctggc                         >  1:2847038/1‑100 (MQ=255)
                                                                       gcgTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCt                        <  1:1539080/100‑1 (MQ=255)
                                                                       gcgTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCAGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCt                        <  2:1142959/100‑1 (MQ=255)
                                                                       cgcgttCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctggc                         >  1:3128059/6‑100 (MQ=37)
                                                                        gcgttCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCt                        >  2:832636/5‑100 (MQ=37)
                                                                         gTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCg                      >  1:1101682/1‑100 (MQ=255)
                                                                         gTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCg                      <  1:479071/100‑1 (MQ=255)
                                                                         gTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCg                      <  1:2720012/100‑1 (MQ=255)
                                                                          tATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGt                     <  1:2083788/100‑1 (MQ=255)
                                                                           aTCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTc                    >  1:1093451/1‑100 (MQ=255)
                                                                                tatTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcc               >  2:2575933/1‑100 (MQ=255)
                                                                                tatTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcc               <  2:445469/100‑1 (MQ=255)
                                                                                 atTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcca              <  2:298993/100‑1 (MQ=255)
                                                                                  tcaatatCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcca              >  1:1120751/4‑100 (MQ=255)
                                                                                   tGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacn            <  1:3145512/100‑2 (MQ=255)
                                                                                   tGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGGCATGccac             >  1:1933410/1‑99 (MQ=255)
                                                                                   caatatCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacc            >  2:3176048/3‑100 (MQ=255)
                                                                                    gATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacca           >  2:1199/1‑100 (MQ=255)
                                                                                     atatCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAg          <  1:1181706/100‑1 (MQ=255)
                                                                                        tCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaa       <  2:1369797/100‑1 (MQ=255)
                                                                                         cGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaaa      >  2:1841902/1‑100 (MQ=255)
                                                                                         cGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaaa      >  1:141589/1‑100 (MQ=255)
                                                                                          gAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaaaa     >  1:1638388/1‑100 (MQ=255)
                                                                                          gAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaaaa     <  1:2155796/100‑1 (MQ=255)
                                                                                              cgcgTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATcn  <  2:816911/100‑2 (MQ=255)
                                                                                              cgcgTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCg  >  1:998635/1‑100 (MQ=255)
                                                                                              cgcgTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCg  <  2:2748459/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: