Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,082,651 C→G 5.8% A539P (GCG→CCG)  alaS ← alanine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,082,6510CG 5.8% 262.2 / 4.0 81A539P (GCG→CCG) alaSalanine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (38/38);  new base G (2/3);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCAGCGAGCCCTTCTGCTGTACATGCTCGCCCAGTACCTGGCGCAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCAAGAAGCCGGTATCACCCACCTGACCACCGGACTCGGCGTAGAACGGCGTGCGGTCGAGCAC  >  NC_002947/5082505‑5082798
                                                                                                                                                  |                                                                                                                                                   
gACCAGCGAGCCCTTCTGCTGTACATGCTCGCCCAGTACCTGGCGCAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCAc                                                                                                                                                   >  1:42267/1‑149 (MQ=255)
gACCAGCGAGCCCTTCTGCTGTACATGCTCGCCCAGTACCTGGCGCAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCAc                                                                                                                                                   >  6:219396/1‑149 (MQ=255)
     gCGAGCCCTTCTGCTGTACATGCTCGCCCAGTACCTGGCGCAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTg                                                                                                                                              <  4:146146/149‑1 (MQ=255)
        aGCCCTTCTGCTGTACATGCTCGCCCAGTACCTGGCGCAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGt                                                                                                                                            >  1:228404/1‑148 (MQ=255)
         gCCCTTCTGCTGTACATGCTCGCCCAGTACCTGGCGCAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGt                                                                                                                                            <  5:158704/147‑1 (MQ=255)
           ccTTCTGCTGTACATGCTCGCCCAGTACCTGGCGCAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTcacc                                                                                                                                        <  4:257421/149‑1 (MQ=255)
           ccTTCTGCTGTACATGCTCGCCCAGTACCTGGCGCAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTcacc                                                                                                                                        >  6:125565/1‑149 (MQ=255)
           ccTTCTGCTGTACATGCTCGCCCAGTACCTGGCGCAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCGGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGGTGCGCTGGtgatc                                                                                                                                        <  1:356425/149‑5 (MQ=255)
             ttCTGCTGTACATGCTCGCCCAGTACCTGGCGCAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTcaccac                                                                                                                                      >  3:7168/1‑149 (MQ=255)
              tCTGCTGTACATGCTCGCCCAGTACCTGGCGCAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACg                                                                                                                                     >  2:123110/1‑149 (MQ=255)
               ctgctgTACATGCTCGCCCAGTACCTGGCGCAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGGTGCGCTGGtgatcgg                                                                                                                                      <  4:324560/147‑7 (MQ=255)
               ctgctgTACATGCTCGCCCAGTACCTGGCGCAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTcaccac                                                                                                                                      <  5:120589/147‑1 (MQ=255)
               ctgctgTACATGCTCGCCCAGTACCTGGCGCAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACAGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGc                                                                                                                                    >  6:258965/1‑149 (MQ=255)
                  ctgTACATGCTCGCCCAGTACCTGGCGCAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGc                                                                                                                                    <  5:125565/146‑1 (MQ=255)
                  ctgTACATGCTCGCCCAGTACCTGGCGCAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCgt                                                                                                                                 <  6:283645/149‑1 (MQ=255)
                   tgTACATGCTCGCCCAGTACCTGGCGCAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCgtg                                                                                                                                <  1:123110/149‑1 (MQ=255)
                    gTACATGCTCGCCCAGTACCTGGCGCAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCgtgg                                                                                                                               >  8:202799/1‑149 (MQ=255)
                         tGCTCGCCCAGTACCTGGCGCAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAg                                                                                                                          <  1:188393/149‑1 (MQ=255)
                          gCTCGCCCAGTACCTGGCGCAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGGTGCGCTGGtgatc                                                                                                                                        >  7:151799/1‑130 (MQ=255)
                          gCTCGCCCAGTACCTGGCGCAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGGTGCGCTGGtgatc                                                                                                                                        <  8:151799/134‑5 (MQ=255)
                          gCTCGCCCAGTACCTGGCGCAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAg                                                                                                                          >  5:341062/1‑148 (MQ=255)
                          gCTCGCCCAGTACCTGGCGCAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAg                                                                                                                          >  2:308355/1‑148 (MQ=255)
                            tCGCCCAGTACCTGGCGCAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAgg                                                                                                                         >  8:61583/1‑147 (MQ=255)
                              gCCCAGTACCTGGCGCAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGCAAgg                                                                                                                     >  8:270084/1‑149 (MQ=255)
                                      ccTGGCGCAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCAccgccg                                                                                                             >  5:370796/1‑149 (MQ=255)
                                       gtGGCGCAATGCTTCGTGCAGCTGGTGGGTAGCGGAGTGATTCAACGTGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCTGCGCAGCGCTGGTCACCACGCGGTGGTGCAGGAAGGCACCGCCgg                                                                                                            <  8:293559/148‑1 (MQ=255)
                                          gcgcAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGt                                                                                                           <  2:228404/147‑1 (MQ=255)
                                          gcgcAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTc                                                                                                          <  4:274607/148‑1 (MQ=255)
                                          gcgcAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTc                                                                                                          <  6:238071/148‑1 (MQ=255)
                                            gcAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGAGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGAGGTGCAGGAAGGCACCGCCGGTCTtg                                                                                                       >  6:244225/1‑149 (MQ=255)
                                            gaaaTGCTTCGTGCAGCAGGTGGGTAGCTGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGAACTTGGCCGCGACCGGCGAACCGATCACAAGCGCACCGCTGGTCACCACGCCGAGGTGCAGGAAGGCACCGCAGCTCTtg                                                                                                       >  6:189416/3‑149 (MQ=255)
                                             cAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTtgg                                                                                                      >  2:341237/1‑149 (MQ=255)
                                               aTGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTtggtg                                                                                                    <  4:111029/149‑1 (MQ=255)
                                               aTGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTtggtg                                                                                                    <  4:7168/149‑1 (MQ=255)
                                                tGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGc                                                                                                                                                <  1:278724/104‑1 (MQ=255)
                                                tGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGc                                                                                                                                                >  2:278724/1‑104 (MQ=255)
                                                         cagcagGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCa                                                                                          <  8:219795/149‑1 (MQ=255)
                                                         cagcagGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCa                                                                                          >  2:298923/1‑149 (MQ=255)
                                                          agcagGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCAc                                                                                         <  7:189562/149‑1 (MQ=255)
                                                               gTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTtg                                                                                                       <  8:132612/130‑1 (MQ=255)
                                                               gTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTtg                                                                                                       >  7:132612/1‑130 (MQ=255)
                                                                tGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCaaa                                                                                   <  5:329037/149‑1 (MQ=255)
                                                                        ggAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAg                                                                           >  6:320755/1‑149 (MQ=255)
                                                                         gAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAgc                                                                          <  3:319665/149‑1 (MQ=255)
                                                                              aTTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAg                                                                     >  4:355593/1‑149 (MQ=255)
                                                                                 cAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCt                                                                  >  7:347328/1‑149 (MQ=255)
                                                                                  aaCGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTg                                                                 <  7:183961/149‑1 (MQ=255)
                                                                                   aCGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTg                                                                 >  6:318998/1‑148 (MQ=255)
                                                                                   aCGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGc                                                                <  5:219396/149‑1 (MQ=255)
                                                                                    cGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGc                                                                >  1:402/1‑148 (MQ=255)
                                                                                    cGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCa                                                               <  5:189416/149‑1 (MQ=255)
                                                                                     gAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCaa                                                              <  5:308057/149‑1 (MQ=255)
                                                                                        gTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCaagaa                                                           <  1:341237/149‑1 (MQ=255)
                                                                                             gTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCAAGAAGCCgg                                                      >  6:178673/1‑149 (MQ=255)
                                                                                              tGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCAAGAAGCCgg                                                      >  6:51178/1‑148 (MQ=255)
                                                                                              tGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCAAGAAGCCGGt                                                     <  7:202799/149‑1 (MQ=255)
                                                                                                cTGTACGTCGGCATCGACCTTTGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCAAGAAGCCGGTAt                                                   >  2:307902/1‑149 (MQ=255)
                                                                                                  gTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTTGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCAAGAAGCCGGTAtac                                                 >  5:87015/1‑147 (MQ=255)
                                                                                                   tACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCAAGAAGCCGGTATcac                                                <  6:69592/149‑1 (MQ=255)
                                                                                                   tACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCAAGAAGCCGGTATca                                                 <  8:310431/148‑1 (MQ=255)
                                                                                                    aCGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCCCCGCCGGTCTTGGTGGTGTCGCGCCCGTCCAAGCGTGCCGCGCCCGCCTGCAAGAAGCCGGTATcacc                                               >  6:264096/1‑149 (MQ=255)
                                                                                                         ggCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGGGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCAAGAAGCCGGTATCACCCACCt                                          <  6:370796/149‑1 (MQ=255)
                                                                                                             tCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTACAAGAAGCCGGTATCACCCACCTGac                                       >  3:307728/1‑148 (MQ=255)
                                                                                                              cGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCAAGAAGCCGGTATCACCCACCTGacca                                     >  7:345680/1‑149 (MQ=255)
                                                                                                               gACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCAAGAAGCCGGTATCACCCACCTGaccac                                    <  5:58944/149‑1 (MQ=255)
                                                                                                                 ccTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTTGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCAAGAAGCCGGTATCACCCACCTGACCACCg                                  <  8:345680/149‑1 (MQ=255)
                                                                                                                  cTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCATGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCATGCAAGAAGCAGGTATCCCCCACCTGACCACCg                                  >  4:208510/1‑148 (MQ=255)
                                                                                                                     ggCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTACAAGAAGCCGGTATCACCCACCTGACCACCGGAc                               <  4:307728/148‑1 (MQ=255)
                                                                                                                      gCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCAAGAAGCCGGTATCACCCACCTGACCACCGGACTc                             >  2:178052/1‑149 (MQ=255)
                                                                                                                        cTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCAAGAAGCCGGTATCACCCACCTGACCACCGGACTCgg                           >  8:308307/1‑149 (MQ=255)
                                                                                                                         tCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCAAGAAGCCGGTATCACCCACCTGACCACCGGACTCGGc                          >  8:164598/1‑149 (MQ=255)
                                                                                                                           gACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCAAGAAGCCGGTATCACCCACCTGACCACCGGACTCGGc                          <  5:318998/147‑1 (MQ=255)
                                                                                                                                gCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCAAGAAGCCGGTATCACCCACCTGACCACCGGACTCGGCGTAGAAc                   >  4:241757/1‑149 (MQ=255)
                                                                                                                                   aaCCGATCACCAGCGCACCGCTGGTCACCACGGCGTTGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCAAGAATCCGGTATCACCCACCTGACCACCGGACTCGGCGTAGAACGGc                <  5:244225/149‑1 (MQ=255)
                                                                                                                                     ccGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCAAGAAGCCGGTATCACCCACCTGACCACCGGACTCGGCGTAGAACGGCGt              <  7:308307/149‑1 (MQ=255)
                                                                                                                                       gATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCAAGAAGCCGGTATCACCCACCTGACCACCGGACTCGGCGTAGAACGGCGTGc            <  3:355593/149‑1 (MQ=255)
                                                                                                                                        gtgaCCAGCGGTGCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCTAAGCGTGCAGCGCCAGCCTGCAAGAAGCCGGTATCACCCACCTGACCACCGGACTCGGCGTAGAACGGCGTGCg           >  4:206383/4‑149 (MQ=255)
                                                                                                                                        aTCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCAAGAAGCCGGTATCACCCACCTGACCACCGGACTCGGCGTAGAACGGCGTGCg           >  1:332729/1‑149 (MQ=255)
                                                                                                                                           aCCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCAAGAAGCCGGTATCACCCACCTGACCACCGGACTCGGCGTAGAACGGCGTGCGGt         <  5:320755/148‑1 (MQ=255)
                                                                                                                                            ccAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCAAGAAGCCGGTATCACCCACCTGACCACCGGACTCGGCGTAGAACGGCGTGCGGTc        >  7:110163/1‑148 (MQ=255)
                                                                                                                                              aGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCAAGAAGCCGGTATCACCCACCTGACCACCGGACTCGGCGTAGAACGGCGTGCGGTCGAg     <  2:402/149‑1 (MQ=255)
                                                                                                                                                 gcACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCAAGAAGCCGGTATCACCCACCTGACCACCGGACTCGGCGTAGAACGGCGTGCGGTCGAGcac  <  5:51178/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GACCAGCGAGCCCTTCTGCTGTACATGCTCGCCCAGTACCTGGCGCAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCAAGAAGCCGGTATCACCCACCTGACCACCGGACTCGGCGTAGAACGGCGTGCGGTCGAGCAC  >  NC_002947/5082505‑5082798

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