Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,767,507 G→A 5.2% intergenic (+18/‑123) gpmI → / → PP_5057 2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase/M23/M37 family peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,767,5070GA 5.2% 680.0 / 16.9 211intergenic (+18/‑123)gpmI/PP_50572,3‑bisphosphoglycerate‑independent phosphoglycerate mutase/M23/M37 family peptidase
Reads supporting (aligned to +/- strand):  ref base G (87/113);  new base A (3/8);  total (90/121)
Fisher's exact test for biased strand distribution p-value = 3.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

Only 100 of 215 total aligned reads displayed.

TACGTGGGCAAGCGCAACGTGAAAGTGCGCGAGGGCGGCGTGCTGGCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGCCAGTCTCCATCCTTGGTACGCCAAACCCCATGCTTCGCGCCCTGATCCTCCTAGCC  >  NC_002947/5767358‑5767656
                                                                                                                                                     |                                                                                                                                                     
tACGTGGGCAAGCGCAACGTGAAAGTGCGCGAGGGCGGCGTGCTGGCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCg                                                                                                                                                       >  1:1255680/1‑150 (MQ=255)
tACGTGGGCAAGCGCAACGTGAAAGTGCGCGAGGGCGGCGTGCTGGCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCg                                                                                                                                                       <  2:31992/150‑1 (MQ=255)
 aCGTGGGCAAGCGCAACGTGAAAGTGCGCGAGGGCGGCGTGCTGGCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGt                                                                                                                                                      >  2:50683/1‑150 (MQ=255)
   gTGGGCAAGCGCAACGTGAAAGTGCGCGAGGGCGGCGTGCTGGCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTg                                                                                                                                                    <  1:3064860/150‑1 (MQ=255)
     gggCAAGCGCAACGTGAAAGTGCGCGAGGGCGGCGTGCTGGCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAACTGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCa                                                                                                                                                  <  1:5583612/150‑1 (MQ=255)
     gggCAAGCGCAACGTGAAAGTGCGCGAGGGCGGCGTGCTGGCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCa                                                                                                                                                  <  2:4377776/150‑1 (MQ=255)
       gCAAGCGCAACGTGAAAGTGCGCGAGGGCGGCGTGCTGGCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGc                                                                                                                                                >  1:5107177/1‑150 (MQ=255)
             gcAACGTGAAAGTGCGCGAGGGCGGCGTGCTGGCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCAACGCAGccccttt                                                                                                                                          <  2:4062389/150‑4 (MQ=255)
             gcAACGTGAAAGTGCGCGAGGGCGGCGTGCTGGCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCAACGCAGccccttt                                                                                                                                          <  2:4029876/150‑4 (MQ=255)
             gcAACGTGAAAGTGCGCGAGGGCGGCGTGCTGGCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCAACGCAGccccttt                                                                                                                                          <  1:3998584/150‑4 (MQ=255)
              cAACGTGAAAGTGCGCGAGGGCGGCGTGCTGGCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAAtt                                                                                                                                         >  1:527108/1‑150 (MQ=255)
               aaCGTGAAAGTGCGCGAGGGCGGCGTGCTGGCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCCATTg                                                                                                                                        >  1:781585/1‑150 (MQ=255)
               aaCGTGAAAGTGCGCGAGGGCGGCGTGCTGGCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTg                                                                                                                                        <  2:3546441/150‑1 (MQ=255)
               aaCGTGAAAGTGCGCGAGGGCGGCGTGCTGGCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTg                                                                                                                                        <  1:108398/150‑1 (MQ=255)
                 cGTGAAAGTGCGCGAGGGCGGCGTGCTGGCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGca                                                                                                                                      <  1:2828193/150‑1 (MQ=255)
                    gAAAGTGCGCGAGGGCGGCGTGCTGGCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACAt                                                                                                                                   <  1:2700241/150‑1 (MQ=255)
                      aaGTGCGCGAGGGCGGCGTGCTGGCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGa                                                                                                                                 <  2:2870973/150‑1 (MQ=255)
                            gcgAGGGCGGCGTGCTGGCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCTGTAAATCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCa                                                                                                                           <  2:1063769/150‑1 (MQ=255)
                            gcgAGGGCGGCGTGCTGGCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCa                                                                                                                           <  2:2547498/150‑1 (MQ=255)
                            gcgAGGGCGGCGTGCTGGCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCa                                                                                                                           <  1:5251102/150‑1 (MQ=255)
                              gAGGGCGGCGTGCTGGCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGa                                                                                                                         <  2:4811744/150‑1 (MQ=255)
                                    ggcgTGCTGGCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGTATGTCGATAGGGGCTGCGTTGCTGCCCCAATTGCACATGAATTCCTGACATACg                                                                                                                   >  2:3443380/1‑150 (MQ=255)
                                     gcgTGCTGGCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGc                                                                                                                  >  2:2608407/1‑150 (MQ=255)
                                       gTGCTGGCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGccc                                                                                                                >  1:3168370/1‑150 (MQ=255)
                                         gCTGGCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGccccc                                                                                                              <  1:4825632/150‑1 (MQ=255)
                                         gCTGGCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGccccc                                                                                                              <  1:1026596/150‑1 (MQ=255)
                                           tGGCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCt                                                                                                            <  2:2271518/150‑1 (MQ=255)
                                           tGGCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCt                                                                                                            <  1:3143375/150‑1 (MQ=255)
                                            ggCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTg                                                                                                           >  2:4286812/1‑150 (MQ=255)
                                             gCCGACGTGGCGCCAACCATGCTGACGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGc                                                                                                          >  2:3390583/1‑150 (MQ=255)
                                               cGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCAt                                                                                                        >  1:87706/1‑150 (MQ=255)
                                                        gcCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGa                                                                                               >  1:1124403/1‑150 (MQ=255)
                                                        gcCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGa                                                                                               <  2:1242599/150‑1 (MQ=255)
                                                        gcCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGa                                                                                               <  2:733609/150‑1 (MQ=255)
                                                         ccAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAg                                                                                              <  1:5081457/150‑1 (MQ=255)
                                                           aaCCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACCTGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGc                                                                                            >  2:2011/1‑150 (MQ=255)
                                                             ccATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGt                                                                                          <  1:5556166/150‑1 (MQ=255)
                                                                  cTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTTTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCTCGCTGCGGATGCAGGGGGGCGtttgtc                                                                                     <  1:5260587/150‑4 (MQ=255)
                                                                  cTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGtttttt                                                                                     <  1:4689161/150‑1 (MQ=255)
                                                                    gAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGtttttttt                                                                                   >  2:4507338/1‑150 (MQ=255)
                                                                      aGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTtgc                                                                                 <  2:3033495/150‑1 (MQ=255)
                                                                      aGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTtgc                                                                                 <  2:137863/150‑1 (MQ=255)
                                                                        ctgctgGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTtgcct                                                                               <  1:2782779/150‑1 (MQ=255)
                                                                           ctgGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTtgcctgcc                                                                            <  2:2093735/150‑1 (MQ=255)
                                                                            tgGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGccc                                                                           <  2:1732202/150‑1 (MQ=255)
                                                                              gggCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAg                                                                         <  2:1459025/150‑1 (MQ=255)
                                                                                gCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGc                                                                       <  2:4263478/150‑1 (MQ=255)
                                                                                  tGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCgg                                                                     <  1:2507629/150‑1 (MQ=255)
                                                                                  tGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCgg                                                                     <  1:759000/150‑1 (MQ=255)
                                                                                  tGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCgg                                                                     <  1:3873849/150‑1 (MQ=255)
                                                                                  tGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCgg                                                                     <  1:2448121/150‑1 (MQ=255)
                                                                                   ggAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCggg                                                                    <  1:330577/150‑1 (MQ=255)
                                                                                   ggAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCggg                                                                    >  2:1179277/1‑150 (MQ=255)
                                                                                    gAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGc                                                                   <  2:5157525/150‑1 (MQ=255)
                                                                                     aacaGCCGGTGGAAATGACCGGTACCTCGATTCGGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAACCGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCa                                                                  <  1:5625193/147‑1 (MQ=255)
                                                                                      aaaGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGcct                                                                 >  2:3946215/1‑148 (MQ=255)
                                                                                       aaGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATa                                                                <  1:4504863/150‑1 (MQ=255)
                                                                                         gCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCCTACt                                                              >  2:5349857/1‑150 (MQ=255)
                                                                                         gCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACt                                                              >  2:3949624/1‑150 (MQ=255)
                                                                                          ccGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTa                                                             >  1:3117242/1‑150 (MQ=255)
                                                                                           cGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAg                                                            <  1:360316/150‑1 (MQ=255)
                                                                                             gTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGc                                                          <  1:1498085/150‑1 (MQ=255)
                                                                                             gTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGc                                                          <  1:3467889/150‑1 (MQ=255)
                                                                                             gTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGc                                                          <  1:3468052/150‑1 (MQ=255)
                                                                                             gTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCCGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGc                                                          >  2:5330566/1‑150 (MQ=255)
                                                                                              tGGAAATGACCGGTACCTCGAGTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGcc                                                         <  2:1349894/150‑1 (MQ=255)
                                                                                               ggAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGCCa                                                        <  2:4221556/150‑1 (MQ=255)
                                                                                                 aaaTGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGCCAGt                                                      <  2:2760420/150‑1 (MQ=255)
                                                                                                      aCCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGCCAGTCTCCa                                                 <  2:1024471/150‑1 (MQ=255)
                                                                                                       ccGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGCCAGTCTCCaa                                                >  2:1952515/1‑149 (MQ=255)
                                                                                                        cGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGCCAGTCTCCATc                                               <  2:5251515/150‑1 (MQ=255)
                                                                                                              cTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTTCCTGCCCCGGCGGGCATACTAGGCCAGTCTCCAACCATGGt                                         >  1:2378024/1‑150 (MQ=255)
                                                                                                              cTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTTCCTGCCCAGGCGGGCATACTAGGCCCGTCTCCAACCTTGGt                                         >  1:999961/1‑150 (MQ=255)
                                                                                                              cTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTTCCTGCCCAGGCGGGCATACAAGGCCAGTCACCAACCTTGGt                                         >  1:2561854/1‑150 (MQ=255)
                                                                                                               tCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTTCCTGCCCAGGCGGGCATACTAGGCCAGTCTCCATCCTTGGTa                                        >  1:4309405/1‑150 (MQ=255)
                                                                                                               tCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGCCAGTCTCCATCCTTGGTa                                        >  2:366842/1‑150 (MQ=255)
                                                                                                               tCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGCCAGTCTCCATCCGTGGTa                                        >  1:5370120/1‑150 (MQ=255)
                                                                                                               tCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGCCAGTCTCCATCCGTGGTa                                        >  1:3799391/1‑150 (MQ=255)
                                                                                                                  aTTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGCCAGTCTCCATCCTTGGTACGc                                     <  1:4872850/150‑1 (MQ=255)
                                                                                                                      tGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGCCAGTCTCCATCCTTGGTACGCCaaa                                 <  2:5103783/150‑1 (MQ=255)
                                                                                                                          cGACGCCTGAATGACGAAAGGCGCTGCGTTGCAGCCCCAATAGCACATGAATTCCACACAAACGCCCCCCACCATCCGCAGCGAGGCGAATTATTTGCCTGCCCAGGCGGGCATACTAGGCCAGTCTCCATCCTTGGTACGCCAAAcccc                             <  2:5561180/150‑1 (MQ=255)
                                                                                                                           gACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGCCAGTCTCCATCCTTGGTACGCCAAACCCCa                            <  2:877024/150‑1 (MQ=255)
                                                                                                                           gACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGCCAGTCTCCATCCTTGGTACGCCAAACCCCa                            <  2:877144/150‑1 (MQ=255)
                                                                                                                           gACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGCCAGTCTCCATCCTTGGTACGCCAAACCCCa                            >  2:3758269/1‑150 (MQ=255)
                                                                                                                           gACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGCCAGTCTCCATCCTTGGTACGCCAAACCCCa                            >  2:3340775/1‑150 (MQ=255)
                                                                                                                            aCGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGCCAGTCTCCATCCATGGTACGCCAAACCCCaa                           >  1:5207347/1‑149 (MQ=255)
                                                                                                                              gCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGCCAGTCTCCATCCTTGGTACGCCAAACCCCATGc                         <  2:1600370/150‑1 (MQ=255)
                                                                                                                              gCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCAGCCCAGGCCGGCAAACAAGGCCCCTCTCCATCCTTGGTACGCCCAACCCCAAGc                         >  1:1961675/1‑150 (MQ=255)
                                                                                                                                cTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGCCAGTCTCCATCCTTGGTACGCCAAACCCCATGCtt                       >  2:5109744/1‑150 (MQ=255)
                                                                                                                                 tGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGCCAGTCTCCATCCTTGGTACGCCAAACCCCATGCTTc                      <  2:5227583/150‑1 (MQ=255)
                                                                                                                                        caattggggCTGCAACGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGCCAGTCTCCATCCTTGGTACGCCAAACCCCATGCTTCGCGCCCt               >  2:4574373/6‑150 (MQ=255)
                                                                                                                                         gAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGCCAGTCTCCATCCTTGGTACGCCAAACCCCATGCTTCGCGCCCTg              >  2:3066677/1‑150 (MQ=255)
                                                                                                                                          attggggCTGCAACGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGCCAGTCTCCATCCTTGGTACGCCAAACCCCATGCTTCGCGCCCTGa             <  1:4574373/147‑1 (MQ=255)
                                                                                                                                          aaaGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGGCCAGGCGGGCATACTAGGCCAGTCTCCATCCTTGGTAAGCCAAAAACCATTCTTCGCGCCCTGa             >  1:2023966/1‑150 (MQ=255)
                                                                                                                                          aaaGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGCCAGTCTCCATCCTTGGTACGCCAAACCCCATGCTTCGCGCCCTGa             <  2:5049516/150‑1 (MQ=255)
                                                                                                                                          aaaGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGCCAGTCTCCATCCTTGGTACGCCAAACCCCATGCTTCGCGCCCTGa             <  1:1664573/150‑1 (MQ=255)
                                                                                                                                               ggCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGCCAGTCTCCATCCTTGGTACGCCAAACCCCATGCTTCGCGCCCTGAtcctc        <  1:4272937/150‑1 (MQ=255)
                                                                                                                                                gCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGCCAGTCTCCATCCTTGGTACGCCAAACCCCATGCTTCGCGCCCTGAtcctcc       >  2:365795/1‑150 (MQ=255)
                                                                                                                                                  tGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGCCAGTCTCCATCCTTGGTACGCCAAACCCCATGCTTCGCGCCCTGATCCTCCTa     >  2:3871305/1‑150 (MQ=255)
                                                                                                                                                  tGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGCCAGTCTCCATCCTTGGTACGCCAAACCCCATGCTTCGCGCCCTGATCCTCCTa     <  1:1061013/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TACGTGGGCAAGCGCAACGTGAAAGTGCGCGAGGGCGGCGTGCTGGCCGACGTGGCGCCAACCATGCTGAAGCTGCTGGGGCTGGAAAAGCCGGTGGAAATGACCGGTACCTCGATTCTGGTCGACGCCTGAATGTCGAAAGGGGCTGCGTTGCAGCCCCAATTGCACATGAATTCCAGACAAACGCCCCCCTGCATCCGCAGCGAGGCGTTTTTTTTGCCTGCCCAGGCGGGCATACTAGGCCAGTCTCCATCCTTGGTACGCCAAACCCCATGCTTCGCGCCCTGATCCTCCTAGCC  >  NC_002947/5767358‑5767656

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