Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,510,720 G→A 14.1% E119K (GAA→AAA) ‡ urtA → urea ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,510,7200GA14.1% 374.0 / 45.9 149E119K (GAA→AAA) ‡urtAurea ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (69/59);  new base A (11/10);  total (80/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

Only 100 of 159 total aligned reads displayed.

CAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGC  >  NC_002947/5510572‑5510869
                                                                                                                                                    |                                                                                                                                                     
cAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCga                                                                                                                                                      <  1:5043307/150‑1 (MQ=255)
 aaGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCgaa                                                                                                                                                     <  1:3925074/150‑1 (MQ=255)
 aaGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCgaa                                                                                                                                                     >  1:3190422/1‑150 (MQ=255)
 aaGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCgaa                                                                                                                                                     <  1:254430/150‑1 (MQ=255)
  agagGCTGGAACCGGGGGCGGTCGCCCCGGCGTCCACCTCGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAATGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCgaag                                                                                                                                                    <  2:3531798/150‑1 (MQ=255)
   gATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCgaaga                                                                                                                                                   >  1:4386524/1‑150 (MQ=255)
     tGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAAc                                                                                                                                                 >  1:776702/1‑150 (MQ=255)
        tGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCa                                                                                                                                              >  2:3436776/1‑150 (MQ=255)
         ggAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAAcaccg                                                                                                                                             <  1:732373/150‑5 (MQ=255)
             ccGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGc                                                                                                                                         <  2:3000484/150‑1 (MQ=255)
                  ggtggtCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCctgct                                                                                                                                    <  2:2605474/150‑1 (MQ=255)
                     ggtCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGtt                                                                                                                                 >  1:2930767/1‑150 (MQ=255)
                       tCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCt                                                                                                                               >  2:573265/1‑150 (MQ=255)
                         gACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTAc                                                                                                                             >  2:591744/1‑150 (MQ=255)
                             cGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCgg                                                                                                                         >  2:3313171/1‑150 (MQ=255)
                                    cgattGGCCGCTCTTCGCTGAAAAGGCGCGCCAGCTGATTTCTCAGGATAAGGTCGCAGCCGTGTTCGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTa                                                                                                                  <  1:2656348/146‑1 (MQ=255)
                                    cAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTa                                                                                                                  <  1:252363/150‑1 (MQ=255)
                                    cAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTa                                                                                                                  <  1:252323/150‑1 (MQ=255)
                                    cAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTa                                                                                                                  <  2:4284947/150‑1 (MQ=255)
                                    cAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTa                                                                                                                  <  2:4221832/150‑1 (MQ=255)
                                     aaCTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACCTCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGGTCTCCCCGGTGccggat                                                                                                                 >  2:1165634/1‑145 (MQ=255)
                                     aaCTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTCCCCGGTGCAGTag                                                                                                                 >  2:3696854/1‑149 (MQ=255)
                                     aaCTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTAt                                                                                                                 >  1:3119883/1‑150 (MQ=255)
                                       cTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGa                                                                                                               <  2:1763034/150‑1 (MQ=255)
                                        tGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAg                                                                                                              <  2:278135/150‑1 (MQ=255)
                                        tGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAg                                                                                                              >  2:1478350/1‑150 (MQ=255)
                                            cGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCg                                                                                                          >  1:710358/1‑150 (MQ=255)
                                            cGCTCTTCGCTGAAAAGGCGCGCCAGCTGATTTCTCAGGATAAGGTCGCAGCCGTGTTCGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAaggtg                                                                                                          >  2:291867/1‑148 (MQ=255)
                                             gCTCTTCGCTGAAAAGGCGCGCCAGCTGATTTCTCAGGATAAGGTCGCAGCCGTGTTCGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAAGGTGa                                                                                                         <  1:4429763/150‑1 (MQ=255)
                                                gTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAaga                                                                                                      <  1:2763395/150‑1 (MQ=255)
                                                gTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAaga                                                                                                      <  2:1884154/150‑1 (MQ=255)
                                                 ttCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAagag                                                                                                     <  1:4130937/150‑1 (MQ=255)
                                                   cGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCa                                                                                                                                              >  2:3749085/1‑107 (MQ=255)
                                                   cGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGAt                                                                                                   <  2:229393/150‑1 (MQ=255)
                                                    gCTGAAAAGAGCCGTTAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATg                                                                                                  >  2:4438672/1‑150 (MQ=255)
                                                      tGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGTCCGGGGTGATCGGCTGCTGGACAACCGTCTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACACAACGGCCTGCTGATCTACCCGGTGCAGTATGAAGGCGAAGAGATgac                                                                                                >  1:3001662/1‑148 (MQ=255)
                                                       gAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTCCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCg                                                                                               <  1:1978327/150‑1 (MQ=255)
                                                        aaaaGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGt                                                                                                                   >  1:2472768/1‑129 (MQ=255)
                                                        aaaaGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGt                                                                                                                   <  2:2472768/129‑1 (MQ=255)
                                                        aaaaGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGc                                                                                              <  2:279065/150‑1 (MQ=255)
                                                         aaaGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGcc                                                                                             <  1:3451639/150‑1 (MQ=255)
                                                                    aaCTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCt                                                                                  >  1:3586969/1‑150 (MQ=255)
                                                                       tgctgACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGc                                                                                                           <  2:66878/122‑1 (MQ=255)
                                                                       tgctgACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGc                                                                                                           >  1:66878/1‑122 (MQ=255)
                                                                       tgctgACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTaca                                                                               >  2:2090652/1‑150 (MQ=255)
                                                                         ctgACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCaa                                                                                                                                             <  1:839555/86‑1 (MQ=255)
                                                                         ctgACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTAcaca                                                                             <  2:4778140/150‑1 (MQ=255)
                                                                         ctgACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTAcaca                                                                             <  2:4777762/150‑1 (MQ=255)
                                                                          tgACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAg                                                                            >  2:3727896/1‑150 (MQ=255)
                                                                          tgACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAg                                                                            >  2:2560809/1‑150 (MQ=255)
                                                                              ccAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGgcgc                                                                        >  1:2967894/1‑150 (MQ=255)
                                                                                   aCAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGccc                                                                   <  2:4386524/150‑1 (MQ=255)
                                                                                     aaaGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCg                                                                 >  2:1088833/1‑150 (MQ=255)
                                                                                       aGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGaa                                                               >  1:2240546/1‑150 (MQ=255)
                                                                                        gTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAAc                                                              <  2:3607468/150‑1 (MQ=255)
                                                                                         tGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAAcc                                                             >  1:1300329/1‑150 (MQ=255)
                                                                                         tGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTAGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAAcc                                                             >  2:5121353/1‑150 (MQ=255)
                                                                                          ggCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCa                                                            <  2:1025797/150‑1 (MQ=255)
                                                                                             cGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCa                                                         >  1:3433527/1‑150 (MQ=255)
                                                                                             agccgtgtTCGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAAGGTGAGGAATCCGAGCGCAACGTTTTCTATACGGGCGCGGCACCGAACCAGCa                                                         <  1:1709933/146‑1 (MQ=255)
                                                                                              gtggtgTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCaa                                                        <  2:3568303/150‑1 (MQ=255)
                                                                                               ccgtgtTCGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAAGGTGAGGAATCCGAGCGCAACGTTTTCTATACGGGCGCGGCACCGAACCAGCagg                                                       >  2:5241592/3‑148 (MQ=255)
                                                                                                cgtgtTCGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAAGGTGAGGAATCCGAGCGCAACGTTTTCTATACGGGCGCGGCACCGAACCAGCAGGc                                                      >  1:4157793/2‑150 (MQ=255)
                                                                                                     tCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATcc                                                 <  1:5193582/150‑1 (MQ=255)
                                                                                                       ggCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTAccc                                                                                                                           <  2:1227765/74‑1 (MQ=255)
                                                                                                       ggCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTAccc                                                                                                                           >  1:1227765/1‑74 (MQ=255)
                                                                                                          tgctgGACATCCGTGTCGCGCTAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGcc                                            <  1:2892864/150‑1 (MQ=255)
                                                                                                          tgctgGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGcc                                            >  2:2208667/1‑150 (MQ=255)
                                                                                                             tgGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTg                                         >  2:1921545/1‑150 (MQ=255)
                                                                                                               gACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAAGGTGAGGAATCCGAGCGCAACGTTTTCTATACGGGCGCGGCACCGAACCAGCAGGCAATCCCTGCTGTTGa                                       <  1:1636560/150‑1 (MQ=255)
                                                                                                               gACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGa                                       >  1:2383576/1‑150 (MQ=255)
                                                                                                               gACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGa                                       >  1:5442699/1‑150 (MQ=255)
                                                                                                               gACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGa                                       >  1:5422223/1‑150 (MQ=255)
                                                                                                                  aTCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGTCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGACGCCCCGAACCAGCAAGCCATCCCGGCCTTGGAATa                                    >  1:786003/1‑150 (MQ=255)
                                                                                                                    cggtgtCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAAGGTGAGGAATCCGAGCGCAACGTTTTCTATACGGGCGCGGCACCGAACCAGCAGGCAATCCCTGCTGTTGACTatc                                  >  1:1093673/3‑148 (MQ=255)
                                                                                                                    cggtgtCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAAGGTGAGGAATCCGAGCGCAACGTTTTCTATACGGGCGCGGCACCGAACCAGCAGGCAATCCCTGCTGTTGACTatc                                  <  1:793379/148‑3 (MQ=255)
                                                                                                                    ccGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATAcc                                  >  1:702789/1‑150 (MQ=255)
                                                                                                                    ccGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATAcc                                  >  1:2448142/1‑150 (MQ=255)
                                                                                                                     cGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCt                                 <  1:1810689/150‑1 (MQ=255)
                                                                                                                       tgtCGCGCAAGTCTGTTATGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTAccc                                                                                                                           <  2:2663956/58‑1 (MQ=255)
                                                                                                                       tgtCGCGCAAGTCTGTTATGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTAccc                                                                                                                           >  1:2663956/1‑58 (MQ=255)
                                                                                                                        gtCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCtgat                              <  1:4450277/150‑1 (MQ=255)
                                                                                                                         tCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCtgatg                             >  1:3567681/1‑150 (MQ=255)
                                                                                                                         tCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCtgatg                             >  1:3567800/1‑150 (MQ=255)
                                                                                                                         tCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCtgatg                             >  1:1081113/1‑150 (MQ=255)
                                                                                                                           gcgcAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGACTACCTGATGAg                           >  2:3932390/1‑150 (MQ=255)
                                                                                                                            cgcAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGc                          <  1:4807774/150‑1 (MQ=255)
                                                                                                                              cAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGa                        <  1:2208667/150‑1 (MQ=255)
                                                                                                                              cAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGACTGCTGTTCTACCCGGTGCAGTATGAGTTCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGa                        <  1:4933492/150‑1 (MQ=255)
                                                                                                                               aaGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAg                       >  1:4122523/1‑150 (MQ=255)
                                                                                                                                  tCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGAc                    >  2:3578625/1‑150 (MQ=255)
                                                                                                                                         tGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGAcggatgtc             <  2:2383576/150‑6 (MQ=255)
                                                                                                                                         tGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCa             >  1:36115/1‑150 (MQ=255)
                                                                                                                                           ccGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAgc           <  2:3406583/150‑1 (MQ=255)
                                                                                                                                                 ttCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAg     >  1:4994043/1‑150 (MQ=255)
                                                                                                                                                 ttCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAg     >  2:794008/1‑150 (MQ=255)
                                                                                                                                                   cGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCa       <  1:1581292/146‑1 (MQ=255)
                                                                                                                                                   cGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCa       >  2:1581292/1‑146 (MQ=255)
                                                                                                                                                    gaagaaCTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAgcgc  >  1:2804754/1‑150 (MQ=255)
                                                                                                                                                    gaagaaCTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAgcgc  >  1:4314188/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGC  >  NC_002947/5510572‑5510869

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: