Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,293,713 G→A 12.3% T225I (ACC→ATC)  mmsA‑II ← methylmalonate‑semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,293,7130GA12.3% 168.2 / 12.6 65T225I (ACC→ATC) mmsA‑IImethylmalonate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (31/26);  new base A (4/4);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01

GTGCGGTTGGCATCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCC  >  NC_002947/5293566‑5293855
                                                                                                                                                   |                                                                                                                                              
gtgcgGTTGGCATCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTggtgggt                                                                                                                                               <  7:306060/149‑1 (MQ=255)
     gTTGGCATCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGTTAACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGATGccat                                                                                                                                          <  1:307130/149‑4 (MQ=255)
        ggCATCAGGCAGCACCACCTGGTGGTTCTTGGCCCCCATCATTGACTGCACACGGTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATg                                                                                                                                        <  5:120461/148‑1 (MQ=255)
        ggCATCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATg                                                                                                                                        <  2:139752/148‑1 (MQ=255)
         gCATCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGATGccatttg                                                                                                                                       >  3:243565/1‑142 (MQ=255)
          cATCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCa                                                                                                                                     >  3:253033/1‑149 (MQ=255)
           aTCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGCGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCAt                                                                                                                                    >  8:272236/1‑149 (MQ=255)
           aTCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCAt                                                                                                                                    <  3:258037/149‑1 (MQ=255)
            tCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATc                                                                                                                                   <  8:72011/149‑1 (MQ=255)
             cAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATc                                                                                                                                   >  8:274583/1‑148 (MQ=255)
               ggCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATccac                                                                                                                                <  4:253033/149‑1 (MQ=255)
                 cagcaCCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGATGCCATTTGCACccacca                                                                                                                              <  4:243565/149‑1 (MQ=255)
                   gcaCCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCGAGCTGGCTGCCGAGGTTGTAAACGAGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCaccacc                                                                                                                            >  1:279242/1‑149 (MQ=255)
                    caccacCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCt                                                                                                                           >  4:242393/1‑149 (MQ=255)
                     accaccGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGATGCCATTTGCACccacca                                                                                                                              <  4:50694/145‑1 (MQ=255)
                     accaccGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGATGCCATTTGCACccacca                                                                                                                              >  3:50694/1‑145 (MQ=255)
                              tGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTgccgcc                                                                                                                 >  5:94497/1‑149 (MQ=255)
                                        cccccATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACg                                                                                                       <  7:272236/149‑1 (MQ=255)
                                          cccATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGtt                                                                                                     >  6:21831/1‑149 (MQ=255)
                                           ccATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTc                                                                                                    <  4:318918/149‑1 (MQ=255)
                                                      tGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACgccggc                                                                                         >  5:168788/1‑149 (MQ=255)
                                                      tGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACgccggc                                                                                         >  8:240707/1‑149 (MQ=255)
                                                      tGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGAGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACgccggc                                                                                         >  4:326354/1‑149 (MQ=255)
                                                        cacaCGCTTGCCATGCTGGCTGCCGATGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCGAATGGCATCCACCACCTGCTTgccgc                                                                                                                  <  2:35278/122‑1 (MQ=255)
                                                        cacaCGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACgccggccg                                                                                       >  6:73795/1‑149 (MQ=255)
                                                          cacGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGCCTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGGTCAGCACGccggccgg                                                                                      >  8:283297/1‑148 (MQ=255)
                                                          cacGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGAGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGt                                                                                     <  2:279242/149‑1 (MQ=255)
                                                           acGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTa                                                                                    >  1:26755/1‑149 (MQ=255)
                                                             gCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACg                                                                                  >  8:43108/1‑149 (MQ=255)
                                                                   cATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCtt                                                                           >  5:87595/1‑149 (MQ=255)
                                                                    aTGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATccac                                                                                                                                >  4:149563/1‑96 (MQ=255)
                                                                    aTGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATccac                                                                                                                                <  3:149563/96‑1 (MQ=255)
                                                                     tGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTc                                                                          >  2:46584/1‑149 (MQ=255)
                                                                       ctggctgCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGCCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTc                                                                          >  1:3233/1‑147 (MQ=255)
                                                                           ctgCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAgcgc                                                                    <  5:73795/149‑1 (MQ=255)
                                                                                     tGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACTACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGa                                                            >  5:335820/1‑147 (MQ=255)
                                                                                     tGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACTACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGa                                                            <  6:335820/147‑1 (MQ=255)
                                                                                      gTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCATATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGAc                                                           >  4:34454/1‑147 (MQ=255)
                                                                                      gTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCTGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGcc                                                                                >  6:99787/1‑126 (MQ=255)
                                                                                      gTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCTGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGcc                                                                                <  5:99787/126‑1 (MQ=255)
                                                                                      gTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGAcc                                                          <  5:21831/148‑1 (MQ=255)
                                                                                      gTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGAcc                                                          <  7:297854/148‑1 (MQ=255)
                                                                                          aCGTGCGTACCGACTTCAGTAGAGCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGc                                                    <  6:87595/149‑1 (MQ=255)
                                                                                                    cGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACgcc                                                                                            <  8:329207/100‑1 (MQ=255)
                                                                                                    cGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACgcc                                                                                            >  7:329207/1‑100 (MQ=255)
                                                                                                     gACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACa                                          <  7:274583/149‑1 (MQ=255)
                                                                                                     gACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACa                                          <  7:240707/149‑1 (MQ=255)
                                                                                                        ttCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCg                                       >  1:213672/1‑149 (MQ=255)
                                                                                                          cAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGccggccgg                                                                                      >  6:45215/1‑100 (MQ=255)
                                                                                                          cAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGccggccgg                                                                                      <  5:45215/100‑1 (MQ=255)
                                                                                                            gTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATc                                   <  3:242393/149‑1 (MQ=255)
                                                                                                               gagCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCa                                                   >  8:342835/1‑130 (MQ=255)
                                                                                                               gagCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCa                                                   <  7:342835/130‑1 (MQ=255)
                                                                                                               gagCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTg                                <  2:26755/149‑1 (MQ=255)
                                                                                                                   ccACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCg                            >  2:251232/1‑149 (MQ=255)
                                                                                                                       gAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATg                        >  7:274731/1‑149 (MQ=255)
                                                                                                                           gAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGttt                    >  6:374411/1‑149 (MQ=255)
                                                                                                                           gAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGttt                    >  3:98392/1‑149 (MQ=255)
                                                                                                                             aaTCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGa                  <  1:46584/149‑1 (MQ=255)
                                                                                                                               tCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGc                <  6:168788/149‑1 (MQ=255)
                                                                                                                                   cTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACg             >  3:102738/1‑148 (MQ=255)
                                                                                                                                         caggtggtggATGCCATTTGCACCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGgtgt      >  6:58516/5‑149 (MQ=39)
                                                                                                                                            cTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGc   <  2:3233/149‑1 (MQ=255)
                                                                                                                                            cTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGc   >  5:242953/1‑149 (MQ=255)
                                                                                                                                             tggtggATGCCATTTGCACCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCg                                                             >  8:20209/1‑90 (MQ=25)
                                                                                                                                             tggtggATGCCATTTGCACCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCg                                                             <  7:20209/90‑1 (MQ=38)
                                                                                                                                               gagggtgCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGcc  <  3:326354/145‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
GTGCGGTTGGCATCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCC  >  NC_002947/5293566‑5293855

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: