Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,510,662 G→C 5.2% V99V (GTG→GTC 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,6620GC 5.2% 510.4 / 13.5 153V99V (GTG→GTCurtAurea ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (70/75);  new base C (4/4);  total (74/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 160 total aligned reads displayed.

CTCAAGGACATGGCGCTGATGACCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCA  >  NC_002947/5510513‑5510811
                                                                                                                                                     |                                                                                                                                                     
   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)
           ggCGCTGATGACCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGtt                                                                                                                                            <  1:859491/150‑1 (MQ=255)
            gcgcTGATGACCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTc                                                                                                                                           <  2:2812089/150‑1 (MQ=255)
            gcgcTGATGACCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCAGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTc                                                                                                                                           <  2:2639592/150‑1 (MQ=255)
               cTGATGACCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGgc                                                                                                                                        >  2:2945204/1‑150 (MQ=255)
                tgatgaCCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGgct                                                                                                                                       >  2:709477/1‑150 (MQ=255)
                 gatgaCCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGgctg                                                                                                                                      >  2:4130937/1‑150 (MQ=255)
                   tgaCCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGgctgct                                                                                                                                    <  2:3603918/150‑1 (MQ=255)
                   tgaCCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGgctgct                                                                                                                                    <  2:2129117/150‑1 (MQ=255)
                   tgaCCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGgctgct                                                                                                                                    <  2:2129080/150‑1 (MQ=255)
                       cATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGAc                                                                                                                                >  1:3406583/1‑150 (MQ=255)
                        aTCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACa                                                                                                                               >  2:2313135/1‑150 (MQ=255)
                           gAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATcc                                                                                                                            <  2:2483879/150‑1 (MQ=255)
                           gAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATcc                                                                                                                            <  2:1072225/150‑1 (MQ=255)
                            aGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCg                                                                                                                           >  2:1714191/1‑150 (MQ=255)
                             gCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCgt                                                                                                                          <  1:3294987/150‑1 (MQ=255)
                              cAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACCAAGTGGCCGTGGTGTTCGgctg                                                                                                                                      >  2:4031858/1‑137 (MQ=255)
                              cAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGgctg                                                                                                                                      <  1:4031858/137‑1 (MQ=255)
                               aGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCgtgt                                                                                                                        <  2:1734871/150‑1 (MQ=255)
                               aGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCgtgt                                                                                                                        <  2:2909572/150‑1 (MQ=255)
                                   cAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTcgcg                                                                                                                    >  1:1668118/1‑150 (MQ=255)
                                    aaCGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTcgc                                                                                                                     <  2:4790061/148‑1 (MQ=255)
                                    aaCGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCgcgc                                                                                                                   >  2:3557888/1‑150 (MQ=255)
                                    aaCGCCAAGGGTGGCGTGAACGGCAAGATGCTCGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCgcgc                                                                                                                   >  2:3557877/1‑150 (MQ=255)
                                      cGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCaa                                                                                                                 <  2:2518827/150‑1 (MQ=255)
                                       gCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAg                                                                                                                >  2:1810689/1‑150 (MQ=255)
                                        ccAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGTGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGt                                                                                                               <  1:1363134/150‑1 (MQ=255)
                                        ccAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGt                                                                                                               <  2:4395536/150‑1 (MQ=255)
                                        ccAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGt                                                                                                               <  2:4395490/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:732347/1‑150 (MQ=255)
                                                 gCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGc                                                                                                      >  2:2763395/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:254430/150‑1 (MQ=255)
                                                            aaGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCgaa                                                                                           >  1:3190422/1‑150 (MQ=255)
                                                            aaGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCgaa                                                                                           <  1:3925074/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)
                                                                     gAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAAc                                                                                  >  1:2605474/1‑150 (MQ=255)
                                                                        ccGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGc                                                                               <  2:3000484/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                                                                   >  2:591744/1‑150 (MQ=255)
                                                                                    gACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTAc                                                                   <  1:2990344/150‑1 (MQ=255)
                                                                                        cGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCgg                                                               >  2:3313176/1‑150 (MQ=255)
                                                                                          gCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCgt                                                                                                                                                  >  1:4599057/1‑65 (MQ=255)
                                                                                          gCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCgt                                                                                                                                                  <  2:4599057/65‑1 (MQ=255)
                                                                                             tCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGgt                                                                                                  >  1:2688440/1‑110 (MQ=255)
                                                                                               cAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTa                                                        <  2:4284947/150‑1 (MQ=255)
                                                                                               cAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTa                                                        <  2:776702/150‑1 (MQ=255)
                                                                                               cAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTa                                                        <  1:252323/150‑1 (MQ=255)
                                                                                               cAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTa                                                        <  1:252363/150‑1 (MQ=255)
                                                                                                   tGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAg                                                    <  2:278135/150‑1 (MQ=255)
                                                                                                     gCCGCTCTTCGCTGAAAAGGCGCGCCAGCTGATTTCTCAGGATAAGGTCGCAGCCGTGTTCGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAagg                                                  >  1:4133602/1‑150 (MQ=255)
                                                                                                       cGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCg                                                >  1:710358/1‑150 (MQ=255)
                                                                                                       cGCTCTTCGCTGAAAAGGCGCGCCAGCTGATTTCTCAGGATAAGGTCGCAGCCGTGTTCGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAaggtg                                                >  1:4971632/1‑148 (MQ=255)
                                                                                                        gCTCTTCGCTGAAAAGGCGCGCCAGCTGATTTCTCAGGATAAGGTCGCAGCCGTGTTCGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAAGGTGa                                               <  2:1111401/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)
                                                                                                              cGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCa                                                                                    >  2:3749085/1‑107 (MQ=255)
                                                                                                              cGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCa                                                                                    <  1:3749085/107‑1 (MQ=255)
                                                                                                              cGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGAt                                         <  2:229393/150‑1 (MQ=255)
                                                                                                              cGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGAt                                         <  2:5179273/150‑1 (MQ=255)
                                                                                                               gCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATg                                        >  1:3583574/1‑150 (MQ=255)
                                                                                                                cTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGt                                       <  2:349119/150‑1 (MQ=255)
                                                                                                                  gAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTCCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCg                                     <  1:1978327/150‑1 (MQ=255)
                                                                                                                   aaaaGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGt                                                         >  1:2472768/1‑129 (MQ=255)
                                                                                                                    aaaGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGcc                                   <  1:3451639/150‑1 (MQ=255)
                                                                                                                               aaCTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCt                        >  1:3586969/1‑150 (MQ=255)
                                                                                                                                aCTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTa                       >  1:3159329/1‑150 (MQ=255)
                                                                                                                                  tgctgACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTCTGAGGGCGAaga                                            >  1:2830197/1‑127 (MQ=255)
                                                                                                                                  tgctgACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGc                                                 >  1:66878/1‑122 (MQ=255)
                                                                                                                                  tgctgACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAaga                                            <  2:2830197/127‑1 (MQ=255)
                                                                                                                                  tgctgACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTaca                     >  2:2090652/1‑150 (MQ=255)
                                                                                                                                    ctgACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCaa                                                                                   <  1:839555/86‑1 (MQ=255)
                                                                                                                                    ctgACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTAcaca                   <  2:4777762/150‑1 (MQ=255)
                                                                                                                                     tgACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAg                  >  2:2560809/1‑150 (MQ=255)
                                                                                                                                     tgACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAg                  >  1:2256943/1‑150 (MQ=255)
                                                                                                                                       aCCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGgc                >  2:1234492/1‑150 (MQ=255)
                                                                                                                                         ccAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGgcgc              >  1:2967894/1‑150 (MQ=255)
                                                                                                                                         ccAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGgcgc              <  2:3716216/150‑1 (MQ=255)
                                                                                                                                         ccAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGgcgc              >  2:1802853/1‑150 (MQ=255)
                                                                                                                                                aaaGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCg       >  2:1088833/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)
                                                                                                                                                     |                                                                                                                                                     
CTCAAGGACATGGCGCTGATGACCATCGAGCAGATCAACGCCAAGGGTGGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCA  >  NC_002947/5510513‑5510811

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: