Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,729,576 T→A 6.7% D87V (GAT→GTT)  sucC ← bifunctional succinyl‑CoA synthetase subunit beta/glutaryl‑CoA synthetase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,729,5760TA 6.7% 164.1 / 3.2 60D87V (GAT→GTT) sucCbifunctional succinyl‑CoA synthetase subunit beta/glutaryl‑CoA synthetase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (28/28);  new base A (2/2);  total (30/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.69e-01

GACTTTCTCGATGTCCACGCCACCTTCGGTGGAAGCCATGAACACGATACGGCGGCTCGAACGATCGACTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTGGCGTCTTCCTTGCTGCGAACCAGCTTTACGCCGCCCGCTTTACCGCGACCACCCGCGTGAACCTGGGCTTTTACAACCCACTCGTTCCCGCCGA  >  NC_002947/4729428‑4729723
                                                                                                                                                    |                                                                                                                                                   
gACTTTCTCGATGTCCACGCCACCTTCGGTGGAAGCCATGAACACGATACGGCGGCTCGAACGATCGACTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGCAGCTTGGAGACGGGTTGACCGTTGGCAt                                                                                                                                                     >  3:204981/1‑149 (MQ=255)
  cTTTCTCGATGTCCACGCCACCTTCGGTGGAAGCCATGAACACGATACGGCGGCTCGAACGATCGACTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCg                                                                                                                                                   <  5:308113/149‑1 (MQ=255)
  cTTTCTCGATGTCCACGCCACCTTCGGTGGAAGCCATGAACACGATACGGCGGCTCGAACGATCGACTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCg                                                                                                                                                   <  3:254926/149‑1 (MQ=255)
     tctcGATGTCCACGCCACCTTCGGTGGAAGCCATGAACACGATACGGCGGCTCGAACGATCGACTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGATGCCAACGGt                                                                                                                                                 >  8:222619/1‑148 (MQ=255)
       tcGATGTCCACGCCACCTTCGGTGGAAGCCATGAACACGATACGGCGGCTCGAACGATCGACTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATc                                                                                                                                                    >  8:265599/1‑143 (MQ=255)
       tcGATGTCCACGCCACCTTCGGTGGAAGCCATGAACACGATACGGCGGCTCGAACGATCGACTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATc                                                                                                                                                    <  7:265599/143‑1 (MQ=255)
         gATGTCCACGCCACCTTCGGTGGAAGCCATGAACACGATACGGCGGCTCGAACGATCGACTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTggt                                                                                                                                            <  1:193068/149‑1 (MQ=255)
         gATGTCCACGCCACCTTCGGTGGAAGCCATGAACACGATACGGCGGCTCGAACGATCGACTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTggt                                                                                                                                            <  7:118642/149‑1 (MQ=255)
         gATGTCCACGCCACCTTCGGTGGAAGCCATGAACACGATACGGCGGCTCGAACGATCGACTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTgg                                                                                                                                             >  2:139540/1‑148 (MQ=255)
         gATGTCCACGCCACCTTCGGTGGAAGCCATGAACACGATACGGCGGCTCGAACGATCGACTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTg                                                                                                                                              <  2:372225/147‑1 (MQ=255)
         gATGTCCACGCCACCTTCGGTGGAAGCCATGAACACGATACGGCGGCTCGAACGATCGACTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGATGCCAACGGTcaacc                                                                                                                                            <  7:222619/149‑5 (MQ=255)
                aCGCCACCTTCGGTGGAAGCCATGAACACGATACGGCGGCTCGAACGATCGACTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAAc                                                                                                                                     >  3:150651/1‑149 (MQ=255)
                    cACCTTCGGTGGAAGCCATGAACACGATACGGCGGCTCGAACGATCGACTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAAc                                                                                                                                     >  2:297684/1‑145 (MQ=255)
                       cTTCGGTGGAAGCCATGAACACGATACGGCGTCTCGAACGATCGACTACAGCGCCCAGGTACAGCTCTTTGTCGATTTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTc                                                                                                                              <  5:159/149‑1 (MQ=255)
                       cTTCGGTGGAAGCCATGAACACGATACGGCGGCTCGAACGATCGACTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTc                                                                                                                              >  5:97390/1‑149 (MQ=255)
                         tCGGTGGAAGCCATGAACACGATACGGCGGCTCGAACGATCGACTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCtt                                                                                                                            <  3:7483/149‑1 (MQ=255)
                           ggTGGAAGCCATGAACACGATACGGCGGCTCGAACGATCGACTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCtt                                                                                                                            >  4:7483/1‑147 (MQ=255)
                           ggTGGAAGCCATGAACACGATACGGCGGCTCGAACGATCGACTACAGCGCCCAGGTACAGCGCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTg                                                                                                                           <  8:102226/148‑1 (MQ=255)
                            gTGGAAGCCATGAACACGATACGGCGGCTCGAACGATCGACTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGcc                                                                                                                         >  6:116894/1‑149 (MQ=255)
                               gAAGCCATGAACACGATACGGCGGCTCGAACGATCGACTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCaa                                                                                                                      >  7:55954/1‑149 (MQ=255)
                                    cATGAACACGATACGGCGGCTCGAACGATCGACTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTggt                                                                                                                                            >  2:348394/1‑122 (MQ=255)
                                    cATGAACACGATACGGCGGCTCGAACGATCGACTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTggt                                                                                                                                            <  1:348394/122‑1 (MQ=255)
                                     aTGAACACGATACGGCGGCTCGAACGATCGACTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCCAGCGGCAGGGTGGAGCCGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACtg                                                                                                                >  1:127241/1‑149 (MQ=255)
                                      ccaaCACGATACGGCGGCTCGAACGATCGACTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACtgt                                                                                                               <  1:139540/147‑1 (MQ=255)
                                         acacGATACGGCGGCTCGAACGATCGACTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACtgtg                                                                                                              >  2:169042/1‑147 (MQ=255)
                                         acacGATACGGCGGCTCGAACGATCGACTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCa                                                                                                            <  6:211395/149‑1 (MQ=255)
                                           acGATACGGCGGCTCGAACGATCGACTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTATTGACCAACAACTGTTCAGc                                                                                                          >  4:16605/1‑149 (MQ=255)
                                                     ggcTCGACCGATCGACTCCAGCGCCCAGGGACAGCTCTTTGGCGACGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTg                                                                                                <  4:336763/149‑1 (MQ=255)
                                                          gAACGATCGACTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCtt                                                                                                                            >  4:326195/1‑116 (MQ=255)
                                                          gAACGATCGACTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTGGCGTc                                                                                           <  1:297684/149‑1 (MQ=255)
                                                             cGATCGACTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTGGCGTCTTc                                                                                        <  5:116894/149‑1 (MQ=255)
                                                                   aCTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTGGCGTCTTCCTtgct                                                                                  >  3:311701/1‑149 (MQ=255)
                                                                       cAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTGGCGTCTTCCTTGCTGCg                                                                               >  4:189929/1‑148 (MQ=255)
                                                                       cAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTGGCGTCTTCCTTGCTGCGa                                                                              <  2:190510/149‑1 (MQ=255)
                                                                       cAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATGCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCATCGAACGCCTTGGCGTCTTCCTtgctgc                                                                                <  6:182111/147‑1 (MQ=255)
                                                                         gcgcCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTGGCGTCTTCCTTGCTGCGAAc                                                                            >  3:189695/1‑149 (MQ=255)
                                                                             ccAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTGGCGTCTTCCTTGCTGCGAACCAGc                                                                        >  3:179351/1‑149 (MQ=255)
                                                                                     aGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTGGCGTCTTCCTTGCTGCGAACCAGCTTTAcgcc                                                                >  1:358349/1‑149 (MQ=255)
                                                                                          tttGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTCGGCGTCTTCCTtgctgc                                                                                >  8:10715/1‑128 (MQ=255)
                                                                                          tttGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTCGGCGTCTTCCTtgctgc                                                                                <  7:10715/128‑1 (MQ=255)
                                                                                          tttGGCGATGGCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTGGCGTCTTCCTTGCTGCGAACCAGCTTTACgccgcc                                                             <  8:64475/147‑1 (MQ=255)
                                                                                              gCGAAGTGAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTCGGTATCCAGATTCTAGCCGACCCACTGTGCAGCGAACGCCTTGGCGTCTTCCTTGCTGCGCACCAGCTTTACGCCGGCCGCtt                                                        >  6:231105/1‑148 (MQ=255)
                                                                                                    tCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTGGCGTCTTCCTTGCTGCGAACCAGCTTTACGCCGCCCGCTTTACcgc                                                  <  1:169042/148‑1 (MQ=255)
                                                                                                      aGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTGGCGTCTTCCTTGCTGCGAACCAGCTTTACGCCGCCCGCTTTACCGCGa                                                <  8:55954/148‑1 (MQ=255)
                                                                                                            ggATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTGGCGTCTTCCTTGCTGCGAACCAGCTTTACGCCGCCCGCTTTACCGCGACCACCcg                                         >  5:1172/1‑149 (MQ=255)
                                                                                                               ttCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTGGCGTCTTCCTTGCTGCGAACCAGCTTTACGCCGCCCGCTTTACCGCGACCACCCGCGt                                      <  2:358349/149‑1 (MQ=255)
                                                                                                               ttCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTGGCGTCTTCCTTGCTGCGAACCAGCTTTACGCCGCCCGCTTTACCGCGACCACCCGCGt                                      <  4:150651/149‑1 (MQ=255)
                                                                                                                   aCCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTGGCGTCTTCCTTGCTGCGAACCAGCTTTACGCCGCCCGCTTTACCGCGACCACCCGCGTGAAc                                  <  4:189695/149‑1 (MQ=255)
                                                                                                                        aaTCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTGGCGTCTTCCTTGCTGCGAACCAGCTTTACGCCGCCCGCTTTACCGCGACCACCCGCGTGAACCTggg                             <  2:127241/149‑1 (MQ=255)
                                                                                                                              ggAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTGGCGTCTTCCTTGCTGCGAACCAGCTTTACGCCGCCCGCTTTACCGCGACCACCCGCGTGAACCTGGGCtttt                        <  6:1172/148‑1 (MQ=255)
                                                                                                                                     accagACCGATGCCAACGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTGGCGTCTTCCTTGCTGCGAACCAGCtt                                                                      <  7:223931/91‑1 (MQ=255)
                                                                                                                                     accagACCGATGCCAACGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTGGCGTCTTCCTTGCTGCGAACCAGCtt                                                                      >  8:223931/5‑95 (MQ=255)
                                                                                                                                      gTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTGGCGTCTTCCTTGCTGCGAACCAGCTTTACGCCGCCCGCTTTACCGCGACCACCCGCGTGAACCTGGGCTTTTACAACCCAc               >  8:311327/1‑149 (MQ=255)
                                                                                                                                         gcccGTTGTCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTGGCGTCTTCCTTGCTGCGAACCAGCTTTACGCCGCCCGCTTTACCGCGACCACCCGCGTGAACCTGGGCTTTTACAACCCACTc             <  4:204981/146‑1 (MQ=255)
                                                                                                                                            cGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTGGCGTCTTCCTTGCTGCGAACCAGCTTTACGCCGCCCGCTTTACCGCGACCACCCGCGTGAACCTGGGCTTTTACAACCCACTCGTTc         >  1:88813/1‑149 (MQ=255)
                                                                                                                                            cGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTGGCGTCTTCCTTGCTGCGAACCAGCTTTACGCCGCCCGCTTTACCGCGACCACCCGCGTGAACCTGGGCTTTTACAACCCACTCGTTc         >  4:133531/1‑149 (MQ=255)
                                                                                                                                              ttGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTGGCGTCTTCCTTGCTGCGAACCAGCTTTACGCCGCCCGCTTTACCGCGACCACCCGCGTGAACCTGGGCTTTTACAACCCACTCGTTccc       >  7:37110/1‑149 (MQ=255)
                                                                                                                                              ttGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTGGCGTCTTCCTTGCTGCGAACCAGCTTTACGCCGCCCGCTTTACCGCGACCACCCGCGTGAACCTGGGCTTTTACAACCCACTCGTTccc       >  1:187083/1‑149 (MQ=255)
                                                                                                                                                   aTCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTGGCGTCTTCCTTGCTGCGAACCAGCTTTACGCCGCCCGCTTTACCGCGACCACCCGCGTGAACCTGGGCtt                          <  7:192728/125‑1 (MQ=255)
                                                                                                                                                   aTCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTGGCGTCTTCCTTGCTGCGAACCAGCTTTACGCCGCCCGCTTTACCGCGACCACCCGCGTGAACCTGGGCtt                          >  8:192728/1‑125 (MQ=255)
                                                                                                                                                   aTCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTGGCGTCTTCCTTGCTGCGAACCAGCTTTACGCCGCCCGCTTTACCGCGACCACCCGCGTGAACCTGGGCTTTTACAACCCACTCGTTCCCGCCGa  >  7:261229/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GACTTTCTCGATGTCCACGCCACCTTCGGTGGAAGCCATGAACACGATACGGCGGCTCGAACGATCGACTACAGCGCCCAGGTACAGCTCTTTGGCGATGTCAGTGCAGGATTCGACCAGAATCTTGGAGACGGGTTGACCGTTGGCATCGGTCTGGTAGGTAACCAGATTCTTGCCCAACCACTGTGCAGCGAACGCCTTGGCGTCTTCCTTGCTGCGAACCAGCTTTACGCCGCCCGCTTTACCGCGACCACCCGCGTGAACCTGGGCTTTTACAACCCACTCGTTCCCGCCGA  >  NC_002947/4729428‑4729723

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑