Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,566,413	T→A	10.5%	Q358L (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_002947	5,566,413	0	T	A	10.5%	153.2 / 8.3	57	Q358L (CAG→CTG)	mutL	DNA mismatch repair protein MutL
Reads supporting (aligned to +/- strand):  ref base T (23/28);  new base A (4/2);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.48e-01

AGCGGCTGCGAGGGACGTGGGGTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAGCGCCCTTCGCGGAACCGTACTTCATGCTTGGT  >  NC_002947/5566269‑5566561                                                                                                                                                  |                                                                                                                                                      aGCGGCTGCGAGGGACGTGGGGTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCctgctg                                                                                                                                                   <  7:309905/149‑1 (MQ=255)     gCTGCGAGGGACGTGGGGTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCg                                                                                                                                               <  5:228239/149‑1 (MQ=255)          gAGGGACGTGGGGTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTggg                                                                                                                                          >  8:132554/1‑149 (MQ=255)           aGGGACGTGGGGTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTgggg                                                                                                                                         <  5:207368/149‑1 (MQ=255)               aCGTGGGGTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGaa                                                                                                                                     <  3:205750/149‑1 (MQ=255)               aCGTGGGGTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGaa                                                                                                                                     <  2:368409/149‑1 (MQ=255)               aCGTGGGGTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGaa                                                                                                                                     <  5:125339/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                                                                                                                                     <  2:322316/134‑6 (MQ=255)                              aTACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTgggttcgg                                                                                                                                     >  1:322316/1‑129 (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‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTg                                                                                                                                            <  4:316610/96‑1 (MQ=255)                                                             gAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTg                                                                                                                                            >  3:316610/1‑96 (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)                                                                                                                                                  |                                                                                                                                                      AGCGGCTGCGAGGGACGTGGGGTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAGCGCCCTTCGCGGAACCGTACTTCATGCTTGGT  >  NC_002947/5566269‑5566561

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