Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,565,909 A→T 5.3% L526Q (CTG→CAG)  mutL ← DNA mismatch repair protein MutL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,565,9090AT 5.3% 243.4 / 2.4 80L526Q (CTG→CAG) mutLDNA mismatch repair protein MutL
Reads supporting (aligned to +/- strand):  ref base A (36/39);  new base T (2/2);  total (39/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.83e-01

AGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCAGCAGCGGCTGGCCGCTCAGGCCCTCGCTGGCCATGGCCACCTT  >  NC_002947/5565768‑5566048
                                                                                                                                              |                                                                                                                                           
agcagcTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTg                                                                                                                                      <  8:139607/149‑1 (MQ=255)
 gcagcTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCTGGGCCCCGAGACCCTGgcgattc                                                                                                                                     <  7:226165/149‑5 (MQ=255)
   agcTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCTGGGCCCCGAGACCCTGgcgattcgc                                                                                                                                   <  1:326654/149‑7 (MQ=255)
     cTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAAtt                                                                                                                                 <  8:154406/149‑1 (MQ=255)
     cTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAAtt                                                                                                                                 >  6:286824/1‑149 (MQ=255)
        aTTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCg                                                                                                                               >  4:113074/1‑148 (MQ=255)
        aTTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCg                                                                                                                               >  4:250519/1‑148 (MQ=255)
             ggTGCGCCTGAATGCGAGCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCGAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCGGGCGCTGCAATTCGAGGccc                                                                                                                         >  6:163815/1‑149 (MQ=255)
              gTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCAGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGcc                                                                                                                          >  5:217978/1‑147 (MQ=255)
                    cTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCt                                                                                                                  >  6:151295/1‑149 (MQ=255)
                    cTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCt                                                                                                                  >  2:23393/1‑149 (MQ=255)
                     tGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCt                                                                                                                  >  7:79253/1‑148 (MQ=255)
                       aaTGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCGTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGa                                                                                                                >  1:232001/1‑148 (MQ=255)
                        aTGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGaa                                                                                                               <  2:132596/148‑1 (MQ=255)
                        aTGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGaa                                                                                                               <  1:163608/148‑1 (MQ=255)
                           cGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCa                                                                                                           <  4:292438/149‑1 (MQ=255)
                             aTCACTGGTGCTGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACt                                                                                                         <  2:2684/149‑1 (MQ=255)
                              tCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACtg                                                                                                        <  2:285859/149‑1 (MQ=255)
                                aCTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACtgt                                                                                                       <  5:307083/148‑1 (MQ=255)
                                 cTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGc                                                                                                     <  5:151295/149‑1 (MQ=255)
                                    gTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATg                                                                                                  <  6:17285/149‑1 (MQ=255)
                                    gTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCAAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATg                                                                                                  <  1:68122/149‑1 (MQ=255)
                                    gTGCCGTACTCCATAAGGTCGGCAAGCAAGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTTGCGAATCGCCAGGGTCTCGGTGCCCAGGCTCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCAt                                                                                                   <  3:245553/148‑1 (MQ=255)
                                        cGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCATCAACCACTGTGCATGctct                                                                                              >  8:86688/1‑149 (MQ=255)
                                          tACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTc                                                                                            <  5:286824/149‑1 (MQ=255)
                                            cTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTGCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCg                                                                                           >  1:67186/1‑149 (MQ=255)
                                               cATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGca                                                                                        >  1:261865/1‑148 (MQ=255)
                                               cATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGca                                                                                        >  2:279512/1‑148 (MQ=255)
                                                   aGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTc                                                                                   <  2:261865/149‑1 (MQ=255)
                                                         gCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCtcc                                                                             >  7:283515/1‑147 (MQ=255)
                                                           aaGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCAc                                                                           >  7:31465/1‑149 (MQ=255)
                                                           aaGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGCCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCAc                                                                           >  2:52841/1‑149 (MQ=255)
                                                           aaGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGCCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCAc                                                                           >  2:57569/1‑149 (MQ=255)
                                                             gCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACgc                                                                         >  2:226772/1‑149 (MQ=255)
                                                              cACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACgct                                                                        <  6:217978/149‑1 (MQ=255)
                                                                cGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACgctgg                                                                      >  6:102528/1‑149 (MQ=255)
                                                                       cACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCTCTGCAATTCCTAGCCCAGGCGCTGAAACCACTGTGCATTCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAg                                                                 <  7:86688/147‑1 (MQ=255)
                                                                        aCCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGc                                                              >  8:293515/1‑149 (MQ=255)
                                                                           aGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGa                                                                  >  4:125076/1‑142 (MQ=255)
                                                                           aGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGa                                                                  <  3:125076/142‑1 (MQ=255)
                                                                           aGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAg                                                           <  3:250519/149‑1 (MQ=255)
                                                                           aGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAg                                                           >  5:17343/1‑149 (MQ=255)
                                                                                 ttGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACgctggct                                                                    <  6:116130/134‑1 (MQ=255)
                                                                                 ttGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACgctggct                                                                    >  5:116130/1‑134 (MQ=255)
                                                                                     cTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCgg                                                  >  2:148783/1‑148 (MQ=255)
                                                                                     cTTCAGCCTGCTTGAGCAAGGCCAGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGAGt                                                        <  2:26571/142‑1 (MQ=255)
                                                                                     cTTCAGCCTGCTTGAGCAAGGCCAGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGAGt                                                        >  1:26571/1‑142 (MQ=255)
                                                                                         cAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCAc                                               <  8:79253/148‑1 (MQ=255)
                                                                                          aGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACca                                             <  7:293515/149‑1 (MQ=255)
                                                                                            ccTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAgg                                                          >  6:54636/1‑134 (MQ=255)
                                                                                            ccTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAgg                                                          <  5:54636/134‑1 (MQ=255)
                                                                                            ccTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTcc                                                    >  4:31747/1‑140 (MQ=255)
                                                                                            ccTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTcc                                                    <  3:31747/140‑1 (MQ=255)
                                                                                               gCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCTCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACca                                             <  8:80503/144‑1 (MQ=255)
                                                                                               gCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACca                                             >  7:80503/1‑144 (MQ=255)
                                                                                                  tGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCAGCAGCGGc                                     >  1:113394/1‑149 (MQ=255)
                                                                                                    aGCAAGGCCGGGATCTGGCGAATCGCCAGGGCCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCAGCAGCGGCTg                                   <  1:57569/149‑1 (MQ=255)
                                                                                                    aGCAAGGCCGGGATCTGGCGAATCGCCAGGGCCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCAGCAGCGGCTg                                   <  1:52841/149‑1 (MQ=255)
                                                                                                     gCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCAGCAGCGGCTgg                                  <  3:216704/149‑1 (MQ=255)
                                                                                                        aGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCAGCAGCGGCTGGCCg                               >  6:218494/1‑149 (MQ=255)
                                                                                                          gCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCAGCAGCGGCTGGCCGCt                             >  3:179430/1‑149 (MQ=255)
                                                                                                             gggATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCAGCAGCGGCTGGCCGCTCAg                          <  5:328769/149‑1 (MQ=255)
                                                                                                              ggATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCAGCAGCGGCTGGCCGCTCAgg                         <  7:210940/149‑1 (MQ=255)
                                                                                                              ggAGCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCAGCAGCGGCTGGCCGCTCAgg                         <  5:163815/149‑1 (MQ=255)
                                                                                                               gATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCAGCAGCGGCTGGCCGCTCAGGc                        <  2:67186/149‑1 (MQ=255)
                                                                                                                aTCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCAGCAGCGGCTGGCCGCTCAGGcc                       >  4:48266/1‑149 (MQ=255)
                                                                                                                   tGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCAGCAGCGGCTGGCCGCTCAGGCCCTc                    >  4:120672/1‑149 (MQ=255)
                                                                                                                     gCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCAGCAGCGGCTGGCCGCTCAGGCCCt                     <  3:113074/146‑1 (MQ=255)
                                                                                                                     gCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCAGCAGCGGCTGGCCGCTCAGGCCCTCGc                  >  7:279959/1‑149 (MQ=255)
                                                                                                                     gCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCAGCAGCGGCCAGCCGCtgctggt                       <  1:148783/144‑7 (MQ=255)
                                                                                                                       gAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGTGTCTCCGGCACCAGCAGCGGCTGGCCGCTCAGGCCCTCGCtg                <  1:226772/149‑1 (MQ=255)
                                                                                                                          tCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCAGCAGCGGCTGGCCGCTCAGGCCCTCGCtggcc             <  5:102528/149‑1 (MQ=255)
                                                                                                                           cGCCTGGGCCCCGAGACCCTGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCAGCAGCGGCTGGCCGCTCAGGCCCTCGCtggcca            >  3:27494/1‑149 (MQ=255)
                                                                                                                           cGCCTGGGCCCCGAGACCCTGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGATTAACGCTGGATGAGGGCCAGGGTCTCCGGCCCCAGCAGCGGCTGGCCGCTCCGGCCAGCGCtggcca            >  4:219449/1‑149 (MQ=255)
                                                                                                                            gCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCAGCAGCGGCTGGCCGCTCAGGCCCTCGCtggccat           <  8:279959/149‑1 (MQ=255)
                                                                                                                            gCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCAGCAGCGGCTGGCCGCTCAGGCCCTCGCtggccat           <  5:218494/149‑1 (MQ=255)
                                                                                                                             ccAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCagcagc                                        >  2:85816/1‑119 (MQ=255)
                                                                                                                              cAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCAGCAGCGGCTGGCCGCTCAGGCCCTCGCtggccatgg         <  2:113394/149‑1 (MQ=255)
                                                                                                                                   tctcGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAAACACAGTGCATGCTCATCGGCACAGTCGGCTTCACCCAGGCTGAGGGCCAGGGTGTCCGGCACCCGCAGCGGCTGGCCGCGAGGGCGCTCGCtgggcgtggggggg    >  4:255614/1‑144 (MQ=255)
                                                                                                                                     tcGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCAGCAGCGGCTGGCCGCTCAGGCCCTCGCTGGCCATGGCCACCtt  >  8:169/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                           
AGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGGGCCCAGGCGCTGCAATTCGAAGCCCAGGCGCTGAAACCACTGTGCATGCTCTTCGGCACAGTCGGCTTCACGCTGGCTGAGGGCCAGGGTCTCCGGCACCAGCAGCGGCTGGCCGCTCAGGCCCTCGCTGGCCATGGCCACCTT  >  NC_002947/5565768‑5566048

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: