Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,553,989 T→A 8.3% intergenic (+117/+64) pydA → / ← rutR NADP‑dependent dihydropyrimidine dehydrogenase subunit PreA/transcriptional regulator, TetR family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,553,9890TA 8.3% 262.4 / 11.7 84intergenic (+117/+64)pydA/rutRNADP‑dependent dihydropyrimidine dehydrogenase subunit PreA/transcriptional regulator, TetR family
Reads supporting (aligned to +/- strand):  ref base T (49/28);  new base A (5/2);  total (54/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 = 8.99e-01

GATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTCTGCACCGCATTCTCGAAGGCCATGTCAGACAGCAC  >  NC_002947/4553840‑4554133
                                                                                                                                                     |                                                                                                                                                
gATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCtt                                                                                                                                                  <  1:3625085/150‑1 (MQ=255)
gATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCtt                                                                                                                                                  <  1:2096255/150‑1 (MQ=255)
   gcgtgatcttgcgggggatcgggctggagccgtgttccggcattccaggtggacgtggtcttttTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTcgc                                                                                                                                               >  1:2504935‑M1/65‑150 (MQ=255)
    cgtgatcttgcgggggatcgggctggagccgtgatccggcattccaggtggacgtggtcttttTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTcgcg                                                                                                                                              >  1:1109144‑M1/64‑150 (MQ=255)
    cgcgTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcg                                                                                                                                              >  2:3391248‑M2/1‑147 (MQ=255)
    cgcgTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcg                                                                                                                                              >  2:346354‑M2/1‑147 (MQ=255)
     gcgTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcgg                                                                                                                                             >  1:767965‑M2/1‑146 (MQ=255)
     gcgTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcgg                                                                                                                                             >  2:4986638‑M2/1‑146 (MQ=255)
       gTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcc                                                                                                                                           >  2:2094882‑M2/1‑144 (MQ=255)
       gTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcc                                                                                                                                           >  1:5236537‑M2/1‑144 (MQ=255)
          cttgcgggggatcgggctggagccgtgatccggcattccaggtggacgtggtcttttTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGa                                                                                                                                        >  1:3222274‑M1/58‑150 (MQ=255)
           ttgcgggggatcgggctggagccgtgatccggcattccaggtggacgtggtcttttTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCGAAGcggcccctc                                                                                                                                       >  1:3986652‑M1/57‑147 (MQ=255)
           ccccTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCGAAgcggcccctc                                                                                                                                       >  1:2502783‑M2/1‑140 (MQ=255)
           ccccTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCGAAgcggcccctc                                                                                                                                       >  1:2502787‑M2/1‑140 (MQ=255)
            cccTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatc                                                                                                                                      >  2:4858568‑M2/1‑139 (MQ=255)
            cccTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatc                                                                                                                                      >  1:720746‑M2/1‑139 (MQ=255)
                aCCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatcgcga                                                                                                                                  >  1:3537332‑M2/1‑135 (MQ=255)
                                  gCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcg                                                                                                                                              <  2:4035707‑M2/120‑4 (MQ=255)
                                  gCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcg                                                                                                                                              >  1:4035707‑M2/1‑117 (MQ=255)
                                            gCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatcgcga                                                                                                                                  >  1:4492935‑M2/1‑107 (MQ=255)
                                            gCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACAACGAGGGGCCGCTTCgcggcccgatcgcga                                                                                                                                  <  2:4492935‑M2/122‑16 (MQ=255)
                                                      cTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggccc                                                                                                                                          >  2:66952‑M2/1‑97 (MQ=255)
                                                      cTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggccc                                                                                                                                          <  1:66952‑M2/104‑8 (MQ=255)
                                                                     tgGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCGAAgcggcccctc                                                                                                                                       <  2:3637275‑M2/92‑11 (MQ=255)
                                                                     tgGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCGAAgcggcccctc                                                                                                                                       >  1:3637275‑M2/1‑82 (MQ=255)
                                                                                     TCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATcccc                                                             <  1:4354438‑M1/150‑1 (MQ=255)
                                                                                     TCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATcccc                                                             <  2:1295271‑M1/150‑1 (MQ=255)
                                                                                      actgctgGGGCCGCTGTGCGGCCCGAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATccccc                                                            <  1:1252848‑M1/143‑1 (MQ=255)
                                                                                        tgctgGGGCCGCTGTGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGc                                                          <  1:2791670‑M1/145‑1 (MQ=255)
                                                                                        tgcTCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGc                                                          <  1:1135857‑M1/147‑1 (MQ=255)
                                                                                        tgcTCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGc                                                          <  1:262103‑M1/147‑1 (MQ=255)
                                                                                         gATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatcgcgacacaaggccgctcccacaggggccgcgtcagctgttcaggcggtgacccgcttcaactggcctcacggtaggg                                                         <  1:4658833‑M2/150‑89 (MQ=255)
                                                                                         gATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatcgcgacacaaggccgctcccacaggggccgcgtcagctgttcaggcggtgacccgcttcaactggcctcacggtaggg                                                         <  2:721187‑M2/150‑89 (MQ=255)
                                                                                         gATCGGGCCGCACAGCGGCCCCAGCAGTCGTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatcgcgacacaaggccgctcccacaggggccgcgtcagctgttcaggcggtgacccgcttcaactggcctcacgggcggg                                                         <  2:3648997‑M2/150‑89 (MQ=255)
                                                                                         GATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCa                                                         <  1:5216457‑M1/150‑1 (MQ=255)
                                                                                          aTCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatcgcgacacaaggccgctcccacaggggccgcgtcagctgttcaggcggtgacccgcttcaactggcctcacggtagggg                                                        <  2:2414742‑M2/150‑90 (MQ=255)
                                                                                             GGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGAt                                                     <  2:695289‑M1/150‑1 (MQ=255)
                                                                                              GGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATc                                                    <  1:1796327‑M1/150‑1 (MQ=255)
                                                                                              GGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATc                                                    <  1:1163057‑M1/150‑1 (MQ=255)
                                                                                               GCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCa                                                   <  2:53784‑M1/150‑1 (MQ=255)
                                                                                                CCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCAc                                                  <  1:3927719‑M1/150‑1 (MQ=255)
                                                                                                 cGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatcgcgacacaaggccgctcccgcaggggccgcgtcagctgttcaggcggtgacccgcttcaactggcctcacggtaggggttacgcg                                                 <  2:1803120‑M2/150‑97 (MQ=255)
                                                                                                 cGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatcgcgacacaaggccgctcccacaggggccgcgtcagctgttcaggcggtgacccgcttcaactggcctcacggtaggggttacgcg                                                 <  2:1825420‑M2/150‑97 (MQ=255)
                                                                                                  gCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatcgcgacacaaggccgctcccacaggggccgcgtcagctgttcaggcggtgacccgcttcaactggcctcacggtaggggttacgcgg                                                <  2:60657‑M2/150‑98 (MQ=255)
                                                                                                   cacaGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatcgcgacacaaggccgctcccacaggggccgcgtcagctgttcaggcggtgacccgcttcaactggc                                                                     >  2:3135414‑M2/1‑52 (MQ=255)
                                                                                                   cacaGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatcgcgacacaaggccgctcccacaggggccgcgtcagctgttcaggcggtgacccgcttcaactggc                                                                     <  1:3135414‑M2/128‑77 (MQ=255)
                                                                                                   CACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCt                                               <  1:877172‑M1/150‑1 (MQ=255)
                                                                                                   CACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCt                                               <  1:877598‑M1/150‑1 (MQ=255)
                                                                                                     caGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCGAAgcggcccctc                                                                                                                                       >  2:1789645‑M2/1‑50 (MQ=255)
                                                                                                     caGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCGAAgcggcccctc                                                                                                                                       <  1:1789645‑M2/60‑11 (MQ=255)
                                                                                                         ggCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACaaa                                                                                                              >  2:1847016/1‑81 (MQ=255)
                                                                                                         ggCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACaaa                                                                                                              <  1:1847016/81‑1 (MQ=255)
                                                                                                            cccAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCt                                               <  2:2332956/141‑1 (MQ=255)
                                                                                                            cccAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCt                                               <  2:2333029/141‑1 (MQ=255)
                                                                                                            cccAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCt                                               >  1:2333029/1‑141 (MQ=255)
                                                                                                            cccAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCt                                               >  1:2332956/1‑141 (MQ=255)
                                                                                                            cccAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGt                                      >  2:4492041/1‑150 (MQ=255)
                                                                                                             ccAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTc                                     >  1:1403971/1‑150 (MQ=255)
                                                                                                             ccAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTc                                     >  1:3296899/1‑150 (MQ=255)
                                                                                                                 cagTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTCTGCa                                 >  2:4282785/1‑150 (MQ=255)
                                                                                                                   gTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatcgcgacacaaggccgctcccacaggggccgcgtcagctgttcaggcggtgacccgcttcaactggcctcacggtaggggttacgcggatcctgcgtccagttca                               >  1:60657‑M2/1‑36 (MQ=255)
                                                                                                                     cTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTCTGCACCGc                             >  2:2802744/1‑150 (MQ=255)
                                                                                                                     cTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTCTGCACCGc                             >  1:850253/1‑150 (MQ=255)
                                                                                                                      tttGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTCTGCACCGcc                            >  2:4286160/1‑149 (MQ=255)
                                                                                                                      tttGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTCTGCACCGCa                            >  2:5158597/1‑150 (MQ=255)
                                                                                                                      tttGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTCTGCACCGCa                            >  2:3661107/1‑150 (MQ=255)
                                                                                                                      tttGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTCTGCACCGCa                            >  2:4286031/1‑150 (MQ=255)
                                                                                                                      tttGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTCTGCACCGCa                            >  2:86723/1‑150 (MQ=255)
                                                                                                                         gCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTCTGCACCGCATtc                         >  1:2611251/1‑150 (MQ=255)
                                                                                                                           cgcgATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTCTGCACCGCATtctc                       >  2:3808706/1‑150 (MQ=255)
                                                                                                                            gcgATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTCTGCACCGCATTCTCg                      >  1:4452194/1‑150 (MQ=255)
                                                                                                                                   aGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTCTGCATCGCATTCTCGAAGGCCa               >  2:1271757/1‑150 (MQ=255)
                                                                                                                                   aGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTCTGCACCGCATTCTCGAAGGCCa               >  1:527279/1‑150 (MQ=255)
                                                                                                                                       aTCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTCTGCACCGCATTCTCGAAGGCCATGTc           >  2:4127801/1‑150 (MQ=255)
                                                                                                                                       aTCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTCTGCACCGCATTCTCGAAGGCCATGTc           >  1:2124304/1‑150 (MQ=255)
                                                                                                                                        tCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTCTGCACCGCATTCTCGAAGGCCATGTCa          >  2:3827650/1‑150 (MQ=255)
                                                                                                                                         cGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTCTGCACCGCATTCTCGAAGGCCATGTCAg         >  1:4225505/1‑150 (MQ=255)
                                                                                                                                         cGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTCTGCACCGCATTCTCGAAGGCCATGTCAg         >  1:4225618/1‑150 (MQ=255)
                                                                                                                                          gAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTCTGCACCGCATTCTCGAAGGCCATGTCAGa        >  2:5029912/1‑150 (MQ=255)
                                                                                                                                          gAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTCTGCACCGCATTCTCGAAGGCCATGTCAGa        >  1:1861216/1‑150 (MQ=255)
                                                                                                                                          gAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTCTGCACCGCATTCTCGAAGGCCATGTCAGa        >  2:1186674/1‑150 (MQ=255)
                                                                                                                                            ggggCCGCTTCgcggcccgatcgcgacacaaggccgctcccacaggggccgcgtcagctgttcaggcggtgacccgcttcaactggcctcacggtaggggttacgcggatcctgcgtccagttcaggtacggcttgccggtatcctgtgc      >  1:721187‑M2/1‑11 (MQ=255)
                                                                                                                                              ggCCGCTTCgcggcccgatcgcgacacaaggccgctcccacaggggccgcgtcagctgttcaggcggtgacccgcttcaactggcctcacggtaggggttacgcggatcctgcgtccagttcaggtacggcttgccggtatcctgtgcct    >  2:4658833‑M2/1‑9 (MQ=255)
                                                                                                                                               gCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTCTGCACCGCATTCTCGAAGGCCATGTCAGACAGCa   >  2:3245774/1‑150 (MQ=255)
                                                                                                                                               gCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTCTGCACCGCATTCTCGAAGGCCATGTCAGACAGCa   >  2:1777336/1‑150 (MQ=255)
                                                                                                                                                ccGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTCTGCACCGCATTCTCGAAGGCCATGTCAGACAGCAc  >  2:2537772/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTCTGCACCGCATTCTCGAAGGCCATGTCAGACAGCAC  >  NC_002947/4553840‑4554133

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: