Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,293,725 T→C 5.4% D221G (GAT→GGT)  mmsA‑II ← methylmalonate‑semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,293,7250TC 5.4% 219.8 / 2.2 74D221G (GAT→GGT) mmsA‑IImethylmalonate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (34/36);  new base C (2/2);  total (36/38)
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.73e-01

ATCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATGGCCATGGGG  >  NC_002947/5293577‑5293874
                                                                                                                                                     |                                                                                                                                                     
aTCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAACCGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCAt                                                                                                                                                       >  5:271266/1‑149 (MQ=255)
 tCATGCAGCACCACCGCGTTGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATc                                                                                                                                                      <  3:248505/149‑1 (MQ=255)
      cagcaCCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATccacca                                                                                                                                                 >  2:225678/1‑149 (MQ=255)
      cagcaCCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGGACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGATGCCATTTGCACccacca                                                                                                                                                 <  5:235003/149‑1 (MQ=255)
      cagcaCCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGATGGCCAAGCTGGCGGCCGAGGTTGTAAACGTGCGTACCGACTTCATTAGAGCCCACGAAGTAAATCGCCTTGATGTCCTG‑GTGGATGCCATTTGCACccacca                                                                                                                                                 <  5:235005/149‑1 (MQ=255)
       agcaCCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCcaccac                                                                                                                                                <  6:79460/149‑1 (MQ=255)
        gcaCCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCaccacc                                                                                                                                               >  2:206587/1‑149 (MQ=255)
          accaccGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATccac                                                                                                                                                   <  5:311924/143‑1 (MQ=255)
          accaccGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATccac                                                                                                                                                   >  6:311924/1‑143 (MQ=255)
             accGCGTGGTTCGTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCtt                                                                                                                                          <  7:232137/149‑1 (MQ=255)
               cgcgTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGACGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTgc                                                                                                                                        <  7:274054/149‑1 (MQ=255)
                  gTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTgccgc                                                                                                                                     >  3:313494/1‑149 (MQ=255)
                     gTTCTTGGCCCCCATCATTGACTGCGCACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGt                                                                                                                                  >  4:262756/1‑149 (MQ=255)
                     gTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGccgccg                                                                                                                                   <  5:37562/148‑1 (MQ=255)
                           ggCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGcacca                                                                                                                            >  2:217217/1‑149 (MQ=255)
                           ggCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGca                                                                                                                               <  6:271266/146‑1 (MQ=255)
                             cccccATCATTTACTTCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGcaccac                                                                                                                           <  3:298995/148‑1 (MQ=255)
                               cccATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGtt                                                                                                                        >  2:21195/1‑149 (MQ=255)
                               cccATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGtt                                                                                                                        <  1:10734/149‑1 (MQ=255)
                                ccATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGAGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGtt                                                                                                                        <  1:79655/149‑1 (MQ=255)
                                ccATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTc                                                                                                                       <  5:28563/149‑1 (MQ=255)
                                 catcatTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTCCa                                                                                                                      <  1:21195/149‑1 (MQ=255)
                                   tcatTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGc                                                                                                                    <  6:46452/149‑1 (MQ=255)
                                      ttGACTGCCCACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCAACGAAGGAAATCGCCTTCATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCATCCCGCCCAGCACCACGATCAGCACg                                                                                                                 >  8:254435/1‑149 (MQ=255)
                                            gCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACgccggcc                                                                                                           >  8:115991/1‑149 (MQ=255)
                                                 cGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTAc                                                                                                      <  1:206587/149‑1 (MQ=255)
                                                         aTGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCGGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCtt                                                                                              >  5:294294/1‑149 (MQ=255)
                                                             tggctgCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAg                                                                                          >  3:307275/1‑149 (MQ=255)
                                                               gctgCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAgcg                                                                                        >  6:176725/1‑149 (MQ=255)
                                                                ctgCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAgcgc                                                                                       <  4:307275/149‑1 (MQ=255)
                                                                ctgCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAgcgc                                                                                       >  8:83101/1‑149 (MQ=255)
                                                                  gCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAgcgc                                                                                       >  6:75825/1‑147 (MQ=255)
                                                                          tGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGAc                                                                              >  4:314870/1‑148 (MQ=255)
                                                                           gTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATACACCACCTGCTTGCAGCCGTGCACCACGTTAAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGAc                                                                              >  7:265979/1‑147 (MQ=255)
                                                                            tAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAg                                                                           <  6:294294/149‑1 (MQ=255)
                                                                             aaaCGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCAACACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAg                                                                           >  4:58/1‑148 (MQ=255)
                                                                              aaCGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAg                                                                           >  2:270413/1‑147 (MQ=255)
                                                                              aaCGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGaa                                                                         <  1:270413/149‑1 (MQ=255)
                                                                                     gTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGt                                                                  <  4:313494/149‑1 (MQ=255)
                                                                                        ccGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGa                                                               >  5:218265/1‑149 (MQ=255)
                                                                                           aCTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAg                                                            >  4:46129/1‑149 (MQ=255)
                                                                                            cTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGc                                                           <  3:250449/149‑1 (MQ=255)
                                                                                                aGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGAt                                                       >  2:224310/1‑149 (MQ=255)
                                                                                                aGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGAt                                                       <  6:136601/149‑1 (MQ=255)
                                                                                                     agCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGc                                                  <  7:83101/149‑1 (MQ=255)
                                                                                                         cACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCg                                               >  6:40343/1‑148 (MQ=255)
                                                                                                         cACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCg                                               >  8:317399/1‑148 (MQ=255)
                                                                                                          aCGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGa                                             <  3:46129/149‑1 (MQ=255)
                                                                                                                gAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGttt                                       <  7:115991/149‑1 (MQ=255)
                                                                                                                 aaaTCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGttt                                       >  7:66232/1‑148 (MQ=255)
                                                                                                                 aaaTCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTg                                      <  5:75825/149‑1 (MQ=255)
                                                                                                                  aaTCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGa                                     <  1:291991/149‑1 (MQ=255)
                                                                                                                  aaTCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGa                                     <  6:218265/149‑1 (MQ=255)
                                                                                                                  aaTCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGa                                     <  3:262756/149‑1 (MQ=255)
                                                                                                                       ccTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGCCAGCGGATCCTGCTCGGATGGTTTGAGCAc                                 >  7:305829/1‑148 (MQ=255)
                                                                                                                           gATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAgg                            >  4:16249/1‑149 (MQ=255)
                                                                                                                             tGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGa                                     <  4:22445/138‑1 (MQ=255)
                                                                                                                             tGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGgtg                          >  5:273059/1‑149 (MQ=255)
                                                                                                                             tGTACTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGa                                     >  3:22445/1‑138 (MQ=255)
                                                                                                                                ggtg‑gtggATGCCATTTGCACCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTg                       >  4:194781/3‑149 (MQ=255)
                                                                                                                                 gtg‑gtggATGCCATTTGCACCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGc                      >  5:238455/2‑149 (MQ=255)
                                                                                                                                  tg‑gtggGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTc                                                <  2:287152/122‑1 (MQ=255)
                                                                                                                                  tg‑gtggGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGcc                     >  6:201213/1‑149 (MQ=255)
                                                                                                                                   g‑gtgggtgCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCg                    >  3:265516/1‑149 (MQ=255)
                                                                                                                                     gtgggtgCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGcc                     <  8:66232/147‑1 (MQ=255)
                                                                                                                                     gtgggtgCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCg                    <  1:217217/148‑1 (MQ=255)
                                                                                                                                      tgggtgCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCg                    >  5:323565/1‑147 (MQ=255)
                                                                                                                                         gtgCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGAc               <  3:16249/149‑1 (MQ=255)
                                                                                                                                          tgCAAATGGTATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGgtgt                         >  6:48136/1‑138 (MQ=255)
                                                                                                                                          tgCAAATGGTATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGgtgt                         <  5:48136/138‑1 (MQ=255)
                                                                                                                                             aaaTGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATg           >  2:168614/1‑149 (MQ=255)
                                                                                                                                             aaaTGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATg           <  5:176725/149‑1 (MQ=255)
                                                                                                                                               aTGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATGGc         <  5:40343/149‑1 (MQ=255)
                                                                                                                                                   cATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCACCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATGGCCATg     >  6:47476/1‑149 (MQ=255)
                                                                                                                                                    aTCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATGGCCATgg    <  1:168614/149‑1 (MQ=255)
                                                                                                                                                     tCCACCACCTGCTTGGCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCATAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATGGCCATggg   <  7:297315/149‑1 (MQ=255)
                                                                                                                                                     tCCACCACCTGCTTGCCGCCGTGCACCACGTTCATCACTCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGTAGCATGGTCGACAGCGGATCCTGCTGGGATGGTTTGAGCACGAAGGTGTTTCCGCAGACGATGGCCATggg   <  7:254435/149‑1 (MQ=255)
                                                                                                                                                     tCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATGGCCATggg   >  8:36631/1‑149 (MQ=255)
                                                                                                                                                     tCCACAACCTGCTTCCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCTCCTTCCAGCGCCAGCTCCACCACAAGCAT‑GTCGACAGCGGATCCTGCTCGCATGCTTTGAGCACGAAGGAGTAGCCGCAGACGATGGCCAAgggc  >  2:316561/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATGGCCATGGGG  >  NC_002947/5293577‑5293874

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