Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,293,717 A→C 5.4% C224G (TGC→GGC)  mmsA‑II ← methylmalonate‑semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,293,7170AC 5.4% 235.6 / 4.4 74C224G (TGC→GGC) mmsA‑IImethylmalonate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (36/34);  new base C (2/2);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

TGGCATCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATG  >  NC_002947/5293573‑5293865
                                                                                                                                                 |                                                                                                                                                    
tGGCATCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATg                                                                                                                                                  >  5:79460/1‑149 (MQ=255)
  gCATCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATgg                                                                                                                                                 <  4:327772/148‑1 (MQ=255)
    aTCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAACCGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCAt                                                                                                                                              >  5:271266/1‑149 (MQ=255)
    aTCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCa                                                                                                                                               >  4:250449/1‑148 (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)
                                 cccccATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAAACGCCTTGATGTCCTG‑GTGGGTGCAAATGGcg                                                                                                                                               >  2:299793/1‑118 (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                                                                              >  8:83101/1‑149 (MQ=255)
                                                                    ctgCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAgcgc                                                                              <  4:307275/149‑1 (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                <  5:48136/138‑1 (MQ=255)
                                                                                                                                              tgCAAATGGTATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGgtgt                >  6:48136/1‑138 (MQ=255)
                                                                                                                                                 aaaTGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATg  >  2:168614/1‑149 (MQ=255)
                                                                                                                                                 aaaTGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATg  <  5:176725/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TGGCATCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATG  >  NC_002947/5293573‑5293865

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