Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,086,998 C→A 6.5% A109S (GCC→TCC)  astB ← N‑succinylarginine dihydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,086,9980CA 6.5% 197.4 / 3.8 62A109S (GCC→TCC) astBN‑succinylarginine dihydrolase
Reads supporting (aligned to +/- strand):  ref base C (30/28);  new base A (2/2);  total (32/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 = 4.11e-01

TATTTCTCGTCATTGAACATGGCGCCCAGCACGCGGCTGGTGGTCGGGTGCTCGATGCTGCGGTGATATTTGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCC  >  NC_002947/5086850‑5087146
                                                                                                                                                    |                                                                                                                                                    
tATTTCTCGTCATTGAACATGGCGCCCAGCACGCGGCTGGTGGTCGGGTGCTCGATGCTGCGGTGATATTTGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGc                                                                                                                                                      <  2:241765/149‑1 (MQ=255)
         tCATTGAACATGGCGCCCAGCACGCGGCTGGTGGTCGGGTGCTCGATGCTGCGGTGATATTTGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACAt                                                                                                                                             <  3:353803/149‑1 (MQ=255)
               aaCATGGCGCCCAGCACGCGGCTGGTGGTCGGGTGCTCGATGCTGCGGTGATATTTGCAGTTGAGGTTGGCGGCGGTGAAGTGCAGGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGa                                                                                                                                       <  3:291916/149‑1 (MQ=255)
                aCATGGCGCCCAGCACGCGGCTGGTGGTCGGGTGCTCGATGCTGCGGTGTTTATTGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGa                                                                                                                                       <  1:338169/148‑1 (MQ=255)
                 cATGGCGCCCAGCACGCGGCTGGTGGTCGGGTGCTCGATGCTGCGGTGATATTTGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAgg                                                                                                                                     <  4:256470/149‑1 (MQ=255)
                  aTGGCGCCCAGCACGCGGCTGGTGGTCGGGTGCTCGATGCTGCGGTGATATTTGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGc                                                                                                                                    >  2:53145/1‑149 (MQ=255)
                          cAGCACGCGGCTGGTGGTCGGGTGCTCGATGCTGCGGTGATATTTGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAAc                                                                                                                            >  6:285654/1‑149 (MQ=255)
                            gCACGCGGCTGGTGGTCGGGTGCTCGATGCTGCGGTGATATTTGCAGTTGAGGTTGGCGGTGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTg                                                                                                                          <  3:119916/149‑1 (MQ=255)
                             cACGCGGCTGGTGGTCGGGTGATCGATGCTGCGGTGATATTTGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTgg                                                                                                                         >  5:300635/1‑149 (MQ=255)
                                gcgGCTGGTGGTCGGGTGCTCGATGCTGCGGTGATATTTGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTg                                                                                                                          >  2:220784/1‑145 (MQ=255)
                                gcgGCTGGTGGTCGGGTGCTCGATGCTGCGGTGATATTTGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTg                                                                                                                          <  1:220784/145‑1 (MQ=255)
                                        tggtCGGGTGCTCGATGCTGCGGTGATATTTGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAAc                                                                                                              >  6:219415/1‑149 (MQ=255)
                                             gggTGCTCGATGCTGCGGTGATATTTGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCAt                                                                                                         <  7:116136/149‑1 (MQ=255)
                                                  cTCGATGCTGCGGTGATATTTGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGc                                                                                                                                    <  4:248281/117‑1 (MQ=255)
                                                  cTCGATGCTGCGGTGATATTTGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGc                                                                                                                                    >  3:248281/1‑117 (MQ=255)
                                                  cTCGATGCTGCGGTGATATTTGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCt                                                                                                    >  4:384564/1‑149 (MQ=255)
                                                   tCGATGCTGCGGTGATATTTGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATgg                                                                                                       <  4:206210/145‑1 (MQ=255)
                                                   tCGATGCTGCGGTGATATTTGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATgg                                                                                                       >  3:206210/1‑145 (MQ=255)
                                                   tCGATGCTGCGGTGATATTTGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTc                                                                                                   <  6:244990/149‑1 (MQ=255)
                                                              gTGATATTTGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGACACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGgcgc                                                                                        >  2:30485/1‑149 (MQ=255)
                                                              gTGATATTTGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGACACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGgcgc                                                                                        >  2:36602/1‑149 (MQ=255)
                                                                gATATTTGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGgcgcgc                                                                                      <  5:285654/149‑1 (MQ=255)
                                                                  tatTTGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGACCTCGCGGGCGGCGCGCTg                                                                                    <  6:300635/149‑1 (MQ=255)
                                                                   atTTGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGccaccca                                                                                                                                                <  4:169855/88‑1 (MQ=255)
                                                                   atTTGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGccaccca                                                                                                                                                >  3:169855/1‑88 (MQ=255)
                                                                    tttGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTTACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCt                                                                                     <  8:174003/146‑1 (MQ=255)
                                                                    tttGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTTACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCt                                                                                     >  7:174003/1‑146 (MQ=255)
                                                                    tttGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGACACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGa                                                                                  <  1:36602/149‑1 (MQ=255)
                                                                    tttGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGACACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGa                                                                                  <  1:30485/149‑1 (MQ=255)
                                                                         aGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACt                                                                             >  3:350972/1‑149 (MQ=255)
                                                                             gAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCAt                                                                                                         >  2:375988/1‑117 (MQ=255)
                                                                             gAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCAt                                                                                                         <  1:375988/117‑1 (MQ=255)
                                                                              aGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGGCTCGCGGGCGGCGCGCTGGATAACTTCgg                                                                         >  5:258073/1‑148 (MQ=255)
                                                                               ggTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCa                                                                       <  7:147101/149‑1 (MQ=255)
                                                                                gTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCAt                                                                      >  7:196695/1‑149 (MQ=255)
                                                                                     cggcggTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAgg                                                                                                                                     >  1:171850/1‑81 (MQ=255)
                                                                                     cggcggTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAgg                                                                                                                                     <  2:171850/81‑1 (MQ=255)
                                                                                      ggcggTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGc                                                                <  5:173853/149‑1 (MQ=255)
                                                                                           tGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACt                                                            >  1:402891/1‑148 (MQ=255)
                                                                                               gTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAAcc                                                       >  5:123404/1‑149 (MQ=255)
                                                                                                     gcgACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAggc                                                 <  1:53145/149‑1 (MQ=255)
                                                                                                     gcgACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGCATAACTTCGGCATCGCTGCCACTGAACCCCAggc                                                 >  7:4597/1‑149 (MQ=255)
                                                                                                           gTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCCGgcggcg                                             >  4:17211/1‑147 (MQ=255)
                                                                                                            tCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAg                                          >  8:236669/1‑149 (MQ=255)
                                                                                                             cagcagTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAgc                                         >  6:225342/1‑149 (MQ=255)
                                                                                                              agcagTGTCGGCACTCGGGCTGACTGGTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAgcg                                        >  1:220867/1‑149 (MQ=255)
                                                                                                               gcagTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAgcgc                                       >  7:53056/1‑149 (MQ=255)
                                                                                                                   tgtCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCgccgcc                                   <  6:123404/149‑1 (MQ=255)
                                                                                                                         cACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTc                                                                           >  4:2464/1‑103 (MQ=255)
                                                                                                                         cACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTc                                                                           <  3:2464/103‑1 (MQ=255)
                                                                                                                          aCTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCa                                                              <  1:246473/115‑1 (MQ=255)
                                                                                                                          aCTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCa                                                              >  2:246473/1‑115 (MQ=255)
                                                                                                                              gggCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTTAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGgc                         <  3:384564/148‑1 (MQ=255)
                                                                                                                              gggCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCa                                                    >  7:21479/1‑121 (MQ=255)
                                                                                                                              gggCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCa                                                    <  8:21479/121‑1 (MQ=255)
                                                                                                                              gggCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGgcg                        >  6:2624/1‑149 (MQ=255)
                                                                                                                               ggCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGgcgc                       <  5:219415/149‑1 (MQ=255)
                                                                                                                                  tGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGctct                    >  1:252599/1‑149 (MQ=255)
                                                                                                                                      tGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTtgtg                >  8:73850/1‑149 (MQ=255)
                                                                                                                                            ggcgTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTtgt                 >  6:247712/1‑142 (MQ=255)
                                                                                                                                            ggcgTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTtgt                 <  5:247712/142‑1 (MQ=255)
                                                                                                                                             gcgTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAg         <  8:196695/149‑1 (MQ=255)
                                                                                                                                                ttGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAac      <  5:225342/149‑1 (MQ=255)
                                                                                                                                                  ggCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAaca     <  5:2624/148‑1 (MQ=255)
                                                                                                                                                  ggCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCACAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAacac    >  6:294374/1‑149 (MQ=255)
                                                                                                                                                    ccacccacATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACAccc  >  1:142437/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TATTTCTCGTCATTGAACATGGCGCCCAGCACGCGGCTGGTGGTCGGGTGCTCGATGCTGCGGTGATATTTGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCC  >  NC_002947/5086850‑5087146

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