Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,566,413 T→A 13.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_0029475,566,4130TA13.5% 147.3 / 12.2 58Q358L (CAG→CTG) mutLDNA mismatch repair protein MutL
Reads supporting (aligned to +/- strand):  ref base T (22/28);  new base A (3/5);  total (25/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCGGCTGCGAGGGACGTGGGGTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAGCGCCCTTCGCGGAACCGTACTTCATGCTT  >  NC_002947/5566268‑5566558
                                                                                                                                                 |                                                                                                                                                 
cAGCGGCTGCGAGGGACGTGGGGTGTATTGATACTGGTAACCGGCGCCGCTGCCAGCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCctgct                                                                                                                                                <  2:241432/149‑1 (MQ=255)
cAGCGGCTGCGAGGGACGTGGGGTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCctgct                                                                                                                                                <  8:302676/149‑1 (MQ=255)
      cTGCGAGGGACGTGGGGTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCgg                                                                                                                                          <  5:40071/149‑1 (MQ=255)
      cTGCGAGGGACGTGGGGTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCgg                                                                                                                                          <  1:85842/149‑1 (MQ=255)
         cGAGGGACGTGGGGTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTgg                                                                                                                                       <  3:78433/149‑1 (MQ=255)
         cGAGGGACGTGGGGTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCAGCAGGCCGGt                                                                                                                                         <  6:188775/147‑1 (MQ=255)
         cGAGGGACGTGGGGTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCAGCAGGCCGGTga                                                                                                                                       <  3:115600/149‑2 (MQ=255)
         cGAGGGACGTGGGGTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACCCACCGGCCAGCAGGCCGGt                                                                                                                                         <  5:7057/147‑1 (MQ=255)
         cGAGGGACGGGGTGTGTATGGATACTGGTAACCGGCGCCGCTGACACCATTGGTAAACGACTGTTGCGGGGCTCTAGGCTTTTCAAGCACCGGCGAGGCCAGGCTCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTg                                                                                                                                        <  8:280241/148‑1 (MQ=255)
           aGGGACGTGGGGTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGGTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTgggg                                                                                                                                     <  3:252370/149‑1 (MQ=255)
               aCGTGGGGTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACCCACCGGCCAGCAGGCCGGTgagttcgg                                                                                                                                 <  5:258328/149‑8 (MQ=255)
                     ggTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCt                                                                                                                           <  6:272108/149‑1 (MQ=255)
                      gtgtATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGCTGCGGGGCGCGAGGCTGTTCCAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTc                                                                                                                          <  5:304859/149‑1 (MQ=255)
                       tgtATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCctgct                                                                                                                                                <  7:296234/126‑1 (MQ=255)
                       tgtATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCctgct                                                                                                                                                >  8:296234/1‑126 (MQ=255)
                       tgtATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCt                                                                                                                           >  5:254336/1‑147 (MQ=255)
                         tATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGGGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAgg                                                                                                                       <  4:322193/149‑1 (MQ=255)
                         tATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAgg                                                                                                                       >  5:323918/1‑149 (MQ=255)
                               tACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGgc                                                                                                                 <  1:268847/149‑1 (MQ=255)
                               tACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGgc                                                                                                                 <  3:284963/149‑1 (MQ=255)
                                aCTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACgg                                                                                                                  >  1:39871/1‑147 (MQ=255)
                                     tAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAgg                                                                                                           >  3:279933/1‑149 (MQ=255)
                                                 cTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTg                                                                                               <  6:83813/149‑1 (MQ=255)
                                                    ccaccaTTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTgg                                                                                              <  2:211375/147‑1 (MQ=255)
                                                     caccaTTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGAGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGAGAACAAGCACAGGTACGGCAGCAGGTGCGGCCAGCTGTTCt                                                                                           >  4:232491/1‑149 (MQ=255)
                                                          ttGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCgg                                                                                      <  6:160973/149‑1 (MQ=255)
                                                           tGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGa                                                                                     >  8:60623/1‑149 (MQ=255)
                                                            ggAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGAc                                                                                    >  4:223921/1‑149 (MQ=255)
                                                             gAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTggg                                                                                                                                      >  6:327614/1‑98 (MQ=255)
                                                             gAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTggg                                                                                                                                      <  5:327614/98‑1 (MQ=255)
                                                             gAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGAc                                                                                    <  2:39871/148‑1 (MQ=255)
                                                                   gACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATc                                                                             <  7:112366/149‑1 (MQ=255)
                                                                    aCTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGTCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCg                                                                            <  4:279933/149‑1 (MQ=255)
                                                                        ttGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGcc                                                                         >  3:231578/1‑148 (MQ=255)
                                                                                   gAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACTGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTg                                                              >  6:282299/1‑148 (MQ=255)
                                                                                   gAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGcagcag                                                                                                            >  5:27346/1‑102 (MQ=255)
                                                                                   gAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGcagcag                                                                                                            <  6:27346/102‑1 (MQ=255)
                                                                                      gCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAg                                                          >  8:3070/1‑149 (MQ=255)
                                                                                      gCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACCGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAg                                                          >  1:140087/1‑149 (MQ=255)
                                                                                       cTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACTGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAgg                                                         <  5:282299/149‑1 (MQ=255)
                                                                                             aaGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCa                                                    >  1:281530/1‑148 (MQ=255)
                                                                                              aGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATa                                                  <  2:281530/149‑1 (MQ=255)
                                                                                                cACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACa                                                <  7:60623/149‑1 (MQ=255)
                                                                                                 aCCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAg                                               <  3:223921/149‑1 (MQ=255)
                                                                                                 aCCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAg                                               >  5:307302/1‑149 (MQ=255)
                                                                                                   cGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAgg                                              >  2:340069/1‑148 (MQ=255)
                                                                                                        aGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCa                                        <  7:43747/149‑1 (MQ=255)
                                                                                                          cccAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATg                                      >  4:7005/2‑149 (MQ=255)
                                                                                                               gcgcATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCg                                 >  4:33282/1‑149 (MQ=255)
                                                                                                                cgcATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGa                                <  6:307302/149‑1 (MQ=255)
                                                                                                                cgcATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGa                                <  6:323918/149‑1 (MQ=255)
                                                                                                                  cATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAgc                              <  6:254336/149‑1 (MQ=255)
                                                                                                                          ccTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAGCGCCCTTcg                      >  3:138677/1‑149 (MQ=255)
                                                                                                                            tGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAGCGCCCTTcgcg                    >  1:282553/1‑149 (MQ=255)
                                                                                                                                       cccACCGGCCAGCAGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATg                                      <  8:3782/118‑1 (MQ=255)
                                                                                                                                       cccACCGGCCAGCAGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATg                                      >  7:3782/3‑120 (MQ=255)
                                                                                                                                         cACCGGCCAGCAGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTTCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAGCGCCCTTCGCGGAACCGTACTTCa       >  3:285926/1‑149 (MQ=255)
                                                                                                                                         cACCGGCCAGCAGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTGCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCCCCGAGCGCCCTTCGCGGAACCGTACTTCa       >  3:281245/1‑149 (MQ=255)
                                                                                                                                          aCCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCa                                                    >  3:252450/1‑103 (MQ=255)
                                                                                                                                          aCCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCa                                                    <  4:252450/103‑1 (MQ=255)
                                                                                                                                          aCCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAGCGCCCTTCGCGGAACCGTACTTCAt      >  4:175180/1‑149 (MQ=255)
                                                                                                                                              gCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACCGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAGCGCCCTTCGCGGAACCGTACTTCATGCtt  <  2:140087/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CAGCGGCTGCGAGGGACGTGGGGTGTATTGATACTGGTAACCGGCGCCGCTGCCACCATTGGAAAACGACTGTTGCGGGGCTCGAGGCTGTTCAAGCACCGGCGAGGCCAGGCGCATTTCACCCTGGGGCCCGAACTCACCGGCCTGCTGGCCGGTGGGGCGAACCAGCTCAGGTACGGCAGCAGGTGCGGCCAGCTGGTCTTCCGGACGTACATCGGCCAAGGCACGGTGCAAGGTGCCATACAGGAAGTCATGCACCGAGCGCCCTTCGCGGAACCGTACTTCATGCTT  >  NC_002947/5566268‑5566558

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