Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,200,983 A→G 5.0% S177P (TCG→CCG) ‡ ybgJ ← allophanate hydrolase subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,200,9830AG 5.0% 1521.9 / 40.0 455S177P (TCG→CCG) ‡ybgJallophanate hydrolase subunit 1
Reads supporting (aligned to +/- strand):  ref base A (201/230);  new base G (11/12);  total (212/243)
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.50e-01

Only 100 of 459 total aligned reads displayed.

CTCCACATCACCACCCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGCAAAC  >  NC_002947/5200834‑5201131
                                                                                                                                                     |                                                                                                                                                    
   cacaTCACCACCCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATAc                                                                                                                                                   >  2:447350/1‑150 (MQ=255)
    acaTCACCACCCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACt                                                                                                                                                  >  1:6778348/1‑150 (MQ=255)
         accaccCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCgg                                                                                                                                             >  1:4929832/1‑150 (MQ=255)
         accaccCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCgg                                                                                                                                             <  1:3782326/150‑1 (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:6499877/1‑150 (MQ=255)
            accCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTa                                                                                                                                          <  1:6997060/150‑1 (MQ=255)
            accCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTa                                                                                                                                          >  2:4888875/1‑150 (MQ=255)
             cccATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCAGTATCGCCGGGtgg                                                                                                                                         <  1:2710176/150‑3 (MQ=255)
             cccATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCAGTATCGCCGGGtgg                                                                                                                                         <  1:3401622/150‑3 (MQ=255)
              ccATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAgg                                                                                                                                        >  1:3142965/1‑150 (MQ=255)
               cctGGCCAGGAATTCGCTGCGCGAAACCGGCACGCAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCAGTATCGCCGGGtggct                                                                                                                                       <  2:321636/148‑5 (MQ=255)
               cATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCTACCCGGCAGTATCGCCGGGtggct                                                                                                                                       <  1:138761/150‑5 (MQ=255)
                  ggCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGgccgc                                                                                                                                    <  1:7428913/150‑1 (MQ=255)
                   gCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGgccgcc                                                                                                                                   <  1:6601019/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)
                          aTTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGgc                                                                                                                            <  2:5091711/150‑1 (MQ=255)
                              gCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTc                                                                                                                        >  1:4098249/1‑150 (MQ=255)
                                tGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCgg                                                                                                                      >  2:3500088/1‑150 (MQ=255)
                                tGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCgg                                                                                                                      <  2:2397421/150‑1 (MQ=255)
                                tGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCgg                                                                                                                      >  1:280796/1‑150 (MQ=255)
                                tGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCgg                                                                                                                      <  1:527535/150‑1 (MQ=255)
                                  cgcgAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCa                                                                                                                    <  2:3697916/150‑1 (MQ=255)
                                    cgAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAAt                                                                                                                  >  1:3496883/1‑150 (MQ=255)
                                    cgAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAAt                                                                                                                  <  1:580557/150‑1 (MQ=255)
                                         ccGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCa                                                                                                             <  2:5954709/150‑1 (MQ=255)
                                         ccGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCa                                                                                                             <  1:6205321/150‑1 (MQ=255)
                                          cGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACCACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATCCTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCACTACCCAc                                                                                                            >  2:1061676/1‑150 (MQ=255)
                                           ggCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGTCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACg                                                                                                           >  1:2385675/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)
                                                             tCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACgcgc                                                                                         <  1:2310712/150‑1 (MQ=255)
                                                                    gTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTg                                                                                  <  1:6359853/150‑1 (MQ=255)
                                                                    gTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTg                                                                                  <  1:190561/150‑1 (MQ=255)
                                                                      tGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTggg                                                                                <  2:2526121/150‑1 (MQ=255)
                                                                       gCAAGAGGGTGTAACCGTCACGCTCACGGTCGAAGATCCGCAGCGGGGTGCGGCCGATCTGGTTCCAGCCACGCGGCGATACTGGGGGGTAGGCCGCCGTCTGGCGGTCGGCAATTCCCTCGCTGACCGCCGCCACGCGCTTGCGTgggt                                                                               >  2:6629392/1‑149 (MQ=255)
                                                                       gCAAGAGGCTGTAACCCTCCCGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTgggg                                                                               <  2:4258207/150‑1 (MQ=255)
                                                                       gCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTgggg                                                                               >  2:6516553/1‑150 (MQ=255)
                                                                         aaGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTg                                                                                                                              <  1:4222109/101‑1 (MQ=255)
                                                                         aaGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTa                                                                             <  2:5635727/150‑1 (MQ=255)
                                                                          agagGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTAc                                                                            >  1:88622/1‑150 (MQ=255)
                                                                             ggCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCa                                                                         <  2:5474382/150‑1 (MQ=255)
                                                                                tGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGAc                                                                      >  2:2755462/1‑150 (MQ=255)
                                                                                tGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGAc                                                                      >  1:1745449/1‑150 (MQ=255)
                                                                                   aaCCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAg                                                                   <  2:262577/150‑1 (MQ=255)
                                                                                    aCCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAgg                                                                  <  2:539744/150‑1 (MQ=255)
                                                                                      ccTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGc                                                                >  1:545762/1‑150 (MQ=255)
                                                                                      ccTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGc                                                                >  1:5012448/1‑150 (MQ=255)
                                                                                      ccTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGc                                                                >  1:4263311/1‑150 (MQ=255)
                                                                                       ctcacgctcacgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCt                                                               >  2:1615295/1‑150 (MQ=255)
                                                                                        tcacgctcacgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTg                                                              >  1:724620/1‑150 (MQ=255)
                                                                                        tcacgctcacgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTg                                                              <  2:1269012/150‑1 (MQ=255)
                                                                                           cgctcacgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGcc                                                           >  1:2514459/1‑150 (MQ=255)
                                                                                             ctcacgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAg                                                         <  1:2169548/150‑1 (MQ=255)
                                                                                             ctcacgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCGGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAg                                                         <  2:6879533/150‑1 (MQ=255)
                                                                                              tcacgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAgg                                                        >  2:7402949/1‑150 (MQ=255)
                                                                                                acgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGgcg                                                      >  1:309092/1‑150 (MQ=255)
                                                                                                acgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGgcg                                                      >  1:950081/1‑150 (MQ=255)
                                                                                                    tCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCtcg                                                  >  1:6978807/1‑150 (MQ=255)
                                                                                                      gAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCtcgtc                                                >  2:2953479/1‑150 (MQ=255)
                                                                                                      gAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCtcgtc                                                <  1:837640/150‑1 (MQ=255)
                                                                                                       aacaacCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCtcgtcg                                               >  1:1204/1‑150 (MQ=255)
                                                                                                       aacaacCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCtcgtcg                                               >  1:5869016/1‑150 (MQ=255)
                                                                                                        acaacCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGa                                              >  2:1818063/1‑150 (MQ=255)
                                                                                                          gggcGCACCCCGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGAcc                                            >  1:5575960/4‑150 (MQ=255)
                                                                                                          aacCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGAcc                                            >  2:3855643/1‑150 (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)
                                                                                                               cACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTgg                                                                                 >  2:1427340/1‑108 (MQ=255)
                                                                                                               cACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGcc                                       >  1:1007889/1‑150 (MQ=255)
                                                                                                                aCCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGccc                                      <  2:4634765/150‑1 (MQ=255)
                                                                                                                aCCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGccc                                      >  1:4808239/1‑150 (MQ=255)
                                                                                                                  cGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCAt                                    <  1:2309021/150‑1 (MQ=255)
                                                                                                                    gggTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGa                                  <  1:5677146/150‑1 (MQ=255)
                                                                                                                       tGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGc                               <  2:379772/150‑1 (MQ=255)
                                                                                                                        gcgcCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGcc                              <  2:7223217/150‑1 (MQ=255)
                                                                                                                                  aGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGt                    >  2:4188258/1‑150 (MQ=255)
                                                                                                                                  aGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGAGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGt                    <  1:368220/150‑1 (MQ=255)
                                                                                                                                   ggTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTg                   >  2:3556315/1‑150 (MQ=255)
                                                                                                                                       ccAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGaa               <  1:447350/150‑1 (MQ=255)
                                                                                                                                        cAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCACGTCGACCAGGCCCAAGAAGCCGAAACCCGGTGCGAAg              <  2:5836368/150‑1 (MQ=255)
                                                                                                                                         aGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGc             >  2:4004416/1‑150 (MQ=255)
                                                                                                                                          gCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGcc            >  1:4716476/1‑150 (MQ=255)
                                                                                                                                           ccACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGccc           <  1:7059057/150‑1 (MQ=255)
                                                                                                                                            taCCCGGCAGTATCGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCa          >  1:4666555/2‑150 (MQ=255)
                                                                                                                                            taCCCGGCAGTATCGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCa          >  2:1954382/2‑150 (MQ=255)
                                                                                                                                            cACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCa          >  1:2905901/1‑150 (MQ=255)
                                                                                                                                             aCCCGGCGATTCTGCCGGGTAGGCCGCAGTCTGGCGCTCGGCAATACCCACGCTGCCCGCAGCCACGCGCTTGCGTGGGGTACTCAGACGTGGGCTGGACAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAg         <  1:3877040/150‑1 (MQ=255)
                                                                                                                                              cccGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAgg        <  2:472300/150‑1 (MQ=255)
                                                                                                                                                cGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGc      >  2:2223317/1‑150 (MQ=255)
                                                                                                                                                 ggCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAg              <  1:1616215/141‑1 (MQ=255)
                                                                                                                                                   cGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCACTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGCaaa   <  1:3988636/150‑1 (MQ=255)
                                                                                                                                                    gATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGCAAac  >  1:4539743/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTCCACATCACCACCCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGGGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGCAAAC  >  NC_002947/5200834‑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: