Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,082,653 C→G 7.2% G538A (GGT→GCT) ‡ alaS ← alanine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,082,6530CG 7.2% 248.8 / 6.9 80G538A (GGT→GCT) ‡alaSalanine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (36/38);  new base G (2/4);  total (38/42)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

GACCAGCGAGCCCTTCTGCTGTACATGCTCGCCCAGTACCTGGCGCAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCAAGAAGCCGGTATCACCCACCTGACCACCGGACTCGGCGTAGAACGGCGTGCGGTCGAGCACCAC  >  NC_002947/5082505‑5082801
                                                                                                                                                    |                                                                                                                                                    
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)
                                                                                                                                                    ccGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCAAGAAGCCGGTATCACCCACCTGACCACCGGACTCGGCGTAGAACGGCGTGCGGTCGAGcaccac  <  6:341062/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GACCAGCGAGCCCTTCTGCTGTACATGCTCGCCCAGTACCTGGCGCAATGCTTCGTGCAGCAGGTGGGTAGCGGAGTGATTCAACGAGGTGGCGTGCTGTACGTCGGCATCGACCTTGGCCTCGACGGGCGAACCGATCACCAGCGCACCGCTGGTCACCACGCCGTGGTGCAGGAAGGCACCGCCGGTCTTGGTGGTGTCGCGCACGTCAAAGCGTGCAGCGCCAGCCTGCAAGAAGCCGGTATCACCCACCTGACCACCGGACTCGGCGTAGAACGGCGTGCGGTCGAGCACCAC  >  NC_002947/5082505‑5082801

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: