Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,566,959 T→A 7.2% Q176L (CAG→CTG)  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,9590TA 7.2% 238.9 / 7.3 83Q176L (CAG→CTG) mutLDNA mismatch repair protein MutL
Reads supporting (aligned to +/- strand):  ref base T (41/36);  new base A (3/3);  total (44/39)
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.77e-02

CATGAAGCCCGGGCCGCAGATGGCACCCACCCGCCGTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCT  >  NC_002947/5566811‑5567105
                                                                                                                                                    |                                                                                                                                                    
cATGAAGCCCGGGCCGCAGATGGCACCCACCCTCCGTGCAAGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCt                                                                                                                                                      >  5:7075/1‑149 (MQ=255)
  tGAAGCCCGGGCCGCAGATGGCACCCACCCGCCGTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCTTTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGc                                                                                                                                                    <  3:225524/149‑1 (MQ=255)
    aaGCCCGGGCCGCAGATGGCACCCACCCGCCGTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGc                                                                                                                                                    >  2:328356/1‑147 (MQ=255)
      gCCCGGGCCGCAGATGGCACCCACCCGCCGTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTGCAGGAAGt                                                                                                                                                <  6:108527/149‑1 (MQ=255)
          gggCCGCAGATGGCACCCACCCGCCGTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATc                                                                                                                                           <  5:45831/149‑1 (MQ=255)
            gCCGCAGATGGCACCCACCCGCCGTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTGCAGGAAGTGATccg                                                                                                                                         <  3:84914/149‑3 (MQ=255)
              cGCAGATGGCACCCACCCGCCGTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTGCAGGAAGTGATccg                                                                                                                                         <  1:110916/147‑3 (MQ=255)
              cGCAGATGGCACCCACCCGCCGTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTGCAGGAAGTGATccg                                                                                                                                         >  2:110916/1‑145 (MQ=255)
              cGCAGATGGCACCCACCCGCCGTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAt                                                                                                                                       >  7:22694/1‑149 (MQ=255)
              cGCAGATGGCACCCACCCGCCGTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAt                                                                                                                                       >  7:280892/1‑149 (MQ=255)
              cGCAGATGGCACCCACCCGCCGTGCACGGGCGATTTCGTCGTGGGCCACGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAt                                                                                                                                       >  2:238472/1‑149 (MQ=255)
               gCAGATGGCACCCACCCGCCGTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAAtt                                                                                                                                      >  3:211874/1‑149 (MQ=255)
                      gCACCCACCCGCCGTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGtttt                                                                                                                              <  6:339451/149‑1 (MQ=255)
                        acccacccGCCGGGCACGGGCGATTTCTGCGTGGGCCTCGGGCAGGCTGAGGATGCGCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGACACGGGCCCGCGCCAGCCGCCGGAGCACTTCCTGCACGTGAGCGAATTCGGTTGtct                                                                                                                            <  6:132528/149‑1 (MQ=255)
                          ccacccGCCGTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGCATCACTTCCTGCAAGTGATCGAAGTCGGTTTTCTCg                                                                                                                          >  8:305205/1‑149 (MQ=255)
                           cacccGCCGTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCgg                                                                                                                         <  1:111549/149‑1 (MQ=255)
                           cacccGCCGTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCgg                                                                                                                         >  4:210570/1‑149 (MQ=255)
                              ccgccgTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCtt                                                                                                                      <  6:19916/149‑1 (MQ=255)
                               cgccgTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCttt                                                                                                                     >  7:29368/1‑149 (MQ=255)
                                 ccgTGCACGGGGGATTTCGTCGTGGGCCTCGTGCAGGCTGTGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCTACATCGAAACGGGCCAGCGCCAGCCGCCGGATTATTTCCTGCAAGTGATTCAATTCGGTTTTCTCGGCTTTGa                                                                                                                   <  7:59622/149‑1 (MQ=255)
                                 ccgTGCACGGGGGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGGCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGa                                                                                                                   <  7:138143/149‑1 (MQ=255)
                                 ccgTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGa                                                                                                                   <  3:210570/149‑1 (MQ=255)
                                         gggCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTg                                                                                                           >  3:320190/1‑149 (MQ=255)
                                         gggCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTg                                                                                                           >  4:13616/1‑149 (MQ=255)
                                            cGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTgcgc                                                                                                        <  1:330412/149‑1 (MQ=255)
                                               tttCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTg                                                                                                     >  6:143456/1‑149 (MQ=255)
                                                   gtcgtGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCg                                                                                                 <  4:211874/149‑1 (MQ=255)
                                                      gtGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGTAACTTGCGCCTGGCCggg                                                                                               <  4:205359/148‑1 (MQ=255)
                                                      gtGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCTCCTGGCCgggg                                                                                              <  6:7075/149‑1 (MQ=255)
                                                      gtGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCgggg                                                                                              <  2:247056/149‑1 (MQ=255)
                                                      gtGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCgggg                                                                                              <  3:187153/149‑1 (MQ=255)
                                                      gtGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCggg                                                                                               <  1:328356/148‑1 (MQ=255)
                                                      gtGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCggg                                                                                               <  1:238472/148‑1 (MQ=255)
                                                          gCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTAtt                                                                                          <  7:56851/149‑1 (MQ=255)
                                                            cTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGATATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCGGGCCGGGGTATTg                                                                                         >  4:284253/1‑148 (MQ=255)
                                                            cTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTga                                                                                        <  2:317817/149‑1 (MQ=255)
                                                              cGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTg                                                                                         >  5:329597/1‑146 (MQ=255)
                                                               gTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTgaag                                                                                      >  8:299304/1‑148 (MQ=255)
                                                                   aGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACaa                                                                                 <  8:22694/149‑1 (MQ=255)
                                                                     gCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAAt                                                                               >  4:168187/1‑149 (MQ=255)
                                                                     gCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAAt                                                                               <  7:299304/149‑1 (MQ=255)
                                                                           gATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCa                                                                         >  8:373572/1‑149 (MQ=255)
                                                                           gATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCa                                                                         >  5:271470/1‑149 (MQ=255)
                                                                           gATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCAGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCgcgc                                                                          >  5:276331/1‑148 (MQ=255)
                                                                                tctTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTc                                                                     >  3:328201/1‑148 (MQ=255)
                                                                                tctTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTc                                                                     >  1:291893/1‑148 (MQ=255)
                                                                                               gcAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTaccgacc                                                     >  7:3831/1‑149 (MQ=255)
                                                                                                   gTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGt                                                 >  7:335856/1‑149 (MQ=255)
                                                                                                    tGGAAGCCGCCATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGGAAGGGATCGCATTCGGGTTTCTCGGCGTTGAGGGACTTGCGGCTGGGCGGGGTATTGAGGAACGAGTCGGGCACTTCCACGGCGGTACCGCCCggggg                                                >  8:216658/1‑147 (MQ=255)
                                                                                                       aaGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTgcgc                                             >  6:151218/1‑149 (MQ=255)
                                                                                                          ccGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGATTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGcgccgc                                          >  4:16644/1‑149 (MQ=255)
                                                                                                           cGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGcgccgc                                          >  1:289585/1‑149 (MQ=255)
                                                                                                           cGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCgccgcc                                         <  1:143874/149‑1 (MQ=255)
                                                                                                           cGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCgccgcc                                         >  1:276311/1‑149 (MQ=255)
                                                                                                                tCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTg                                    >  3:283051/1‑149 (MQ=255)
                                                                                                                  gAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTtc                                                                                                                             <  4:379759/58‑1 (MQ=255)
                                                                                                                  gAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTtc                                                                                                                             >  3:379759/1‑58 (MQ=255)
                                                                                                                   aaaCGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCa                                  >  1:96659/1‑149 (MQ=255)
                                                                                                                    aaCGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCAcc                                >  3:315073/1‑149 (MQ=255)
                                                                                                                    aaCGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCAcc                                >  5:272425/1‑149 (MQ=255)
                                                                                                                      cGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCcg                              >  7:230060/1‑149 (MQ=255)
                                                                                                                         gccagcgccagcCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCgcgg                           >  7:215128/1‑149 (MQ=255)
                                                                                                                          ccagcgccagcCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCgcggc                          <  6:271470/149‑1 (MQ=255)
                                                                                                                          ccagcgccagcCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCgcggc                          <  6:276331/149‑1 (MQ=255)
                                                                                                                          ccagcgccagcCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCgcggc                          <  4:328201/149‑1 (MQ=255)
                                                                                                                          ccagcgccagcCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCgcggc                          >  5:255146/1‑149 (MQ=255)
                                                                                                                            agcgccagcCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGt                        >  7:243094/1‑149 (MQ=255)
                                                                                                                                  agcCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTc                  <  1:226951/149‑1 (MQ=255)
                                                                                                                                    ccgccgGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTcgc                >  1:273734/1‑149 (MQ=255)
                                                                                                                                     cgccgGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTcgcg               <  8:3831/149‑1 (MQ=255)
                                                                                                                                       ttcgATCACTTGCAGGAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGc             >  1:212185/4‑149 (MQ=255)
                                                                                                                                       ccgGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGc             <  4:315073/149‑1 (MQ=255)
                                                                                                                                          gATCACTTGCAGGAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTAAAGAACAAATCGTGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCtt          >  6:12304/1‑149 (MQ=255)
                                                                                                                                            tCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCg        >  7:38278/1‑149 (MQ=255)
                                                                                                                                            tCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCg        <  4:283051/149‑1 (MQ=255)
                                                                                                                                             cACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCgg       >  6:86044/1‑149 (MQ=255)
                                                                                                                                              aCTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGt      <  2:291893/149‑1 (MQ=255)
                                                                                                                                               cTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGtt     >  8:249797/1‑149 (MQ=255)
                                                                                                                                                  ccTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTcgcg               <  7:214222/136‑1 (MQ=255)
                                                                                                                                                  ccTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTcgcg               >  8:214222/1‑136 (MQ=255)
                                                                                                                                                  ccTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGttt    <  6:255146/147‑1 (MQ=255)
                                                                                                                                                  ccTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGttt    <  3:168187/147‑1 (MQ=255)
                                                                                                                                                  ccTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTc   <  5:143456/148‑1 (MQ=255)
                                                                                                                                                   cTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCt  <  2:96659/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CATGAAGCCCGGGCCGCAGATGGCACCCACCCGCCGTGCACGGGCGATTTCGTCGTGGGCCTCGTGCAGGCTGAGGATGCTCTTGCCATTGTGGCGCAGGTGGAAGCCGACATCGAAACGGGCCAGCGCCAGCCGCCGGATCACTTCCTGCAAGTGATCGAATTCGGTTTTCTCGGCTTTGAGGAACTTGCGCCTGGCCGGGGTATTGAAGAACAAATCGCGCACTTCCACCGAGGTACCGACCGGGTGCGCCGCCGGCTGCACCCGCGGCGTCATGTCGCGGCCTTCGGTTTCT  >  NC_002947/5566811‑5567105

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: