Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,200,982 G→A 5.1% S177L (TCG→TTG) ‡ ybgJ ← allophanate hydrolase subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,200,9820GA 5.1% 1487.7 / 42.4 454S177L (TCG→TTG) ‡ybgJallophanate hydrolase subunit 1
Reads supporting (aligned to +/- strand):  ref base G (198/233);  new base A (11/12);  total (209/245)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01

Only 100 of 461 total aligned reads displayed.

CCTCCACATCACCACCCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGCAAAC  >  NC_002947/5200833‑5201131
                                                                                                                                                     |                                                                                                                                                      
  tCCACATCACCACCCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGAt                                                                                                                                                      >  1:5506626/1‑150 (MQ=255)
  tCCACATCACCACCCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGAt                                                                                                                                                      >  1:563237/1‑150 (MQ=255)
  tCCACATCACCACCCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGAt                                                                                                                                                      >  1:3721720/1‑150 (MQ=255)
          accaccCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCgg                                                                                                                                              >  1:4929832/1‑150 (MQ=255)
          accaccCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCgg                                                                                                                                              <  1:3565295/150‑1 (MQ=255)
          accaccCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGAGCAGGTTCCAGCCACCCGGCGATACTGCCgg                                                                                                                                              >  2:5543000/1‑150 (MQ=255)
           ccaccCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCggg                                                                                                                                             >  1:207811/1‑150 (MQ=255)
            caccCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGt                                                                                                                                            >  1:231559/1‑150 (MQ=255)
             accCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTa                                                                                                                                           >  2:4888875/1‑150 (MQ=255)
              cccATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCAGTATCGCCGGGtgg                                                                                                                                          <  1:3401622/150‑3 (MQ=255)
              cccATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCAGTATCGCCGGGtgg                                                                                                                                          <  1:2710176/150‑3 (MQ=255)
               ccATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAgg                                                                                                                                         >  1:3142965/1‑150 (MQ=255)
                cATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCTACCCGGCAGTATCGCCGGGtggct                                                                                                                                        <  1:3786370/150‑5 (MQ=255)
                cATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCTACCCGGCAGTATCGCCGGGtggct                                                                                                                                        <  1:138761/150‑5 (MQ=255)
                    gCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGgccgcc                                                                                                                                    <  1:901846/150‑1 (MQ=255)
                     ccAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGccgccg                                                                                                                                   >  2:333931/1‑150 (MQ=255)
                       aGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTc                                                                                                                                 >  1:3268934/1‑150 (MQ=255)
                        ggAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCt                                                                                                                                <  2:5218359/150‑1 (MQ=255)
                        ggAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCt                                                                                                                                <  2:2343843/150‑1 (MQ=255)
                           aTTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGgc                                                                                                                             <  2:5091711/150‑1 (MQ=255)
                               gCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTc                                                                                                                         >  1:4098249/1‑150 (MQ=255)
                                 tGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCgg                                                                                                                       <  1:527535/150‑1 (MQ=255)
                                 tGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCgg                                                                                                                       >  1:280796/1‑150 (MQ=255)
                                  gcgcGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGc                                                                                                                      <  1:4550616/150‑1 (MQ=255)
                                   cgcgAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCa                                                                                                                     <  1:2977103/150‑1 (MQ=255)
                                      gAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATa                                                                                                                  <  1:6529406/150‑1 (MQ=255)
                                          ccGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCa                                                                                                              <  2:5954709/150‑1 (MQ=255)
                                           cGGCACGAAGCGCACCCGGTCCCCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGCACAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCAc                                                                                                             <  2:3132431/150‑1 (MQ=255)
                                            ggCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGTCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACg                                                                                                            >  1:2385675/1‑150 (MQ=255)
                                                 gAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGcc                                                                                                       >  2:141520/1‑150 (MQ=255)
                                                     cgcACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCccgcc                                                                                                   <  2:7008846/150‑1 (MQ=255)
                                                       cACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCcgccgc                                                                                                 <  2:5773974/150‑1 (MQ=255)
                                                        aCCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCgccgcc                                                                                                >  2:2093743/1‑150 (MQ=255)
                                                          ccGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCAc                                                                                              <  1:4731054/150‑1 (MQ=255)
                                                                     gTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTg                                                                                   <  1:190561/150‑1 (MQ=255)
                                                                     gTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTg                                                                                   <  1:4264502/150‑1 (MQ=255)
                                                                       tGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTggg                                                                                 <  2:2526121/150‑1 (MQ=255)
                                                                       tGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTggg                                                                                 <  1:5045475/150‑1 (MQ=255)
                                                                       tGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTggg                                                                                 <  1:3522670/150‑1 (MQ=255)
                                                                        gCAAGAGGGTGTAACCGTCACGCTCACGGTCGAAGATCCGCAGCGGGGTGCGGCCGATCTGGTTCCAGCCACGCGGCGATACTGGGGGGTAGGCCGCCGTCTGGCGGTCGGCAATTCCCTCGCTGACCGCCGCCACGCGCTTGCGTgggt                                                                                >  2:6629392/1‑149 (MQ=255)
                                                                        gCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTgggg                                                                                >  2:6516553/1‑150 (MQ=255)
                                                                         cAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGt                                                                               >  1:968018/1‑150 (MQ=255)
                                                                          aaGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTa                                                                              <  2:5635727/150‑1 (MQ=255)
                                                                          aaGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTa                                                                              >  1:6727699/1‑150 (MQ=255)
                                                                            gagGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACt                                                                            >  1:484250/1‑150 (MQ=255)
                                                                            gagGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACt                                                                            >  1:6370399/1‑150 (MQ=255)
                                                                              ggCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCa                                                                          <  2:5474382/150‑1 (MQ=255)
                                                                                 tGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGAc                                                                       >  2:2755462/1‑150 (MQ=255)
                                                                                 tGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGAc                                                                       >  1:6879054/1‑150 (MQ=255)
                                                                                 tGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGAc                                                                       >  1:1745540/1‑150 (MQ=255)
                                                                                 tGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGAc                                                                       >  1:7441868/1‑150 (MQ=255)
                                                                                    aaCCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAg                                                                    <  2:262471/150‑1 (MQ=255)
                                                                                        ctcacgctcacgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTa                                                                              <  1:580857/136‑1 (MQ=255)
                                                                                        ctcacgctcacgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTa                                                                              >  2:580857/1‑136 (MQ=255)
                                                                                        ctcacgctcacgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCt                                                               >  2:3460641/1‑150 (MQ=255)
                                                                                            cgctcacgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGcc                                                           >  1:2514459/1‑150 (MQ=255)
                                                                                             gctcacgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCa                                                          >  1:6072180/1‑150 (MQ=255)
                                                                                              ctcacgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAg                                                         <  1:2169548/150‑1 (MQ=255)
                                                                                               tcacgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAgg                                                        >  2:7402949/1‑150 (MQ=255)
                                                                                                 acgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGgcg                                                      >  1:309092/1‑150 (MQ=255)
                                                                                                    gTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGCGCCTGGCCAgccctc                                                    <  2:4236872/150‑6 (MQ=255)
                                                                                                       gaccacccGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATCCTGCCGGGTCGGCCGCCGGCTGGCGCGCGGCAATACACGCGCTGCCCGCCGCCCCGCGCGTGCGGGGGGTACGCCGACGCGCCCAGGACAGGCGCtcgtc                                                >  1:7124779/7‑150 (MQ=255)
                                                                                                       gAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCtcgtc                                                <  1:837647/150‑1 (MQ=255)
                                                                                                       gAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCtcgtc                                                >  2:1784206/1‑150 (MQ=255)
                                                                                                       gAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCtcgtc                                                <  2:2858175/150‑1 (MQ=255)
                                                                                                        aacaacCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCtcgtcg                                               >  1:1204/1‑150 (MQ=255)
                                                                                                        aacaacCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCtcgtcg                                               >  1:7303984/1‑150 (MQ=255)
                                                                                                           aacCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGAcc                                            >  2:3855643/1‑150 (MQ=255)
                                                                                                            acCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTggg                                                                                 >  2:4707438/1‑113 (MQ=255)
                                                                                                               gCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGc                                        >  2:514451/1‑150 (MQ=255)
                                                                                                               gCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGc                                        >  2:853211/1‑150 (MQ=255)
                                                                                                               gCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGc                                        <  2:6118913/150‑1 (MQ=255)
                                                                                                                cACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTgg                                                                                  <  1:1427500/108‑1 (MQ=255)
                                                                                                                cACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGcc                                       >  1:1007889/1‑150 (MQ=255)
                                                                                                                 aCCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGccc                                      <  2:4634765/150‑1 (MQ=255)
                                                                                                                  ccGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCa                                     <  1:3989676/150‑1 (MQ=255)
                                                                                                                   cGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCAt                                    <  1:2309021/150‑1 (MQ=255)
                                                                                                                     gggTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGa                                  <  1:5677146/150‑1 (MQ=255)
                                                                                                                         gcgcCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGcc                              <  2:7223217/150‑1 (MQ=255)
                                                                                                                             ccGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGaaa                          >  2:1547307/1‑150 (MQ=255)
                                                                                                                             ccGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGaaa                          <  1:7002527/150‑1 (MQ=255)
                                                                                                                              cGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAAc                         <  2:3618765/150‑1 (MQ=255)
                                                                                                                              cGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAAc                         >  2:4465751/1‑150 (MQ=255)
                                                                                                                              cGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAAc                         >  2:2117144/1‑150 (MQ=255)
                                                                                                                                 tCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCg                      <  2:379453/150‑1 (MQ=255)
                                                                                                                                  cAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCgg                     <  2:1297388/150‑1 (MQ=255)
                                                                                                                                   aGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGt                    >  2:4188258/1‑150 (MQ=255)
                                                                                                                                    ggTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTg                   >  2:3556315/1‑150 (MQ=255)
                                                                                                                                        ccAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGaa               <  1:447350/150‑1 (MQ=255)
                                                                                                                                          aGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGc             <  1:6207910/150‑1 (MQ=255)
                                                                                                                                            ctaCCCGGCAGTATCGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGccc           >  2:1099199/3‑150 (MQ=255)
                                                                                                                                            ccACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGccc           <  2:5383371/150‑1 (MQ=255)
                                                                                                                                             taCCCGGCAGTATCGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCa          >  1:4666555/2‑150 (MQ=255)
                                                                                                                                             cACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCa          <  1:6981795/150‑1 (MQ=255)
                                                                                                                                              aCCCGGCGATTCTGCCGGGTAGGCCGCAGTCTGGCGCTCGGCAATACCCACGCTGCCCGCAGCCACGCGCTTGCGTGGGGTACTCAGACGTGGGCTGGACAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAg         <  1:3877040/150‑1 (MQ=255)
                                                                                                                                               cccGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAgg        <  2:3995433/150‑1 (MQ=255)
                                                                                                                                                ccGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAggg       >  1:4231375/1‑150 (MQ=255)
                                                                                                                                                ccGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAggg       >  2:1946912/1‑150 (MQ=255)
                                                                                                                                                  ggCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAg              <  1:1616215/141‑1 (MQ=255)
                                                                                                                                                     gATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGCAAac  >  2:6830784/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CCTCCACATCACCACCCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGCAAAC  >  NC_002947/5200833‑5201131

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: