Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,566,416	T→A	10.5%	Q357L (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,416	0	T	A	10.5%	146.5 / 7.8	57	Q357L (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 = 5.52e-01

AGCGGCTGCGAGGGACGTGGGGTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAGCGCCCTTCGCGGAACCGTACTTCATGCTTGGTT  >  NC_002947/5566269‑5566562                                                                                                                                                     |                                                                                                                                                    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                                                                                                                                             >  3:316610/1‑96 (MQ=255)                                                             gAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTg                                                                                                                                             <  4:316610/96‑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                                                            <  2:352730/147‑1 (MQ=255)                                                                                          ttCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGt                                                            <  7:343972/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)                                                                                                                                                     |                                                                                                                                                    AGCGGCTGCGAGGGACGTGGGGTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCA‑CCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCA‑GGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAGCGCCCTTCGCGGAACCGTACTTCATGCTTGGTT  >  NC_002947/5566269‑5566562

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