Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,510,677 C→T 8.3% G104G (GGC→GGT 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,6770CT 8.3% 424.7 / 19.9 144G104G (GGC→GGTurtAurea ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (65/67);  new base T (6/6);  total (71/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 154 total aligned reads displayed.

CTGATGACCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCG  >  NC_002947/5510528‑5510824
                                                                                                                                                     |                                                                                                                                                   
cTGATGACCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGgc                                                                                                                                                     >  2:2945204/1‑150 (MQ=255)
 tgatgaCCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGgct                                                                                                                                                    >  2:709477/1‑150 (MQ=255)
  gatgaCCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGgctg                                                                                                                                                   >  2:4130937/1‑150 (MQ=255)
    tgaCCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGgctgct                                                                                                                                                 <  2:2129080/150‑1 (MQ=255)
    tgaCCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGgctgct                                                                                                                                                 <  1:2407579/150‑1 (MQ=255)
    tgaCCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGgctgct                                                                                                                                                 <  2:3603918/150‑1 (MQ=255)
     gaCCATCGAGCAGATAAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGctgctg                                                                                                                                                >  2:2892864/1‑150 (MQ=255)
         aTCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACa                                                                                                                                            >  2:2313135/1‑150 (MQ=255)
           cGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATc                                                                                                                                          <  1:324773/150‑1 (MQ=255)
            gAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATcc                                                                                                                                         >  2:2106457/1‑150 (MQ=255)
             aGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCg                                                                                                                                        >  2:1714191/1‑150 (MQ=255)
             aGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCg                                                                                                                                        >  1:3607468/1‑150 (MQ=255)
              gCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCgt                                                                                                                                       <  1:3294987/150‑1 (MQ=255)
                   tCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTcgc                                                                                                                                  <  2:2131470/150‑1 (MQ=255)
                    cAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTcgcg                                                                                                                                 >  1:3568342/1‑150 (MQ=255)
                     aaCGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCgcgc                                                                                                                                >  2:3557888/1‑150 (MQ=255)
                     aaCGCCAAGGGTGGCGTGAACGGCAAGATGCTCGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCgcgc                                                                                                                                >  2:3557877/1‑150 (MQ=255)
                       cGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCaa                                                                                                                              <  2:2518827/150‑1 (MQ=255)
                         ccAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGTGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGt                                                                                                                            <  1:1363134/150‑1 (MQ=255)
                         ccAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGt                                                                                                                            <  2:4395490/150‑1 (MQ=255)
                           aaGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCt                                                                                                                          <  2:637202/150‑1 (MQ=255)
                                  gCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGc                                                                                                                   >  2:2763395/1‑150 (MQ=255)
                                  gCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGc                                                                                                                   >  2:732347/1‑150 (MQ=255)
                                    gTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCg                                                                                                                 >  1:4195718/1‑150 (MQ=255)
                                        aCGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGgtgt                                                                                                             <  1:4980670/150‑1 (MQ=255)
                                          ggCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTc                                                                                                           <  2:4960444/150‑1 (MQ=255)
                                             aaGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCgaa                                                                                                        <  1:3925074/150‑1 (MQ=255)
                                             aaGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCgaa                                                                                                        <  1:254430/150‑1 (MQ=255)
                                                 tGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAAc                                                                                                    >  1:776702/1‑150 (MQ=255)
                                                     ggAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAAcaccg                                                                                                <  1:732373/150‑5 (MQ=255)
                                                         ccGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGc                                                                                            <  1:344936/150‑1 (MQ=255)
                                                         ccGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGc                                                                                            <  1:1714191/150‑1 (MQ=255)
                                                              ggtggtCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCctgct                                                                                       <  2:2605474/150‑1 (MQ=255)
                                                                   tCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCTACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGtt                                                                                                                       <  2:5204281/113‑1 (MQ=255)
                                                                   tCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCGACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGtt                                                                                                                       >  1:5204281/1‑113 (MQ=255)
                                                                     gACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTAc                                                                                <  1:2990344/150‑1 (MQ=255)
                                                                     gACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTAc                                                                                >  2:591744/1‑150 (MQ=255)
                                                                         cGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCgg                                                                            >  2:3313171/1‑150 (MQ=255)
                                                                         cGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCgg                                                                            >  2:3313176/1‑150 (MQ=255)
                                                                              tCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGgt                                                                                                               >  1:2688440/1‑110 (MQ=255)
                                                                                cgattGGCCGCTCTTCGCTGAAAAGGCGCGCCAGCTGATTTCTCAGGATAAGGTCGCAGCCGTGTTCGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTa                                                                     <  1:2656348/146‑1 (MQ=255)
                                                                                cAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTa                                                                     <  1:252323/150‑1 (MQ=255)
                                                                                cAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTa                                                                     <  1:252363/150‑1 (MQ=255)
                                                                                cAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTa                                                                     <  2:4221832/150‑1 (MQ=255)
                                                                                cAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTa                                                                     <  2:776702/150‑1 (MQ=255)
                                                                                 aaCTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTCCCCGGTGCAGTag                                                                    >  2:3696854/1‑149 (MQ=255)
                                                                                 aaCTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTAt                                                                    >  2:2343329/1‑150 (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)
                                                                                      gCCGCTCTTCGCTGAAAAGGCGCGCCAGCTGATTTCTCAGGATAAGGTCGCAGCCGTGTTCGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAagg                                                               >  1:4133602/1‑150 (MQ=255)
                                                                                        cGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCg                                                             >  1:710358/1‑150 (MQ=255)
                                                                                        cGCTCTTCGCTGAAAAGGCGCGCCAGCTGATTTCTCAGGATAAGGTCGCAGCCGTGTTCGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAaggtg                                                             >  2:291867/1‑148 (MQ=255)
                                                                                        cGCTCTTCGCTGAAAAGGCGCGCCAGCTGATTTCTCAGGATAAGGTCGCAGCCGTGTTCGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAaggtg                                                             >  1:4971632/1‑148 (MQ=255)
                                                                                         gCTCTTCGCTGAAAAGGCGCGCCAGCTGATTTCTCAGGATAAGGTCGCAGCCGTGTTCGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAAGGTGa                                                            <  1:4429763/150‑1 (MQ=255)
                                                                                         gCTCTTCGCTGAAAAGGCGCGCCAGCTGATTTCTCAGGATAAGGTCGCAGCCGTGTTCGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAAGGTGa                                                            <  2:3851963/150‑1 (MQ=255)
                                                                                         gCTCTTCGCTGAAAAGGCGCGCCAGCTGATTTCTCAGGATAAGGTCGCAGCCGTGTTCGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAAGGTGa                                                            <  2:1111401/150‑1 (MQ=255)
                                                                                            gTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAaga                                                         <  1:2763395/150‑1 (MQ=255)
                                                                                            gTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAaga                                                         <  2:1884154/150‑1 (MQ=255)
                                                                                             ttCGCTGAAAAGGCGCGCCAGCTGATTTCTCAGGATAAGGTCGCAGCCGTGTTCGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAAGGTGAggaa                                                        >  1:4420799/1‑149 (MQ=255)
                                                                                             ttCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAagag                                                        <  1:4130937/150‑1 (MQ=255)
                                                                                               cGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCa                                                                                                 >  2:3749085/1‑107 (MQ=255)
                                                                                                gCTGAAAAGAGCCGTTAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATg                                                     >  2:4438672/1‑150 (MQ=255)
                                                                                                 cTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGt                                                    <  2:349119/150‑1 (MQ=255)
                                                                                                  tGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTc                                                   <  2:3568342/150‑1 (MQ=255)
                                                                                                    aaaaGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGt                                                                      >  1:2472768/1‑129 (MQ=255)
                                                                                                    aaaaGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGt                                                                      <  2:2472768/129‑1 (MQ=255)
                                                                                                         agCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAAc                                            >  1:1520418/1‑150 (MQ=255)
                                                                                                           ccGACAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACgt                                          >  1:3979186/1‑150 (MQ=255)
                                                                                                                aaCTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCt                                     >  1:3586969/1‑150 (MQ=255)
                                                                                                                 aCTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTa                                    >  1:3159329/1‑150 (MQ=255)
                                                                                                                 aCTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTa                                    >  1:4316618/1‑150 (MQ=255)
                                                                                                                   tgctgACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTCTGAGGGCGAaga                                                         >  1:2830197/1‑127 (MQ=255)
                                                                                                                   tgctgACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAaga                                                         <  2:2830197/127‑1 (MQ=255)
                                                                                                                     ctgACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCaa                                                                                                <  1:839555/86‑1 (MQ=255)
                                                                                                                     ctgACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTAcaca                                <  2:4777762/150‑1 (MQ=255)
                                                                                                                     ctgACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTAcaca                                <  2:4778140/150‑1 (MQ=255)
                                                                                                                      tgACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAg                               >  1:2256943/1‑150 (MQ=255)
                                                                                                                      tgACCCAGGACAAAGTGGCCGTGGTGGTCGGCTGCTGGACAGCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAg                               >  2:2560753/1‑150 (MQ=255)
                                                                                                                          ccAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGgcgc                           <  2:3716216/150‑1 (MQ=255)
                                                                                                                          ccAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGgcgc                           >  2:1802853/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)
                                                                                                                                     tGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAAcc                >  1:1300329/1‑150 (MQ=255)
                                                                                                                                     tGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTAGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAAcc                >  2:5121353/1‑150 (MQ=255)
                                                                                                                                      ggCCGTGGTTTTCGGCTGCTGGACATCCGTGTCGAGCAAGTCTGTTCTGAAGGTGTTAGAAGAACTAAACTGAATGCTGTTCTACAAGGTGCAGTATGAGGTCGAATAGATGTCGCCTAACGTGTTCTACAAAGGAGACGCCCCTAAACa               >  2:2630234/1‑150 (MQ=255)
                                                                                                                                      ggCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCa               <  2:1025797/150‑1 (MQ=255)
                                                                                                                                         cGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCa            >  1:3433527/1‑150 (MQ=255)
                                                                                                                                         cGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCa            >  2:2079276/1‑150 (MQ=255)
                                                                                                                                         agccgtgtTCGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAAGGTGAGGAATCCGAGCGCAACGTTTTCTATACGGGCGCGGCACCGAACCAGCa            <  1:1709933/146‑1 (MQ=255)
                                                                                                                                          gtggtgTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCaa           <  1:5083129/150‑1 (MQ=255)
                                                                                                                                           ccgtgtTCGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAAGGTGAGGAATCCGAGCGCAACGTTTTCTATACGGGCGCGGCACCGAACCAGCagg          >  2:5241592/3‑148 (MQ=255)
                                                                                                                                                 tCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATcc    <  1:591744/150‑1 (MQ=255)
                                                                                                                                                  cGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAAGGTGAGGAATCCGAGCGCAACGTTTTCTATACGGGCGCGGCACCGAACCAGCAGGCAATccc   <  2:2283227/150‑1 (MQ=255)
                                                                                                                                                   ggCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTAccc                                                                              <  2:1227765/74‑1 (MQ=255)
                                                                                                                                                   ggCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTAccc                                                                              >  1:1227765/1‑74 (MQ=255)
                                                                                                                                                   ggCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCg  <  1:709477/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CTGATGACCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCG  >  NC_002947/5510528‑5510824

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: