Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,510,659 A→G 5.2% K98K (AAA→AAG 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,6590AG 5.2% 488.2 / 9.7 154K98K (AAA→AAGurtAurea ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (70/75);  new base G (4/4);  total (74/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.28e-01

Only 100 of 158 total aligned reads displayed.

TCGCTCAAGGACATGGCGCTGATGACCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAA  >  NC_002947/5510510‑5510808
                                                                                                                                                     |                                                                                                                                                     
   cTCAAGGACATGGCGCTGATGACCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACACAGTg                                                                                                                                                    >  2:4807774/1‑150 (MQ=255)
   cTCAAGGACATGGCGCTGATGACCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTg                                                                                                                                                    >  2:3451639/1‑150 (MQ=255)
   cTCAAGGACATGGCGCTGATGACCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTg                                                                                                                                                    <  1:695628/150‑1 (MQ=255)
      aaGGACATGGCGCTGATGACCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTggac                                                                                                                                                 >  2:4933492/1‑148 (MQ=255)
      aaGGACATGGCGCTGATGACCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGcc                                                                                                                                                 <  2:1587360/150‑1 (MQ=255)
       aGGACATGGCGCTGATGACCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCg                                                                                                                                                >  1:172501/1‑150 (MQ=255)
        ggACATGGCGCTGATGACCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCgt                                                                                                                                               <  1:3077441/150‑1 (MQ=255)
          aCATGGCGCTGATGACCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCgtgg                                                                                                                                             <  1:1347366/150‑1 (MQ=255)
           cATGGCGCTGATGACCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCgtggt                                                                                                                                            <  2:3815460/150‑1 (MQ=255)
            aTGGCGCTGATGACCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCgtggtg                                                                                                                                           >  2:5083129/1‑150 (MQ=255)
               gcgcTGATGACCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTc                                                                                                                                        <  2:2812089/150‑1 (MQ=255)
                  cTGATGACCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGgc                                                                                                                                     >  2:2945204/1‑150 (MQ=255)
                   tgatgaCCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGgct                                                                                                                                    >  2:709477/1‑150 (MQ=255)
                      tgaCCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGgctgct                                                                                                                                 <  1:2407579/150‑1 (MQ=255)
                      tgaCCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGgctgct                                                                                                                                 <  2:2129117/150‑1 (MQ=255)
                       gaCCATCGAGCAGATAAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGctgctg                                                                                                                                >  2:2892864/1‑150 (MQ=255)
                          cATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGAc                                                                                                                             >  1:3406583/1‑150 (MQ=255)
                           aTCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACa                                                                                                                            <  1:2722541/150‑1 (MQ=255)
                           aTCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACa                                                                                                                            >  2:2313135/1‑150 (MQ=255)
                             cGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATc                                                                                                                          <  1:324773/150‑1 (MQ=255)
                              gAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATcc                                                                                                                         <  2:1072225/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)
                                 cAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACCAAGTGGCCGTGGTGTTCGgctg                                                                                                                                   >  2:4031858/1‑137 (MQ=255)
                                  aGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCgtgt                                                                                                                     <  2:1734871/150‑1 (MQ=255)
                                      cAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTcgcg                                                                                                                 >  1:3568342/1‑150 (MQ=255)
                                       aaCGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTcgc                                                                                                                  >  1:4790061/1‑148 (MQ=255)
                                       aaCGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCgcgc                                                                                                                >  2:3557888/1‑150 (MQ=255)
                                         cGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCaa                                                                                                              <  2:2518827/150‑1 (MQ=255)
                                          gCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAg                                                                                                             >  2:1810689/1‑150 (MQ=255)
                                           ccAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGTGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGt                                                                                                            <  1:1362904/150‑1 (MQ=255)
                                           ccAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGt                                                                                                            <  1:251909/150‑1 (MQ=255)
                                           ccAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGt                                                                                                            <  2:4395490/150‑1 (MQ=255)
                                           ccAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGt                                                                                                            <  2:4395536/150‑1 (MQ=255)
                                            cAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCgtgg                                                                                                                                             <  1:1446243/116‑1 (MQ=255)
                                            cAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCgtgg                                                                                                                                             >  2:1446243/1‑116 (MQ=255)
                                            cAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTc                                                                                                           >  2:240438/1‑150 (MQ=255)
                                             aaGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCt                                                                                                          <  2:637202/150‑1 (MQ=255)
                                                    gCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGc                                                                                                   >  2:2763395/1‑150 (MQ=255)
                                                    gCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGc                                                                                                   >  2:732347/1‑150 (MQ=255)
                                                              cAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCga                                                                                         <  1:5043307/150‑1 (MQ=255)
                                                               aaGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCgaa                                                                                        <  1:254430/150‑1 (MQ=255)
                                                               aaGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCgaa                                                                                        >  1:3190422/1‑150 (MQ=255)
                                                               aaGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCgaa                                                                                        <  1:3925074/150‑1 (MQ=255)
                                                                agagGCTGGAACCGGGGGCGGTCGCCCCGGCGTCCACCTCGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAATGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCgaag                                                                                       <  2:3531798/150‑1 (MQ=255)
                                                                      tGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCa                                                                                 >  2:3436776/1‑150 (MQ=255)
                                                                       ggAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAAcaccg                                                                                <  1:732373/150‑5 (MQ=255)
                                                                       ggAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAAcaccg                                                                                <  1:732347/150‑5 (MQ=255)
                                                                           ccGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGc                                                                            <  1:1714191/150‑1 (MQ=255)
                                                                           ccGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGc                                                                            <  1:344936/150‑1 (MQ=255)
                                                                                ggtggtCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCctgct                                                                       <  2:2605474/150‑1 (MQ=255)
                                                                                   ggtCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGtt                                                                    >  1:2930767/1‑150 (MQ=255)
                                                                                     tCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCTACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGtt                                                                                                       <  2:5204281/113‑1 (MQ=255)
                                                                                       gACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTAc                                                                >  2:591744/1‑150 (MQ=255)
                                                                                       gACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTAc                                                                <  1:2990344/150‑1 (MQ=255)
                                                                                       gACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCAGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTAc                                                                <  1:2539139/150‑1 (MQ=255)
                                                                                           cGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCgg                                                            >  2:3313171/1‑150 (MQ=255)
                                                                                           cGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCgg                                                            >  2:3313176/1‑150 (MQ=255)
                                                                                             gCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCgt                                                                                                                                               <  2:4599057/65‑1 (MQ=255)
                                                                                                tCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGgt                                                                                               <  2:2688440/110‑1 (MQ=255)
                                                                                                  cgattGGCCGCTCTTCGCTGAAAAGGCGCGCCAGCTGATTTCTCAGGATAAGGTCGCAGCCGTGTTCGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTa                                                     <  1:2656348/146‑1 (MQ=255)
                                                                                                  cAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTa                                                     <  2:4284947/150‑1 (MQ=255)
                                                                                                  cAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTa                                                     <  1:252323/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                                                    >  1:3119883/1‑150 (MQ=255)
                                                                                                        gCCGCTCTTCGCTGAAAAGGCGCGCCAGCTGATTTCTCAGGATAAGGTCGCAGCCGTGTTCGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAagg                                               >  1:4133602/1‑150 (MQ=255)
                                                                                                          cGCTCTTCGCTGAAAAGGCGCGCCAGCTGATTTCTCAGGATAAGGTCGCAGCCGTGTTCGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAaggtg                                             >  2:291867/1‑148 (MQ=255)
                                                                                                          cGCTCTTCGCTGAAAAGGCGCGCCAGCTGATTTCTCAGGATAAGGTCGCAGCCGTGTTCGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAaggtg                                             >  1:4971632/1‑148 (MQ=255)
                                                                                                           gCTCTTCGCTGAAAAGGCGCGCCAGCTGATTTCTCAGGATAAGGTCGCAGCCGTGTTCGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAAGGTGa                                            <  2:1111401/150‑1 (MQ=255)
                                                                                                           gCTCTTCGCTGAAAAGGCGCGCCAGCTGATTTCTCAGGATAAGGTCGCAGCCGTGTTCGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAAGGTGa                                            <  1:4429763/150‑1 (MQ=255)
                                                                                                            cTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGaa                                           <  1:240438/150‑1 (MQ=255)
                                                                                                              gTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAaga                                         <  1:2763395/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)
                                                                                                                 cGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGAt                                      <  2:229393/150‑1 (MQ=255)
                                                                                                                 cGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGAt                                      <  2:5179273/150‑1 (MQ=255)
                                                                                                                    tGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTc                                   <  2:3568342/150‑1 (MQ=255)
                                                                                                                      aaaaGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGt                                                      <  2:2472768/129‑1 (MQ=255)
                                                                                                                      aaaaGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGc                                 <  2:279065/150‑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:4316618/1‑150 (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:4777762/150‑1 (MQ=255)
                                                                                                                                       ctgACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTAcaca                <  2:4778140/150‑1 (MQ=255)
                                                                                                                                        tgACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAg               >  1:820805/1‑150 (MQ=255)
                                                                                                                                        tgACCCAGGACAAAGTGGCCGTGGTGGTCGGCTGCTGGACAGCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAg               >  2:2560753/1‑150 (MQ=255)
                                                                                                                                          aCCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGgc             >  2:1234492/1‑150 (MQ=255)
                                                                                                                                            ccAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGgcgc           >  2:1802853/1‑150 (MQ=255)
                                                                                                                                            ccAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGgcgc           >  1:2967894/1‑150 (MQ=255)
                                                                                                                                            ccAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGgcgc           <  2:3716216/150‑1 (MQ=255)
                                                                                                                                                   aaaGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCg    >  2:1088833/1‑150 (MQ=255)
                                                                                                                                                     aGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGaa  >  1:2240546/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGCTCAAGGACATGGCGCTGATGACCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAA  >  NC_002947/5510510‑5510808

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: