Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,798,803 A→T 6.3% intergenic (+22/‑110) PP_t43 → / → PP_2449 tRNA‑Val/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,798,8030AT 6.3% 520.0 / 17.6 159intergenic (+22/‑110)PP_t43/PP_2449tRNA‑Val/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (80/68);  new base T (7/3);  total (88/71)
Fisher's exact test for biased strand distribution p-value = 5.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

Only 100 of 162 total aligned reads displayed.

AAAAGAGCTTTACAGGTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCGCCC  >  NC_002947/2798654‑2798949
                                                                                                                                                     |                                                                                                                                                  
 aaaGAGCTTTACAGGTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGaaa                                                                                                                                                   <  2:533188/150‑1 (MQ=255)
      gCTTTACAGGTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATc                                                                                                                                              <  2:3673833/150‑1 (MQ=255)
       cTTTACAGGTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATc                                                                                                                                              <  1:4680694/149‑1 (MQ=255)
       cTTTACAGGTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATc                                                                                                                                              >  2:4680694/1‑149 (MQ=255)
       cTTTACAGGTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCg                                                                                                                                             >  1:256545/1‑150 (MQ=255)
          tACAGGTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGccc                                                                                                                                          <  1:3349788/150‑1 (MQ=255)
          tACAGGTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGGTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGccc                                                                                                                                          >  2:1021195/1‑150 (MQ=255)
             aGGTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCATTTCTGATCGCCCtt                                                                                                                                        >  2:5190390/1‑149 (MQ=255)
               gTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATc                                                                                                                                              >  2:2495543/1‑141 (MQ=255)
               gTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCttttt                                                                                                                                     <  1:5012893/150‑1 (MQ=255)
               gTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCatttt                                                                                                                                     <  2:5169010/150‑1 (MQ=255)
                     aTCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCAttt                                                                                                                               <  1:2173045/150‑1 (MQ=255)
                     aTCGGGCCACTAATATGCGCACCACACGAAAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATc                                                                                                                                              <  2:5264823/135‑1 (MQ=255)
                     aTCGGGCCACTAATATGCGCACCACACGAAAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATc                                                                                                                                              >  1:5264823/1‑135 (MQ=255)
                       cGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTgc                                                                                                                             <  1:30924/150‑1 (MQ=255)
                       cGGGCCACTAATAGGCGCCCCACACGACAGGCACAGAGCTCAGTTGGTGAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTgc                                                                                                                             <  1:4907135/150‑1 (MQ=255)
                           ccACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGccc                                                                                                                         <  1:1651801/150‑1 (MQ=255)
                              cTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGaa                                                                                                                      >  2:2422381/1‑150 (MQ=255)
                                 atatGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAagca                                                                                                                   <  1:1340717/150‑1 (MQ=255)
                                     gcgcACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCgt                                                                                                               <  1:171863/150‑1 (MQ=255)
                                      cgcACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCCTTTGCGCCCGAAGCAGCgtg                                                                                                              >  1:2973103/1‑150 (MQ=255)
                                      cgcACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCgtg                                                                                                              <  2:369900/150‑1 (MQ=255)
                                          ccacACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAAc                                                                                                          >  1:2567579/1‑150 (MQ=255)
                                          ccacACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAAc                                                                                                          >  2:3888395/1‑150 (MQ=255)
                                             cacGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGtt                                                                                                       <  2:3300442/150‑1 (MQ=255)
                                               cGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTg                                                                                                     <  1:3810395/150‑1 (MQ=255)
                                               cGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTGTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTg                                                                                                     >  2:5380860/1‑150 (MQ=255)
                                                gACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTgg                                                                                                    >  1:1978171/1‑150 (MQ=255)
                                                  cAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTc                                                                                                  <  1:3875364/150‑1 (MQ=255)
                                                     gCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTcgcg                                                                                               <  2:4014183/150‑1 (MQ=255)
                                                       acaTAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGt                                                                                             >  1:1460594/1‑150 (MQ=255)
                                                         aTAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTcc                                                                                           >  2:4271347/1‑150 (MQ=255)
                                                           aGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCgg                                                                                         <  2:2679576/150‑1 (MQ=255)
                                                           aGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCgg                                                                                         <  1:5006097/150‑1 (MQ=255)
                                                           aGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCgg                                                                                         >  2:3598980/1‑150 (MQ=255)
                                                            gCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCggg                                                                                        <  2:684871/150‑1 (MQ=255)
                                                            gCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTATTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCggg                                                                                        <  2:1905317/150‑1 (MQ=255)
                                                             cTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGc                                                                                       >  2:133258/1‑150 (MQ=255)
                                                                aGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATa                                                                                    <  1:2515115/150‑1 (MQ=255)
                                                                  ttGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATAcc                                                                                  >  1:1151487/1‑150 (MQ=255)
                                                                    ggTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATAcccc                                                                                >  1:4549785/1‑150 (MQ=255)
                                                                     gTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCg                                                                               >  2:5021335/1‑150 (MQ=255)
                                                                         gagCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGt                                                                           <  2:20404/150‑1 (MQ=255)
                                                                         gagCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGt                                                                           <  2:2296759/150‑1 (MQ=255)
                                                                         gagCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAAATGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCGTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCgg                                                                           >  1:3590812/1‑149 (MQ=255)
                                                                           gCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAag                                                                                                                     <  2:5033783/106‑1 (MQ=255)
                                                                             accaccTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTa                                                                       <  1:4554762/150‑1 (MQ=255)
                                                                             accaccTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTa                                                                       <  2:1486677/150‑1 (MQ=255)
                                                                             accaccTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTa                                                                       <  1:59000/150‑1 (MQ=255)
                                                                                  cTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATc                                                                  >  1:1475264/1‑150 (MQ=255)
                                                                                   ttGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATcc                                                                 <  2:1927173/150‑1 (MQ=255)
                                                                                     gACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGc                                                               <  1:1211502/150‑1 (MQ=255)
                                                                                     gACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGc                                                               >  2:1287419/1‑150 (MQ=255)
                                                                                       cATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGccc                                                             >  1:3721268/1‑150 (MQ=255)
                                                                                        aTGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTTGTCGCGGTCCGGGCATAACCCGGCGTGCTAGAATCCGcccc                                                            >  1:1374/1‑150 (MQ=255)
                                                                                        aTGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGcccc                                                            >  1:3423735/1‑150 (MQ=255)
                                                                                        aTGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCAATCAGAAATGATCGCCCTGTTTGCATTTGCGCCCGAAGCAGCGTGTAACGGTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGcccc                                                            >  2:1579379/1‑150 (MQ=255)
                                                                                               gggTCGTTGGTTCGAGTCCAATTTTGCCAACCAAAAACAAAAATGGCGATCAGAAAATAACGCCCATTTTGCAATTTCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGCCGTGCTAGAATACGccccaaccaaa                                                     >  1:2889341/1‑143 (MQ=255)
                                                                                               gggTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAAt                                                     >  1:1160552/1‑150 (MQ=255)
                                                                                                     ttGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAg                                               <  1:1268665/150‑1 (MQ=255)
                                                                                                        gtccGAGTCCAAGTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTa                                            <  1:4329803/147‑1 (MQ=255)
                                                                                                        gTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTa                                            >  1:4487472/1‑150 (MQ=255)
                                                                                                           cGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTacac                                         <  1:1817693/150‑1 (MQ=255)
                                                                                                            gAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACAcg                                        >  1:1125574/1‑150 (MQ=255)
                                                                                                                  aaTTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTgcg                                  <  1:5143444/150‑1 (MQ=255)
                                                                                                                   aTTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTgcgc                                 <  2:2082601/150‑1 (MQ=255)
                                                                                                                    ttGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCa                                >  1:4858015/1‑150 (MQ=255)
                                                                                                                          cTACCAAATCGAAAAAGGGTGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAAACCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTgg                          <  1:411540/150‑1 (MQ=255)
                                                                                                                           tACCAAATCGAAAAGGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGc                         >  1:5346972/1‑150 (MQ=255)
                                                                                                                           tACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGc                         >  2:2446643/1‑150 (MQ=255)
                                                                                                                           tACCAAATCGAAAAAGGGCGATCAGAAATGATAGCCCTTTTTGCATTTGCTACCGAAGCAGAGTGTAACGTTTTGTCGCGGTACGGGCATACACCGGCGTGCTAGAATCCGAACCTTCGAATATGGAGCAACACGCGAGCGCAAACTgga                         >  2:564700/1‑149 (MQ=255)
                                                                                                                             gcAAATGCAAAAAGGGCGATCATTTCTGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAAt                                                     <  1:269154/119‑1 (MQ=37)
                                                                                                                             gcAAATGCAAAAAGGGCGATCATTTCTGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAAt                                                     >  2:269154/2‑120 (MQ=37)
                                                                                                                             ccAAATCGAAAAAGGGTGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTgg                          <  1:3140388/147‑1 (MQ=255)
                                                                                                                             ccAAATCGAAAAAGGGTGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTgg                          >  2:3140388/1‑147 (MQ=255)
                                                                                                                             ccAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGaa                                                                    >  2:4158522/1‑105 (MQ=255)
                                                                                                                             ccAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGaa                                                                    <  1:4158522/105‑1 (MQ=255)
                                                                                                                              cAAATCGAAACAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCtgt                      <  2:2226818/150‑1 (MQ=255)
                                                                                                                                 atgcaaaaaGGGCGATCATTTCTGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGt                   >  2:356849/5‑150 (MQ=39)
                                                                                                                                    gAAAAAGGGCGATCAGAGATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAAt                >  1:4381930/1‑150 (MQ=255)
                                                                                                                                    gAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAAt                >  1:2385540/1‑150 (MQ=255)
                                                                                                                                     aaaaaGGGCGATCATTTCTGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATg               >  2:739693/1‑150 (MQ=255)
                                                                                                                                     aaaaaGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATg               >  1:5080825/1‑150 (MQ=255)
                                                                                                                                        aaGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCa            <  2:4676889/150‑1 (MQ=255)
                                                                                                                                         aGGGCGATCATTTCTGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAg           >  1:691310/1‑150 (MQ=255)
                                                                                                                                         aGGGCGATCATTTCTGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAg           >  1:690626/1‑150 (MQ=255)
                                                                                                                                           ggCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGa                                                <  2:2493477/111‑1 (MQ=255)
                                                                                                                                           ggCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGa                                                >  1:2493436/1‑111 (MQ=255)
                                                                                                                                           ggCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGa                                                >  1:2493477/1‑111 (MQ=255)
                                                                                                                                           ggCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTg         >  1:3531892/1‑150 (MQ=255)
                                                                                                                                           ggCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTg         >  2:3425091/1‑150 (MQ=255)
                                                                                                                                            gCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGc        >  1:3608232/1‑150 (MQ=255)
                                                                                                                                            gCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGc        >  1:512143/1‑150 (MQ=255)
                                                                                                                                              gATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTCCACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTc      >  1:778528/1‑150 (MQ=255)
                                                                                                                                              gATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTc      >  1:937046/1‑150 (MQ=255)
                                                                                                                                               aTCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCg     <  1:3907058/150‑1 (MQ=255)
                                                                                                                                                tCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCGc    >  2:4987574/1‑150 (MQ=255)
                                                                                                                                                tCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCGc    >  2:1626097/1‑150 (MQ=255)
                                                                                                                                                acAGAAATGATCGCCCTATTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCGc    <  2:1011910/149‑1 (MQ=255)
                                                                                                                                                  aGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCGccc  <  2:795290/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AAAAGAGCTTTACAGGTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCGCCC  >  NC_002947/2798654‑2798949

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: