Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,566,425 G→A 9.1% P354L (CCC→CTC)  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,4250GA 9.1% 149.4 / 3.9 54P354L (CCC→CTC) mutLDNA mismatch repair protein MutL
Reads supporting (aligned to +/- strand):  ref base G (23/26);  new base A (3/2);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GAGGGACGTGGGGTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAGCGCCCTTCGCGGAACCGTACTTCATGCTTGGTTGGGTGCACGT  >  NC_002947/5566278‑5566572
                                                                                                                                                    |                                                                                                                                                    
gAGGGACGTGGGGTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTggg                                                                                                                                                     >  8:132554/1‑149 (MQ=255)
 aGGGACGTGGGGTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTgggg                                                                                                                                                    <  5:207368/149‑1 (MQ=255)
     aCGTGGGGTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGaa                                                                                                                                                <  5:125339/149‑1 (MQ=255)
     aCGTGGGGTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGaa                                                                                                                                                <  2:368409/149‑1 (MQ=255)
     aCGTGGGGTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGaa                                                                                                                                                <  3:205750/149‑1 (MQ=255)
          gggTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGc                                                                                                                                           >  7:153765/1‑149 (MQ=255)
            gtgtATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTc                                                                                                                                         <  8:345765/149‑1 (MQ=255)
             tgtATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCa                                                                                                                                        <  7:339789/149‑1 (MQ=255)
               tATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑gg                                                                                                                                     >  3:87254/1‑149 (MQ=255)
               tATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGACGGTGGGGCGAACCAGCTCA‑gg                                                                                                                                     >  4:292222/1‑149 (MQ=255)
                  tGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTAc                                                                                                                                  >  2:231353/1‑149 (MQ=255)
                    aTACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTgggttcgg                                                                                                                                                >  1:322316/1‑129 (MQ=255)
                    aTACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTgggttcgg                                                                                                                                                <  2:322316/134‑6 (MQ=255)
                    aTACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACgg                                                                                                                                >  7:34931/1‑149 (MQ=255)
                     tACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACgg                                                                                                                                <  8:34931/148‑1 (MQ=255)
                               cGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCg                                                                                                                     <  7:132554/149‑1 (MQ=255)
                                    ccgcTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAg                                                                                                                >  4:296957/1‑149 (MQ=255)
                                         gCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGt                                                                                                           <  8:153765/149‑1 (MQ=255)
                                           caccaTTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCt                                                                                                         >  3:221324/1‑149 (MQ=255)
                                              caTTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTcc                                                                                                      <  4:221324/149‑1 (MQ=255)
                                                   gAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACg                                                                                                 >  2:179910/1‑149 (MQ=255)
                                                    aaaaCGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACg                                                                                                 >  3:57237/1‑148 (MQ=255)
                                                    aaaaCGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACg                                                                                                 >  1:352730/1‑148 (MQ=255)
                                                          aCTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCg                                                                                          <  2:60011/149‑1 (MQ=255)
                                                              ttGCGGGGCTCGAGGCTGTTCAAGCAGCCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGc                                                                                        >  1:201442/1‑148 (MQ=255)
                                                              ttGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGc                                                                                        <  5:226456/147‑1 (MQ=255)
                                                                gCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAg                                                                                    >  4:615/1‑149 (MQ=255)
                                                                 cGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAg                                                                                    <  5:270965/148‑1 (MQ=255)
                                                                       tCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGt                                                                             >  4:8126/1‑149 (MQ=255)
                                                                         gAGGATGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGCACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCACGTGCGGCCAGCTGGTCTACCGGACGTACATCGGCCAACGCACGGTg                                                                            >  6:17646/1‑148 (MQ=255)
                                                                              tGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGt                                                                      <  4:57237/149‑1 (MQ=255)
                                                                                ttCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGt                                                                      <  7:343972/147‑1 (MQ=255)
                                                                                ttCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGt                                                                      <  2:352730/147‑1 (MQ=255)
                                                                                ttCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGc                                                                    <  1:231353/149‑1 (MQ=255)
                                                                                   aaGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCa                                                                  >  6:339820/1‑148 (MQ=255)
                                                                                     gCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACCTTGCACCGTGccttgg                                                               <  5:254923/149‑5 (MQ=255)
                                                                                     gCA‑CCGCCGAGCCCAGGCGCATATCACCCTGGGGCCCCAACTCACCGGCATGCTGGCCGGTCGGGCGAACCAGCTCA‑GGTACGGCACCAGGAGCGCCCACCTGGTCTTCCGGACCAACATCCGCCAAGGAACGCTGCAACGAGCCaaa                                                                >  8:259335/1‑146 (MQ=255)
                                                                                      cA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATa                                                                >  6:6964/1‑147 (MQ=255)
                                                                                           ggCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGaa                                                          <  4:87254/149‑1 (MQ=255)
                                                                                               aGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCa                                                      <  3:8126/149‑1 (MQ=255)
                                                                                                 gCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATg                                                    >  4:368840/1‑149 (MQ=255)
                                                                                                                  cTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAGCGCCCTTcgc                                   <  2:201442/149‑1 (MQ=255)
                                                                                                                     gggCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAGCGCCCTTcgcg                                  <  1:179910/147‑1 (MQ=255)
                                                                                                                      ggCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAGCGCCCTTCGCGGaa                               <  5:6964/149‑1 (MQ=255)
                                                                                                                      ggCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAGCGCCCTTCGCGGaa                               <  3:296957/149‑1 (MQ=255)
                                                                                                                           cgccccACCGGCCAGCAGGCCGGTGAGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAGCGCCCTTCGCGGa                                >  4:332081/6‑143 (MQ=255)
                                                                                                                           cgccccACCGGCCAGCAGGCCGGTGAGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAGCGCCCTTCGCGGa                                <  3:332081/138‑1 (MQ=255)
                                                                                                                            aaCTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAGCGCCCTTCGCGGAACCGTAc                         <  3:292222/149‑1 (MQ=255)
                                                                                                                             ccccACCGGCCAGCAGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGACTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAGCGCCCTTCGCGGAACCGTACt                        >  8:75713/4‑149 (MQ=255)
                                                                                                                             ccccACCGGCCAGCAGGCCGGTGAGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAgcgc                                          >  2:163984/4‑131 (MQ=255)
                                                                                                                             ccccACCGGCCAGCAGGCCGGTGAGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAgcgc                                          <  1:163984/128‑1 (MQ=255)
                                                                                                                               tCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAGCGCCCTTCGCGGAACCGTACtt                       >  5:79790/1‑148 (MQ=255)
                                                                                                                                cACCGGCCAGCAGGCCGGTGAGGCGAACCAGCTCAGGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAGCGCCCTTCGCGGAACCGTACtt                       >  1:69168/1‑148 (MQ=255)
                                                                                                                                 aCCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAGCGCCCTTCGCGGAACCGTACTTCAt                    >  3:284780/1‑149 (MQ=255)
                                                                                                                                        tgctgGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAGCGCCCTTCGCGGAACCGTACTTCATGCttggt             >  3:376706/1‑149 (MQ=255)
                                                                                                                                           tgGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAGCGCCCTTCGCGGAACCGTACTTCATGCttggtt            >  4:139155/1‑147 (MQ=255)
                                                                                                                                           tgGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAGCGCCCTTCGCGGAACCGTACTTCATGCttggtt            <  3:139155/147‑1 (MQ=255)
                                                                                                                                             gCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAGCGCCCTTCGCGGAACCGTACTTCATGCTTGGTTGGGt        <  6:281263/149‑1 (MQ=255)
                                                                                                                                                ggtggtGCGAACCAGCTCA‑GGGACGGCAGCAGGGGCGGCCAGCTGGGCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAGCGCCCTTCGCGGAACCGTACTTCATGCTTGGTTGGGTGCa     <  6:116522/149‑1 (MQ=255)
                                                                                                                                                    gggCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACAGAGCGCCCTTAGCTGAACCGTACTTCATTCTTGTTTGGGTGCACGt  >  4:239447/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAGGGACGTGGGGTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAGCGCCCTTCGCGGAACCGTACTTCATGCTTGGTTGGGTGCACGT  >  NC_002947/5566278‑5566572

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: ‑