Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,567,178 A→T 6.2% L103* (TTG→TAG) ‡ mutL ← DNA mismatch repair protein MutL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,567,1780AT 6.2% 192.4 / 3.0 64L103* (TTG→TAG) ‡mutLDNA mismatch repair protein MutL
Reads supporting (aligned to +/- strand):  ref base A (28/32);  new base T (2/2);  total (30/34)
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.75e-01

CACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTGCTGGTGGCGTGTCGAGCCAGGGCCAGCGGCAGGTCGTCGGCAGAAATGCCGCTGCCATCATCCCGCACCTTCAGCAGTTTGACGCC  >  NC_002947/5567033‑5567326
                                                                                                                                                    |                                                                                                                                                    
cACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAgg                                                                                                                                                   <  5:210063/149‑1 (MQ=255)
   ttCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCt                                                                                                                                                >  1:291927/1‑149 (MQ=255)
   ttCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCt                                                                                                                                                >  7:373070/1‑149 (MQ=255)
      cACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCCACCGAACTGATCGAGGCCAAGGCCTc                                                                                                                                               <  1:6255/148‑1 (MQ=255)
       aCCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGc                                                                                                                                             >  7:296404/1‑148 (MQ=255)
        ccGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCa                                                                                                                                           >  2:241273/1‑149 (MQ=255)
         cGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCTTGGCCTCgatca                                                                                                                                          <  3:155456/149‑5 (MQ=255)
         cGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCAc                                                                                                                                          <  3:182314/149‑1 (MQ=255)
           aGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCCCCGAACTGATCGAGGCCAAGGCCTCGCCACg                                                                                                                                         >  4:154323/1‑148 (MQ=255)
              tACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCTTGGCCTCgatc                                                                                                                                           <  3:293326/143‑4 (MQ=255)
              tACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCTTGGCCTCgatc                                                                                                                                           >  4:293326/1‑140 (MQ=255)
                ccgacctGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCTCTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAccc                                                                                                                                   <  1:30386/149‑1 (MQ=255)
                  gaccgGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTTCCAGGCCTCACTGGCACTGGCGCTACGCGAGGTCAGGGTACGCCGCGCCACCGCAAAGAACGAGGCCAAGGCCTCCACACGGACGCCCa                                                                                                                                 >  6:65079/1‑149 (MQ=255)
                  gaccgGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGcc                                                                                                                                   >  8:37662/1‑147 (MQ=255)
                    ccgGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAgg                                                                                                                               <  4:88455/149‑1 (MQ=255)
                     cgGGTGCAGCCGGCGGCGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGCTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGc                                                                                                                              >  6:291707/1‑149 (MQ=255)
                      gggTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCt                                                                                                                             <  7:37662/149‑1 (MQ=255)
                       ggTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTc                                                                                                                            <  4:344388/149‑1 (MQ=255)
                           cgccgcCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATa                                                                                                                        <  1:235506/149‑1 (MQ=255)
                              ccgccgGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATAcc                                                                                                                      >  3:330787/1‑149 (MQ=255)
                               cgccgGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGTCCTCGCCACGGAAGCCCAGGCTCAATAccc                                                                                                                     <  1:241273/149‑1 (MQ=255)
                                      tGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAg                                                                                                              >  3:346408/1‑149 (MQ=255)
                                            cgcgGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGACTCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCtt                                                                                                         <  1:105172/149‑1 (MQ=255)
                                            cgcgCCGTCAGGTCGCCGACTTCGGGATCTACTTCCCAGGCCTCACTGGCACTGGCGGTACGCGACGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTc                                                                                                        <  4:100682/149‑1 (MQ=255)
                                              cggcgTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCa                                                                                                      <  8:242811/149‑1 (MQ=255)
                                                   tCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCaa                                                                                                                          >  1:119449/1‑124 (MQ=255)
                                                   tCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCaa                                                                                                                          <  2:119449/124‑1 (MQ=255)
                                                   tCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTcg                                                                                                 >  6:81930/1‑149 (MQ=255)
                                                     aTGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTcgcg                                                                                               <  7:373544/149‑1 (MQ=255)
                                                     aTGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTcgc                                                                                                <  8:373070/148‑1 (MQ=255)
                                                        tCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGAt                                                                                            <  4:346408/149‑1 (MQ=255)
                                                          gcgGCCTTCGGTTTCTACTTGAGAGTCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATAACCATTCTAGGTCTTCCAGCTCGCGGATCt                                                                                          <  6:141530/149‑1 (MQ=255)
                                                            ggCCTTCGTTTTCTACTTGCCAGGCCTCACTGTCACTGGCGGTACGCTAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCAGCTCGCGGATCTTGc                                                                                       <  6:285891/149‑1 (MQ=255)
                                                                ttCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTGCtgg                                                                                    >  6:127155/1‑149 (MQ=255)
                                                                     tttCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTGCtgg                                                                                    <  2:396909/144‑1 (MQ=255)
                                                                     tttCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTGCtgg                                                                                    >  1:396909/1‑144 (MQ=255)
                                                                         tACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTGCTGGTGGCGTGTc                                                                           >  8:232266/1‑149 (MQ=255)
                                                                         tACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTGCTGGTGGCGTGTc                                                                           >  3:398709/1‑149 (MQ=255)
                                                                              gCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTGCTGGTGGCGTGTCGAg                                                                        <  8:296404/147‑1 (MQ=255)
                                                                              gCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTGCTGGTGGCGTGTCGAGc                                                                       >  4:299897/1‑148 (MQ=255)
                                                                              gCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTGCTGGTGGCGTGTCGAGc                                                                       <  5:81930/148‑1 (MQ=255)
                                                                               ccAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTcgc                                                                                                >  1:149442/1‑122 (MQ=255)
                                                                               ccAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTcgc                                                                                                <  2:149442/122‑1 (MQ=255)
                                                                                cAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTGCTGGTGGCGTGTCGAGCCAg                                                                    <  2:263741/149‑1 (MQ=255)
                                                                                        actggcactggcGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTGCTGGTGGCGTGTCGAGCCAGGGCCAg                                                              <  2:131405/147‑1 (MQ=255)
                                                                                          tggcactggcGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTGCtgg                                                                                    >  2:362316/1‑123 (MQ=255)
                                                                                          tggcactggcGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTGCtgg                                                                                    <  1:362316/123‑1 (MQ=255)
                                                                                          tggcactggcGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTGCTGGTGGCGTGTCGAGCCAGGGCCAGCgg                                                           >  2:121789/1‑148 (MQ=255)
                                                                                            gcactggcGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTGCTGGTGGCGTGTCGAGCCAGGGCCAGCGGCAg                                                        <  5:291707/149‑1 (MQ=255)
                                                                                              actggcGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTGCTGGTGGCGTGTCGAGCCAGGGCCAGCGGCAGgt                                                      >  6:25775/1‑149 (MQ=255)
                                                                                                tggcGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTGCTGGTGGCGTGTCGAGc                                                                       <  2:198366/130‑1 (MQ=255)
                                                                                                tggcGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTGCTGGTGGCGTGTCGAGc                                                                       >  1:198366/1‑130 (MQ=255)
                                                                                                               tcagggtcaggCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTGCTGGTGGCGTGTCGAGCCAGGGCCAGCGGCAGGTCGTCGGCAGAAATGCCg                                     <  5:25775/149‑1 (MQ=255)
                                                                                                               tcagggtcaggCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTGCTGGTGGCGTGTCGAGCCAGGGCCAGCGGCAGGTCGTCGGCAGAAATGCCg                                     >  2:342033/1‑149 (MQ=255)
                                                                                                                 agggtcaggCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTg                                                                                        >  4:11294/1‑96 (MQ=255)
                                                                                                                 agggtcaggCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTg                                                                                        <  3:11294/96‑1 (MQ=255)
                                                                                                                                 aCCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTGCTGGTGGCGTGTCGAGCCAGGGCCAGCGGCAGGTCGTCGGCAGAAATGCCGCTGCCATCATCCCGCAc                    <  2:291927/149‑1 (MQ=255)
                                                                                                                                          aTCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTGCTGGTGGCGTGTCGAGCCAGGGCCAGCGGCAGGTCGTCGGCAGAAATGCCGCTGCCATCATCCCGCACCTTcagcag           >  3:117098/1‑149 (MQ=255)
                                                                                                                                           ggcGAGGCCTTGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTGCTGGTGGCGTGTCGAGCCAGGGCCAGCGGCAGGTCGTCGGCAGAAATGCCGCTGCCATCATCCCGCACCTTcagcag           >  2:169828/3‑149 (MQ=255)
                                                                                                                                               aGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTGCTGGTGGCGTGTCGAGCCAGGGCCAGCGGCAGGTCGTCGGCAGAAATGCCGCTGCCATCATCCCGCACCTTc                <  7:130629/140‑1 (MQ=255)
                                                                                                                                               aGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTGCTGGTGGCGTGTCGAGCCAGGGCCAGCGGCAGGTCGTCGGCAGAAATGCCGCTGCCATCATCCCGCACCTTc                >  8:130629/1‑140 (MQ=255)
                                                                                                                                                   cAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTGCTGGTGGCGTGTCGAGCCAGGGCCAGCGGCAGGTCGTCGGCAGAAATGCCGCTGc                                 <  4:246051/119‑1 (MQ=255)
                                                                                                                                                   cAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTGCTGGTGGCGTGTCGAGCCAGGGCCAGCGGCAGGTCGTCGGCAGAAATGCCGCTGc                                 >  3:246051/1‑119 (MQ=255)
                                                                                                                                                   cAACGCCTCGCCACGGAAGCCCACGCTCAATACCCCTTCTAGGTCTACCAGCTCGCGGATCTAGCTGGTACCGTCTCGAGCCACGGCCAGCGGCAGGTCCTCCGCAGAAATCCCGCTGCCATCCTCCCGCCCCATCAGCAGTTTGAcgc   >  6:295671/1‑149 (MQ=255)
                                                                                                                                                    aaGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTGCTGGTGGCGTGTCGAGCCAGGGCCAGCGGCACGTCGTCCGCAGAAATGCCGCTGCCATCATCCCGCACCTTCAGCAGTTTGAcgcc  >  2:118573/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCTACTTGCCAGGCCTCACTGGCACTGGCGGTACGCGAGGTCAGGGTCAGGCGCGCCACCGAACTGATCGAGGCCAAGGCCTCGCCACGGAAGCCCAGGCTCAATACCCCTTCTAGGTCTTCCAGCTCGCGGATCTTGCTGGTGGCGTGTCGAGCCAGGGCCAGCGGCAGGTCGTCGGCAGAAATGCCGCTGCCATCATCCCGCACCTTCAGCAGTTTGACGCC  >  NC_002947/5567033‑5567326

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