Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,566,961 C→G 6.1% L175F (TTG→TTC)  mutL ← DNA mismatch repair protein MutL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,566,9610CG 6.1% 205.9 / 3.6 66L175F (TTG→TTC) mutLDNA mismatch repair protein MutL
Reads supporting (aligned to +/- strand):  ref base C (30/32);  new base G (2/2);  total (32/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 = 6.82e-01

AAGCCCGGGCCGCAGATGGCACCCACCCGCCGTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTT  >  NC_002947/5566815‑5567102
                                                                                                                                                   |                                                                                                                                              
aaGCGCGGGCCGCAGATGGCACCCACCCGCCGTGCACGGGCGATTTAGTCGTGGGCCTGGTGCAGGCTTAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCTGCGCCAGCCGCCGGATCACTTCCTGCat                                                                                                                                              <  2:309876/149‑2 (MQ=255)
 aGCCCGGGCCGCAGATGGCACCCACCCGCCGTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAg                                                                                                                                             >  1:40282/1‑149 (MQ=255)
  gCCCGGGCCGCAGATGGCACCCACCCGCCGTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGt                                                                                                                                            >  4:385604/1‑149 (MQ=255)
      gggCCGCAGATGGCACCCACCCGCCGTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATc                                                                                                                                        <  6:211145/149‑1 (MQ=255)
      gggCCGCAGATGGCACCCACCCGCCGTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATc                                                                                                                                        <  2:40282/149‑1 (MQ=255)
      gggCCGCAGATGGCACCCACCCGCCGTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATc                                                                                                                                        >  2:35920/1‑149 (MQ=255)
        gCCGCAGATGGCACCCACCCGCCGTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCg                                                                                                                                       >  3:385009/1‑148 (MQ=255)
         ccGCAGATGGCACCCACCCGCCGTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGaa                                                                                                                                     <  6:62053/149‑1 (MQ=255)
          cGCAGATGGCACCCACCCGCCGTGCACGGGCGATTTCGTCGTGGGTCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAt                                                                                                                                    <  3:403225/149‑1 (MQ=255)
           gCAGATGGCACCCACCCGCCGTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCCGATCACTTGCAGGAAGTGATccggcg                                                                                                                                   <  3:261495/149‑6 (MQ=255)
            cAGATGGCACCCAGCCGCCGTGCACGGTCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑Tc                                                                                                                                 <  7:199492/149‑1 (MQ=255)
            cAGATGGCACCCACCCGCCGTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑Tc                                                                                                                                 >  3:270110/1‑149 (MQ=255)
                tGGCACCCACCCGCCGTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGtt                                                                                                                             <  5:31612/149‑1 (MQ=255)
                       cacccGCCGTGCACGGGCGATTTCGTCGTGGGCCTAGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATAACTTCCTGCAAGTGATCGAAT‑TAGGttttatagg                                                                                                                      >  2:306707/1‑144 (MQ=255)
                        acccGCCGTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGc                                                                                                                     >  1:360574/1‑149 (MQ=255)
                             ccgTGCACGGGCGATTTCGTCGTGGGCCTCGTGCGGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGa                                                                                                                <  5:86368/149‑1 (MQ=255)
                             ccgTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGa                                                                                                                <  2:360574/149‑1 (MQ=255)
                             ccgTGCACGGGCGATTTCGTCGTGGGCCTCGTACAGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGa                                                                                                                <  2:155069/149‑1 (MQ=255)
                              cgTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAg                                                                                                               <  8:31099/149‑1 (MQ=255)
                               gTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGACAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTATTGGCTTTGAg                                                                                                               >  7:152583/1‑148 (MQ=255)
                                   aCGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACt                                                                                                          <  1:306707/149‑1 (MQ=255)
                                    cGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACtt                                                                                                         <  7:195499/149‑1 (MQ=255)
                                     gggCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTg                                                                                                        <  3:344652/149‑1 (MQ=255)
                                         gATTTCGTCGTGGGCCTCGAGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTACCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGcc                                                                                                    >  5:46489/1‑149 (MQ=255)
                                               gtcgtGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGtttt                                                                                                                           >  6:355894/1‑120 (MQ=255)
                                               gtcgtGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGtttt                                                                                                                           <  5:355894/120‑1 (MQ=255)
                                                  gtGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCggg                                                                                            <  2:223375/148‑1 (MQ=255)
                                                       ccTCGTGCAGGCTGAGGTTGCTCTTGCCATTGTGGCGCA‑GGTGGAAGGCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTg                                                                                      <  6:46489/149‑1 (MQ=255)
                                                           gTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTgaaga                                                                                  <  3:385604/149‑1 (MQ=255)
                                                           gTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTgaa                                                                                    <  3:202399/147‑1 (MQ=255)
                                                             gCAGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTgaagaa                                                                                 >  7:111457/1‑148 (MQ=255)
                                                                  cTGAGGATGCTCTTGCCATTGTGGCGCAGGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAAt                                                                            >  1:24937/1‑149 (MQ=255)
                                                                     aGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATcgc                                                                         <  3:303627/148‑1 (MQ=255)
                                                                                    ttGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCTGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAgg                                                          >  1:287721/1‑148 (MQ=255)
                                                                                           gcA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATc                                                                           <  1:77897/124‑1 (MQ=255)
                                                                                           gcA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATc                                                                           >  2:77897/1‑124 (MQ=255)
                                                                                           gcA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTaccgacc                                                  >  7:162635/1‑149 (MQ=255)
                                                                                                gTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCCCTTACACCGAGGTACCGACCGGGt                                              >  7:242811/1‑149 (MQ=255)
                                                                                                gTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGt                                              <  5:370555/149‑1 (MQ=255)
                                                                                                   gAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTgcg                                           <  8:111457/149‑1 (MQ=255)
                                                                                                     aGCCGACATCGAAACGGGCCAGAGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTTGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGcgcc                                         <  8:152583/149‑1 (MQ=255)
                                                                                                        cGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCgccgcc                                      >  2:30386/1‑149 (MQ=255)
                                                                                                        cGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCgccgcc                                      <  3:912/149‑1 (MQ=255)
                                                                                                         gACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGccgccg                                     >  5:344708/1‑149 (MQ=255)
                                                                                                           cATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGCCCGGGTGCGCCGCCGGc                                   >  3:344388/1‑149 (MQ=255)
                                                                                                            aTCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCt                                  <  6:344708/149‑1 (MQ=255)
                                                                                                            aTCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCt                                  >  5:38892/1‑149 (MQ=255)
                                                                                                             tCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTg                                 >  5:130124/1‑149 (MQ=255)
                                                                                                             tCGAAACGCGCCACCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATAGCGGTTTTCTCGGCTTTGAGGAACATTCGCCTGGCCGGGGTATTCAAGAACAACTCGCCCACTTCCACCGAGGTACCGACCGGGTGAGCCGCCGGCt                                  >  1:100889/1‑149 (MQ=255)
                                                                                                                aaaCGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCGCCGAGGTACCGACCGGGTGCGCCGCCGGCTGCAc                              >  5:174715/1‑149 (MQ=255)
                                                                                                                aaaCGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCAc                              >  6:277470/1‑149 (MQ=255)
                                                                                                                aaaCGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCAc                              >  1:68195/1‑149 (MQ=255)
                                                                                                                    gggCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTg                                             >  6:106279/1‑130 (MQ=255)
                                                                                                                    gggCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTg                                             <  5:106279/130‑1 (MQ=255)
                                                                                                                       ccagcgccagcCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCgcggc                       <  3:30883/149‑1 (MQ=255)
                                                                                                                          gcgccagcCGCGGGATCACTTCCTGCTAGTGATGGAAT‑TCGGTTTTCTCGGCTTTGAGGTTGTTGCGCCTGGCCGGGGTATTGAAGAACATTTCGCGCTCTTCCACCGAGGTACGGACGGGGTGCGCCGCCGGCTGCACCCGCGGCGTc                    <  2:100889/149‑1 (MQ=255)
                                                                                                                           cgccagcCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCCACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCa                   >  1:263741/1‑149 (MQ=255)
                                                                                                                            gccagcCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCAt                  >  6:61821/1‑149 (MQ=255)
                                                                                                                               agcCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTc               >  2:6255/1‑149 (MQ=255)
                                                                                                                                  cgccgGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTcgcg            <  6:38892/149‑1 (MQ=255)
                                                                                                                                      ggATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCCGCGACATGAcgccgcgg        <  1:35920/149‑3 (MQ=255)
                                                                                                                                      cgATCACTTGCAGGAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTcgcg            <  1:139021/144‑1 (MQ=255)
                                                                                                                                      cgATCACTTGCAGGAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTcgcg            >  2:139021/2‑145 (MQ=255)
                                                                                                                                        aTCACTTGCAGGAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTc      >  8:39922/1‑149 (MQ=255)
                                                                                                                                           aCTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGt   <  5:277470/149‑1 (MQ=255)
                                                                                                                                             ttCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGtt  >  2:235506/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
AAGCCCGGGCCGCAGATGGCACCCACCCGCCGTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCA‑GGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAT‑TCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTT  >  NC_002947/5566815‑5567102

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