Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,293,719 A→T 5.4% I223N (ATT→AAT)  mmsA‑II ← methylmalonate‑semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,293,7190AT 5.4% 233.6 / 2.3 74I223N (ATT→AAT) mmsA‑IImethylmalonate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (36/34);  new base T (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 = 8.82e-01

TGGCATCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATGGC  >  NC_002947/5293573‑5293867
                                                                                                                                                   |                                                                                                                                                    
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                              <  3:262756/149‑1 (MQ=255)
                                                                                                                      aaTCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGa                              <  6:218265/149‑1 (MQ=255)
                                                                                                                      aaTCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGa                              <  1:291991/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)
                                                                                                                                                   |                                                                                                                                                    
TGGCATCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTG‑GTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATGGC  >  NC_002947/5293573‑5293867

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