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.6% 484.9 / ‑1.4 138G49V (GGT→GTT) mutLmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base T (73/63);  total (75/63)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

Only 100 of 138 total aligned reads displayed.

GGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCG  >  NC_000913/4397460‑4397652
                                                                                                  |                                                                                               
ggtgaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTg                                                                                               <  1:1803673/100‑1 (MQ=255)
  tgaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGCTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTggg                                                                                             >  2:2294435/1‑100 (MQ=255)
   gaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTgggg                                                                                            <  1:3241739/100‑1 (MQ=255)
    aggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGc                                                                                           >  2:1928962/1‑100 (MQ=255)
        ggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGaaa                                                                                       >  2:1430996/1‑100 (MQ=255)
        ggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGaaa                                                                                       <  2:438320/100‑1 (MQ=255)
        ggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGaaa                                                                                       <  2:2166739/100‑1 (MQ=255)
         gtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAAc                                                                                      <  2:2932499/100‑1 (MQ=255)
          tCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACt                                                                                     >  2:66867/1‑100 (MQ=255)
          tCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACt                                                                                     <  2:1448389/100‑1 (MQ=255)
           cGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACtt                                                                                    >  1:441675/1‑100 (MQ=255)
           cGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACtt                                                                                    >  2:1102133/1‑100 (MQ=255)
              cgatcttGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATc                                                                                 <  2:2758099/94‑1 (MQ=255)
                gACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCg                                                                               <  1:2705781/100‑1 (MQ=255)
                 aCCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGc                                                                              >  1:217187/1‑100 (MQ=255)
                 aCCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGc                                                                              >  1:902032/1‑100 (MQ=255)
                  ccTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCa                                                                             >  1:2183287/1‑100 (MQ=255)
                    tGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCAtt                                                                           >  1:1121425/1‑100 (MQ=255)
                        cgatctgTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTg                                                                       <  1:3158709/94‑1 (MQ=255)
                            tAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATaa                                                                   >  2:1111065/1‑100 (MQ=255)
                             aGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAAc                                                                  <  1:955400/100‑1 (MQ=255)
                              gTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGagatcg                                                                 >  2:3158697/1‑95 (MQ=255)
                              gTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATagat                                                                 >  2:668501/1‑97 (MQ=255)
                              gTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACg                                                                 >  1:3403142/1‑100 (MQ=255)
                               tCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACgg                                                                <  1:3110946/100‑1 (MQ=255)
                                cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGc                                                               >  2:2778877/1‑100 (MQ=255)
                                  aaGAACGAGTGGAAAACAGCCTCGAGGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTg                                                             <  2:3103603/100‑1 (MQ=255)
                                   atctcTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGc                                                            <  2:229434/96‑1 (MQ=255)
                                   aGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGc                                                            <  1:1045402/100‑1 (MQ=255)
                                    gAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCg                                                           >  1:471642/1‑100 (MQ=255)
                                     aaCTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCgg                                                          <  1:364626/100‑1 (MQ=255)
                                      aCTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGt                                                         <  1:878662/100‑1 (MQ=255)
                                       cTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGcagat                                                        >  1:229435/1‑96 (MQ=255)
                                       cTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTa                                                        >  1:503353/1‑100 (MQ=255)
                                         aGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATc                                                      >  2:3182892/1‑100 (MQ=255)
                                          gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCa                                                     >  2:3265827/1‑100 (MQ=255)
                                          gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCa                                                     >  2:121034/1‑100 (MQ=255)
                                               aaaCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaaaaa                                                <  2:3403129/100‑1 (MQ=255)
                                               aaaCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaaaaa                                                <  1:2591693/100‑1 (MQ=255)
                                               aaaCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaaaaa                                                <  1:3265840/100‑1 (MQ=255)
                                               aaaCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaaaaa                                                <  1:1102135/100‑1 (MQ=255)
                                                aaCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAg                                               <  1:1928966/100‑1 (MQ=255)
                                                aaCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAg                                               <  1:1377665/100‑1 (MQ=255)
                                                aaCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGGATCAAAAAAg                                               >  2:1652200/1‑100 (MQ=255)
                                                  cAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGAt                                             >  2:2410862/1‑100 (MQ=255)
                                                   cgatctgATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATg                                            <  1:3100695/94‑1 (MQ=255)
                                                   aGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATg                                            <  1:66866/100‑1 (MQ=255)
                                                     ccTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAg                                          >  1:2166889/1‑100 (MQ=255)
                                                     ccTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAg                                          <  2:503352/100‑1 (MQ=255)
                                                        tgATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctg                                       <  1:1788682/99‑1 (MQ=255)
                                                         gATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctga                                      >  2:1788678/1‑99 (MQ=255)
                                                         gATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgct                                        <  2:660352/98‑1 (MQ=255)
                                                         gATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgatcg                                      >  2:3100683/1‑96 (MQ=255)
                                                            gCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgc                                   <  1:1516822/100‑1 (MQ=255)
                                                            gCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgc                                   >  2:1516819/1‑100 (MQ=255)
                                                             cAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgct                                  <  2:217186/100‑1 (MQ=255)
                                                             cAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgct                                  >  2:907184/1‑100 (MQ=255)
                                                             cAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgccagcgcc                                  >  1:758211/1‑93 (MQ=255)
                                                              aGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctg                                 >  1:454554/1‑100 (MQ=255)
                                                               ggTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctgg                                <  2:1879089/100‑1 (MQ=255)
                                                               ggTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctgg                                >  1:1845101/1‑100 (MQ=255)
                                                               ggTGCGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTAGTGATAACGGCTGCGGTATCAAAAAAGATgagctggagctg                                 >  2:3267750/1‑99 (MQ=255)
                                                                ccccaaCGCGT‑TCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggc                               >  1:3325499/6‑100 (MQ=25)
                                                                   cGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgct                            >  2:3274557/1‑100 (MQ=255)
                                                                   cGACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgct                            <  2:2504844/100‑1 (MQ=255)
                                                                    gACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctg                           <  2:892272/100‑1 (MQ=255)
                                                                    gACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctg                           >  2:3025801/1‑100 (MQ=255)
                                                                    gACGCGTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctg                           <  1:1032592/100‑1 (MQ=255)
                                                                      cgcgTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctggc                         <  1:267343/100‑1 (MQ=255)
                                                                      cgcgTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctggc                         >  1:1939078/1‑100 (MQ=255)
                                                                        gcgttCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCt                        >  1:1271677/5‑100 (MQ=37)
                                                                        gcgttCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCt                        >  1:556030/5‑100 (MQ=37)
                                                                         gTATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCg                      <  2:1939074/100‑1 (MQ=255)
                                                                          tATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGt                     >  2:1762661/1‑100 (MQ=255)
                                                                          tATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGt                     >  1:1470119/1‑100 (MQ=255)
                                                                          tATCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGt                     >  2:1520691/1‑100 (MQ=255)
                                                                           aTCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTc                    >  1:1310027/1‑100 (MQ=255)
                                                                           aTCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTc                    >  1:1561305/1‑100 (MQ=255)
                                                                            tCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCa                   >  2:1583576/1‑100 (MQ=255)
                                                                            tCGATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCa                   <  1:1111067/100‑1 (MQ=255)
                                                                             cGATATTGATATCGA‑ACGCGTTGGGGCGAAACCTATCCGTATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCAt                  <  2:2390442/100‑1 (MQ=255)
                                                                              gATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATg                 >  1:285161/1‑100 (MQ=255)
                                                                              gATATTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATg                 <  2:445887/100‑1 (MQ=255)
                                                                                tatTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcc               <  2:2183282/100‑1 (MQ=255)
                                                                                tatTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcc               <  1:2622119/100‑1 (MQ=255)
                                                                                 atTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcca              >  2:313346/1‑100 (MQ=255)
                                                                                 atTGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcc               >  1:1082651/1‑99 (MQ=255)
                                                                                  ttGATGTCGT‑AGGTGTTGGGGTGAAAGTTATCCGCATTCGTGATAGCGGCTGTGGGATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGgca              <  2:3016527/99‑1 (MQ=255)
                                                                                  ttGATATCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccac             >  2:366335/1‑100 (MQ=255)
                                                                                     atatCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAg          <  2:1121423/100‑1 (MQ=255)
                                                                                      tatCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGt         >  1:1103950/1‑100 (MQ=255)
                                                                                      tatCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGt         >  2:119853/1‑100 (MQ=255)
                                                                                      tatCGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGt         >  1:2008423/1‑100 (MQ=255)
                                                                                        tCGA‑ACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaa       >  1:2576803/1‑100 (MQ=255)
                                                                                         cGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaaa      >  1:2974598/1‑100 (MQ=255)
                                                                                         cGA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaaa      >  1:2882627/1‑100 (MQ=255)
                                                                                              cgcgTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCg  <  2:1197882/100‑1 (MQ=255)
                                                                                              cgcgTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCg  >  2:1675198/1‑100 (MQ=255)
                                                                                              cgcgTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCg  <  1:119852/100‑1 (MQ=255)
                                                                                              cgcgGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCGGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATc   >  2:1494356/1‑99 (MQ=255)
                                                                                                  |                                                                                               
GGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGA‑ACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCG  >  NC_000913/4397460‑4397652

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