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

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

Only 100 of 134 total aligned reads displayed.

GGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCG  >  NC_000913/4397460‑4397652
                                                                                                   |                                                                                               
ggtgaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTg                                                                                               <  1:3010698/100‑1 (MQ=255)
 gtgaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTgg                                                                                              <  1:2941928/100‑1 (MQ=255)
   gaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTgggg                                                                                            <  2:2486650/100‑1 (MQ=255)
    aggtgGTCGAGCGACCTGCGTCGGTAGGCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGc                                                                                           <  1:1482277/100‑1 (MQ=255)
     ggtggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCg                                                                                          >  1:464780/1‑100 (MQ=255)
     ggtggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCg                                                                                          >  2:2964715/1‑100 (MQ=255)
      gtggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGa                                                                                         >  1:2488252/1‑100 (MQ=255)
       tggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGaa                                                                                        <  2:2540620/100‑1 (MQ=255)
         gtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAAc                                                                                      <  2:193112/100‑1 (MQ=255)
           cGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACtt                                                                                    >  1:127646/1‑100 (MQ=255)
            gAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTa                                                                                   <  1:2821716/100‑1 (MQ=255)
             aGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTAt                                                                                  >  1:644259/1‑100 (MQ=255)
             aGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTAt                                                                                  >  2:980555/1‑100 (MQ=255)
              gCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATc                                                                                 >  2:2382558/1‑100 (MQ=255)
              gCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATc                                                                                 >  2:1398358/1‑100 (MQ=255)
               cGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATcc                                                                                <  1:1133534/100‑1 (MQ=255)
                gACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCg                                                                               <  1:1332006/100‑1 (MQ=255)
                        tCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTg                                                                       <  2:205611/100‑1 (MQ=255)
                        tCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTg                                                                       <  2:2824373/100‑1 (MQ=255)
                          ggTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGAt                                                                     >  1:19275/1‑100 (MQ=255)
                          ggTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGAt                                                                     <  1:1631105/100‑1 (MQ=255)
                             aGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAAc                                                                  <  1:982577/100‑1 (MQ=255)
                             aGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAAc                                                                  <  1:2720085/100‑1 (MQ=255)
                              gTCAAAGGAGTAGTGAAAAACATCCTCGATGCAGGTGCGACGCATATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACg                                                                 <  1:556556/100‑1 (MQ=255)
                              gTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACg                                                                 >  1:1122267/1‑100 (MQ=255)
                               tCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACgg                                                                >  2:1875703/1‑100 (MQ=255)
                                cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGc                                                               <  2:2794784/100‑1 (MQ=255)
                                cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGc                                                               >  2:2113071/1‑100 (MQ=255)
                                cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGc                                                               >  1:615015/1‑100 (MQ=255)
                                 aaaGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCt                                                              >  2:674321/1‑100 (MQ=255)
                                    gAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCg                                                           <  1:1007478/100‑1 (MQ=255)
                                    gAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCg                                                           <  2:1227938/100‑1 (MQ=255)
                                        tAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTAt                                                       <  2:1414315/100‑1 (MQ=255)
                                        tAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTAt                                                       <  1:2261471/100‑1 (MQ=255)
                                         aGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATc                                                      >  2:2319126/1‑100 (MQ=255)
                                          gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCa                                                     <  1:275247/100‑1 (MQ=255)
                                          gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCa                                                     >  1:277683/1‑100 (MQ=255)
                                          gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCa                                                     >  1:2449718/1‑100 (MQ=255)
                                          gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCa                                                     >  1:665108/1‑100 (MQ=255)
                                          gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCa                                                     >  1:556566/1‑100 (MQ=255)
                                           tGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaa                                                    >  1:151739/1‑100 (MQ=255)
                                            ggAAAACGGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaa                                                   >  2:164749/1‑100 (MQ=255)
                                             gatgagaGCGTCGATGCAGGTGCGACGGGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaa                                                   <  1:2550595/93‑1 (MQ=255)
                                                 aCAGCCTCCATGCAGGTGCGACGCGTATCGATATTGATCTCG‑A‑ACGCGTTGGGGAGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGa                                              <  2:562535/100‑1 (MQ=255)
                                                  cAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGAt                                             <  2:631433/100‑1 (MQ=255)
                                                  cAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGAt                                             <  2:1701298/100‑1 (MQ=255)
                                                  cAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGAt                                             <  1:297850/100‑1 (MQ=255)
                                                   aGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATG                                            <  1:1916173‑M1/100‑1 (MQ=255)
                                                   aGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATG                                            <  2:844029‑M1/100‑1 (MQ=255)
                                                   aGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATG                                            <  2:2544455‑M1/100‑1 (MQ=255)
                                                    gCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGA                                           >  1:1246151‑M1/1‑100 (MQ=255)
                                                    gCCTCGATGCAGGGGCGACGCGTATGGATATTGATATCG‑A‑GCGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATG                                            <  1:268789‑M1/99‑1 (MQ=255)
                                                     ccTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAG                                          >  2:283100‑M1/1‑100 (MQ=255)
                                                     ccTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAG                                          >  2:1601858‑M1/1‑100 (MQ=255)
                                                      cTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGC                                         >  2:192659‑M1/1‑100 (MQ=255)
                                                       tCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGC                                         <  1:2257696‑M1/99‑1 (MQ=255)
                                                          aTGCAGGTGCGACCCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggc                                     <  1:2880666/100‑1 (MQ=255)
                                                           tGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcg                                    >  1:331503/1‑100 (MQ=255)
                                                               ggTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCACAAAAGATGAgctggcgctgg                                >  2:1782001/1‑100 (MQ=255)
                                                               ggTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctgg                                >  1:716640/1‑100 (MQ=255)
                                                                gTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggc                               <  2:1387124/100‑1 (MQ=255)
                                                                gTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggc                               <  2:2136253/100‑1 (MQ=255)
                                                                ccccaaCGCGT‑TCGATATCAATATCG‑ATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggc                               >  2:625156/6‑100 (MQ=25)
                                                                 tGCGACGCGTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcg                              >  2:797571/1‑100 (MQ=255)
                                                                      cgcgTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGgc                         <  1:2227064‑M2/100‑3 (MQ=255)
                                                                      cgcgT‑TCGATATCAATATCG‑ATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGgc                         <  1:812250‑M2/100‑3 (MQ=255)
                                                                        cgTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGgctc                       >  1:2568179‑M2/1‑96 (MQ=255)
                                                                         gTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCg                      <  2:2042026/100‑1 (MQ=255)
                                                                         gTATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCg                      <  1:283098/100‑1 (MQ=255)
                                                                          tATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGt                     <  2:2255636/100‑1 (MQ=255)
                                                                          tATCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGt                     <  2:918775/100‑1 (MQ=255)
                                                                            tCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCcagcgccagcgccagctcatct                   <  2:1809721‑M1/100‑23 (MQ=255)
                                                                            tCGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCa                   <  2:1727999/100‑1 (MQ=255)
                                                                             cGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTcct                  <  2:2060215/100‑3 (MQ=255)
                                                                             cGATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCAt                  <  2:2923806/100‑1 (MQ=255)
                                                                              gATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATg                 >  2:1094055/1‑100 (MQ=255)
                                                                              gATATTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCAGAAGAATAAAAAAAGATGAGCcagcgccagcgccagctcatct                   <  1:1817335‑M1/98‑23 (MQ=255)
                                                                              gATATTGATATCG‑A‑ACGCGCTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCcagcgccagcgccagctcatcttt                 <  1:2674195‑M1/100‑25 (MQ=255)
                                                                                tatcaatatCGAA‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGc                >  1:2295506/6‑100 (MQ=255)
                                                                                tatTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCcagcgccagcgccagctcatcttttt               <  1:1398352‑M1/100‑27 (MQ=255)
                                                                                tatTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcc               <  2:728325/100‑1 (MQ=255)
                                                                                tatTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcc               <  1:2178667/100‑1 (MQ=255)
                                                                                tatTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcc               <  2:1246157/100‑1 (MQ=255)
                                                                                 atTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcca              >  1:969168/1‑100 (MQ=255)
                                                                                 atTGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTCATAACGGCTGCGGTATCAAAAAAGATGAGCcagcgccagcgccagctcatctttttt              <  1:2152759‑M1/100‑28 (MQ=255)
                                                                                 atTGATATCG‑A‑ACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCcagcgccagcgccagctcatctttttt              <  2:2987154‑M1/100‑28 (MQ=255)
                                                                                  ttGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCcagcgccagcgccagctcatcttttttg             <  2:3030665‑M1/100‑29 (MQ=255)
                                                                                   tGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacc            <  2:1644013/100‑1 (MQ=255)
                                                                                   tGATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacc            >  2:375765/1‑100 (MQ=255)
                                                                                    gATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacca           >  1:2948918/1‑100 (MQ=255)
                                                                                    gATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacca           >  2:3016751/1‑100 (MQ=255)
                                                                                    gATATCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGAAGAGCTGGCGCTGGCGCTGGCTCGTCATGccacca           >  1:1613823/1‑100 (MQ=255)
                                                                                     atatCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAg          <  2:1725890/100‑1 (MQ=255)
                                                                                     atatCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAg          >  2:217693/1‑100 (MQ=255)
                                                                                      tatCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGt         >  2:1776388/1‑100 (MQ=255)
                                                                                       atCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTa        >  1:121834/1‑100 (MQ=255)
                                                                                        tCG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaa       >  1:2418335/1‑100 (MQ=255)
                                                                                         cG‑A‑ACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTa        <  1:650479/98‑1 (MQ=255)
                                                                                               cgcgTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATc   >  2:49889/1‑99 (MQ=255)
                                                                                               cgcgTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCg  <  2:331505/100‑1 (MQ=255)
                                                                                                   |                                                                                               
GGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCG‑A‑ACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCG  >  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: ‑