Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 4,022,806 C→G 14.0% intergenic (‑501/+167) PP_3547 ← / ← emrB putative NAD(P)‑binding oxidoreductase with NAD(P)‑binding Rossmann‑fold domain/multidrug efflux transport system ‑ permease subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154514,022,8060CG14.0% 129.9 / 14.7 43intergenic (‑501/+167)PP_3547/emrBputative NAD(P)‑binding oxidoreductase with NAD(P)‑binding Rossmann‑fold domain/multidrug efflux transport system ‑ permease subunit
Reads supporting (aligned to +/- strand):  ref base C (24/13);  new base G (3/3);  total (27/16)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.94e-01

TCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCGGGATCTTGCCTTGATTTTAGCTCAATGCCCCGCCGATGCCGGGCCCGCCTTAGCGGTAAAGGGCGGCTTGGCGAACCACACCAGCAG  >  AE015451/4022605‑4023038
                                                                                                                                                                                                          |                                                                                                                                                                                                                                        
tcattcaCGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCt                                                                                                                                                                                                        >  1:406279/1‑236 (MQ=255)
     caCGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCg                                                                                                                                                                                                                     >  1:492627/1‑218 (MQ=255)
          aGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCg                                                                                                                                                                                                                     >  1:838619/1‑213 (MQ=255)
                  agcgcgacgcggtccaagcacggctttcaaacccagcgctcttgatctcgatcttgatcttgatctcgatcttgatcttgatcttgatcttgatcttgatcttgatcttgatcttgatcttgatcttgatcttcgcgacttcaggaggccgaacgcaggccttgcgaaggcaggtgacgggcatCCATGCCCGTCACCTCCCTTCGCAAGCCCTACGCTCGGCCTCCTGAAGTCg                                                                                                                                                                                        >  1:183839‑M2/185‑235 (MQ=255)
                                    cTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCtt                                                                                                                                                                                                                       >  2:717715/1‑185 (MQ=255)
                                    cTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCtt                                                                                                                                                                                                                       <  1:717718/185‑1 (MQ=255)
                                     ttACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCt                                                                                                                                                                    >  1:636628/1‑235 (MQ=255)
                                         ggcggGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAg                                                                                                                                                                                                                 <  1:122869/186‑1 (MQ=255)
                                           ttcaaacccagcgctcttgatctcgatcttgatcttgatctcgatcttgatcttgatcttgatcttgatcttgatcttgatcttgatcttgatcttgatcttgatcttcgcgacttcaggaggccgaacgcaggccttgcgaaggcaggtgacgggcatGGATGCCCGTCACCTGCCTTCGCaa                                                                                                                                                                                                                  >  2:645188‑M2/160‑184 (MQ=255)
                                           ttcaaacccagcgctcttgatctcgatcttgatcttgatctcgatcttgatcttgatcttgatcttgatcttgatcttgatcttgatcttgatcttgatcttgatcttcgcgacttcaggaggccgaacgcaggccttgcgaaggcaggtgacgggcatGGATGCCCGTCACCTGCCTTCGCaa                                                                                                                                                                                                                  <  1:645191‑M2/25‑1 (MQ=255)
                                             ggTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCtt                                                                                                                                                                                                                       <  2:651554/176‑1 (MQ=255)
                                                                            cttgatctcgatcttgatcttgatcttgatcttgatcttgatcttgatcttgatcttgatcttgatcttgatcttcgcgacttcaggaggccgaacgcaggccttgcgaaggcaggtgacgggcatCCATGCCCGTCACCTCCCTTCGCAAGCCCTACGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTATCGCGCGCTTAGAAgcaaca                                                                                                                                               >  1:197625‑M2/127‑216 (MQ=255)
                                                                                    aCGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGcgaa                                                                                                                                                                                    >  2:451763/1‑170 (MQ=11)
                                                                                    aCGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGaa                                                                                                                                                                                            >  2:634312/1‑164 (MQ=40)
                                                                                           cTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGaa                                                                                                                                   >  2:387240/1‑214 (MQ=255)
                                                                                              ccatcatcCGGCCCCTGCGCTGCGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagagcgagagcgaga                                                                                                           <  1:649794/233‑1 (MQ=255)
                                                                                              ccatcatcCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCaa                                                                                                                                                                                                                  <  1:339956/130‑1 (MQ=18)
                                                                                                                 cTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCCAGCCCTGCGCTCGGCCTCCTGAAGGCGCAATCCGCGGCGCCTg                                                                                                                                                                        >  2:88047/1‑155 (MQ=255)
                                                                                                                      cttgatcttgatcttgatcttgatcttgatcttcgcgacttcaggaggccgaacgcaggccttgcggaggcaggtgacgggcgtCCATGCCCGTCACCTCCCTCCGCACGCCCTGCGCGCGACCTCCTGAAGTc                                                                                                                                                                                         >  1:431932‑M2/85‑134 (MQ=255)
                                                                                                                          ggggTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGaa                                                                                                                                                     >  2:757456/1‑165 (MQ=255)
                                                                                                                                cccTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCGGGATCTTGc                                                                                <  2:926090/228‑1 (MQ=255)
                                                                                                                                atcttgatcttgatcttgatcttcgcgacttcaggaggccgaacgcaggccttgcggaggcaggtgacgggcatCCAGGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGaa                                                                                                                                                     <  2:649037‑M2/86‑1 (MQ=255)
                                                                                                                                    cGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCGGGATCTTGCCTTGATTTTAg                                                                     >  1:907156/1‑235 (MQ=255)
                                                                                                                                    cGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCGGGATCTTGCCTTGATTTTAg                                                                     >  1:848866/1‑235 (MQ=255)
                                                                                                                                     gatcttgatcttgatcttcgcgacttcaggaggccgaacgcaggccttgcggaggcaggtgacgggcatCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACaa                                                                                                                                                     <  2:383571‑M2/86‑1 (MQ=255)
                                                                                                                                                    ccGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCGGGATCTTGCCTTGATTTTAGCTCAATGCCCCGCCGAt                                                    >  1:666174/1‑236 (MQ=255)
                                                                                                                                                     ttcgcgacttcaggaggccgaacgcaggccttgcggaggcaggtgacgggcatGGATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCGGGATCTTGCCTTGATTTTAGCTCAATGCCccgccg                                                      <  2:749465‑M2/181‑1 (MQ=255)
                                                                                                                                                          gacttcaggaggccgaacgcaggccttgcggaggcaggtgacgggcatCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCg                                                                                         >  2:743018‑M2/49‑194 (MQ=255)
                                                                                                                                                           gaCCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCGGGATCTTGCCTTGATTTTAGCTCAATGCCCCGCCGATGCCggg                                              <  1:365217/235‑1 (MQ=255)
                                                                                                                                                             ccGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCGGGATCTTGCCTTGATTTTAGCTCACTGCCCCGCCGTTGCCGGGcc                                            <  2:492624/235‑1 (MQ=255)
                                                                                                                                                                 gcTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCGGGATCTTGCCTTGATTTTAGATCAATGCCCCGCCGATGCCGGGCCCGCCt                                       >  1:564366/1‑236 (MQ=255)
                                                                                                                                                                  gaggccgaacgcaggccttgcggaggcaggtgacgggcatGGATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCGGGATCTTGCCTTGATTTTAGCTCAATGCCCCGCCGATGcc                                                 <  2:434872‑M2/186‑1 (MQ=255)
                                                                                                                                                                   cgaaGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATGGATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCCTAGACGCAAAATCAAAAGCCGAagc                                                                                                                                 >  2:38635/4‑145 (MQ=17)
                                                                                                                                                                      aGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCGGGATCTTGCCTTGATTTTAGCTCAATGCCCCGCCGATGCCGGGCCCGCCTTAGCg                                  <  1:432682/236‑1 (MQ=255)
                                                                                                                                                                        gCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagagcgagagcgagagcgagagcatcgg                                                                                              >  2:433908/1‑169 (MQ=255)
                                                                                                                                                                        gCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCGGGATCTTGCCTTGATTTTAGCTCAATGCCCCGCCGATGCCGGGCCCGCCTTAGCg                                  >  2:432679/1‑234 (MQ=255)
                                                                                                                                                                          gcgcagggcttgcgaagggaggtgacgggcatCCATGCCCGTCCCCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAg                                                                                                                                        >  2:867478‑M2/33‑131 (MQ=255)
                                                                                                                                                                                           gcgcTTGACGGGCATGGATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCGGGATCTTGCCTTGATTTTAGCTCAATGCCCCGCCGATGCCGGGCCCGCCTTAGCGGTAAAGGGCGGCTTGGCGa               >  1:862825/1‑235 (MQ=255)
                                                                                                                                                                                            cgcTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGCGCGAGAGCGAGAGCAt                                                                                           >  2:291704/1‑158 (MQ=255)
                                                                                                                                                                                             gcTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCGGGATCTTGCCTTGATTTTAGCTCAATGCCCCGCCGATGCCGGGCCCGCCTTAGCGGTAAAGGGCGGCTTGGCGacc             >  1:167875/1‑233 (MQ=255)
                                                                                                                                                                                              cTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATc                                                                                                                                                                 <  1:36559/86‑1 (MQ=255)
                                                                                                                                                                                              cTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATc                                                                                                                                                                 >  2:36559/1‑86 (MQ=255)
                                                                                                                                                                                                        atCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCGGGATCTTGCCTTGATTTTAGCTCAATGCCCCGCCGATGCCGGGCCCGCCTTAGCGGTAAAGGGCGGCTTGGCGAACCACACcagcag  <  1:828702‑M2/233‑1 (MQ=255)
                                                                                                                                                                                                          |                                                                                                                                                                                                                                        
TCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCGGGATCTTGCCTTGATTTTAGCTCAATGCCCCGCCGATGCCGGGCCCGCCTTAGCGGTAAAGGGCGGCTTGGCGAACCACACCAGCAG  >  AE015451/4022605‑4023038

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