Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,196,858 C→A 14.6% intergenic (+30/+23) PP_4574 → / ← PP_4575 hypothetical protein/urea amidolyase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,196,8580CA14.6% 366.6 / 48.5 138intergenic (+30/+23)PP_4574/PP_4575hypothetical protein/urea amidolyase‑like protein
Reads supporting (aligned to +/- strand):  ref base C (63/54);  new base A (11/9);  total (75/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

Only 100 of 140 total aligned reads displayed.

CCCAAGCAGGCCTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGCCAGCGCCAGCGGCGTCA  >  NC_002947/5196709‑5197002
                                                                                                                                                     |                                                                                                                                                 
cccAAGCAGGCCTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTc                                                                                                                                                   >  1:2187873/1‑150 (MQ=255)
 ccAAGCAGGCCTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCg                                                                                                                                                  >  1:655533/1‑150 (MQ=255)
  cAAGCAGGCCTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGc                                                                                                                                                 <  1:4326522/150‑1 (MQ=255)
  cAAGCAGGCCTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGc                                                                                                                                                 <  1:3137273/150‑1 (MQ=255)
      cAGGCCTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCcgc                                                                                                                                             <  1:2678375/150‑1 (MQ=255)
      cAGGCCTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCGAGCCcgc                                                                                                                                             <  2:696714/150‑1 (MQ=255)
          ccTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGaaa                                                                                                                                         <  1:1200224/150‑1 (MQ=255)
          ccTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGaaa                                                                                                                                         <  2:3546286/150‑1 (MQ=255)
          ccTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGaaa                                                                                                                                         <  2:3052403/150‑1 (MQ=255)
          ccTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCGAGCCCGCGaaa                                                                                                                                         <  2:3300787/150‑1 (MQ=255)
          ccTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCGAGCCCGCGaaa                                                                                                                                         <  1:2623712/150‑1 (MQ=255)
          ccTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGCGAGGTTGCCTGTACCGGCCCTTTCGCGGGCGAGCCCGCGaaa                                                                                                                                         <  1:646047/150‑1 (MQ=255)
          ccTGGGAGTACACCGCCTCGACCCATCTGGGGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGaaa                                                                                                                                         <  2:5302689/150‑1 (MQ=255)
           cTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAg                                                                                                                                        >  2:5116112/1‑150 (MQ=255)
           cTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAg                                                                                                                                        >  1:3175540/1‑150 (MQ=255)
           cTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCGAGCCCGCGAAAg                                                                                                                                        <  2:3866447/150‑1 (MQ=255)
            tGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAgg                                                                                                                                       <  1:2146849/150‑1 (MQ=255)
            tGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCGAGCCCGCGAAAgg                                                                                                                                       <  1:4072040/150‑1 (MQ=255)
            tGGGAGTACACCGCCACGACCCATCTGGTGCGGGTCCTGATCGAGGAGAAGATGTATCGGCGGACGCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGGTGCCTGGACCGGCGCTTGCTCGGGCTCGCCCGCGCAAgg                                                                                                                                       <  2:4027385/150‑1 (MQ=255)
                   acacCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCcagca                                                                                                                                >  2:4476746/1‑150 (MQ=255)
                   acacCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGAAATGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGTTGCGCTACCTGTCCATGTGAGGTTGCCTTTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCcagca                                                                                                                                <  1:1422308/150‑1 (MQ=255)
                     acCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCccgccgg                                                                                                                              >  2:1837520/1‑147 (MQ=255)
                     acCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGt                                                                                                                              >  2:1561832/1‑150 (MQ=255)
                     acCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCACTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGt                                                                                                                              >  1:2585560/1‑150 (MQ=255)
                       cGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGCAAGGGCCAGCAGTCt                                                                                                                            >  1:2872399/1‑150 (MQ=255)
                           tCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACg                                                                                                                        >  1:1262948/1‑150 (MQ=255)
                           tCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCGAGCCCGCGAAAg                                                                                                                                        <  1:3334388/134‑1 (MQ=255)
                            cGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGTGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGa                                                                                                                       <  1:4335695/150‑1 (MQ=255)
                                   tCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGGCGTACGACCAACGc                                                                                                                >  1:1869575/1‑150 (MQ=255)
                                   tCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTAGCGGGCTCGCCCGCGAAAGGGCCAGCACTCTTAAGACCAACGc                                                                                                                >  1:1072712/1‑150 (MQ=255)
                                          cgcgACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTTCCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGt                                                                                                         <  2:4714743/150‑1 (MQ=255)
                                            cgACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCCGCAGTCTTACGACCAACCCCCCAGGTaa                                                                                                       >  1:985515/1‑150 (MQ=255)
                                               ccTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCCGCAGTCGTACGACCAACGCCCCAGGTAAAgc                                                                                                    >  1:2417195/1‑150 (MQ=255)
                                                 tGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAgctg                                                                                                  >  1:5380579/1‑150 (MQ=255)
                                                 tGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCCGGCTCGCCCGCGTAAGGGCCAGCAGTCTTAAGACCAAAGCCCCACGTaaacatg                                                                                                  >  1:245938/1‑146 (MQ=255)
                                                  gATCGAGGACACGATGGATCCGCTGACCCTGCATGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAgctgc                                                                                                 <  2:3667215/150‑1 (MQ=255)
                                                  gATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAgctgc                                                                                                 <  2:3922178/150‑1 (MQ=255)
                                                  gATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAgctgc                                                                                                 <  2:1661711/150‑1 (MQ=255)
                                                   aTCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAgctgct                                                                                                >  1:4039745/1‑150 (MQ=255)
                                                    tCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAgctgctg                                                                                               >  1:1308578/1‑150 (MQ=255)
                                                     cGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTCAAgctgcagc                                                                                              >  1:4966225/1‑150 (MQ=255)
                                                        ggACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGctct                                                                                           >  1:3020824/1‑150 (MQ=255)
                                                         gACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCtctc                                                                                          >  2:2775066/1‑150 (MQ=255)
                                                         gACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCtctc                                                                                          >  2:3685093/1‑150 (MQ=255)
                                                         gACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCtctc                                                                                          >  1:4937619/1‑150 (MQ=255)
                                                         gACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCtctc                                                                                          <  1:3624598/150‑1 (MQ=255)
                                                         gACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCGCGCCCGCGGAAGGGCCAGCCGTCGTAAGACCAACGCCCCCGGTAAAGCTGCTGGtctc                                                                                          >  1:1939658/1‑150 (MQ=255)
                                                           cAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGc                                                                                        <  2:3880070/150‑1 (MQ=255)
                                                            aaGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGcc                                                                                       >  2:4236809/1‑150 (MQ=255)
                                                              gATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAg                                                                                     <  2:544484/150‑1 (MQ=255)
                                                               aTGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAgg                                                                                    >  1:3077782/1‑150 (MQ=255)
                                                                   aTCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTc                                                                                <  1:3770643/150‑1 (MQ=255)
                                                                    tCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCGAGCCCGCGAAAgg                                                                                                                                       >  1:1739684/1‑94 (MQ=255)
                                                                    tCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCGAGCCCGCGAAAgg                                                                                                                                       <  2:1739684/94‑1 (MQ=255)
                                                                        cTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCt                                                                           <  1:3883365/150‑1 (MQ=255)
                                                                        cTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCt                                                                           <  2:3591339/150‑1 (MQ=255)
                                                                                gCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTGGCGGGCTCGCCCGCGGAAGGGCCAGCACTGTTAAGAACAACGCCCCAGGTAAAGCAGCTGCGATCGCCCGGCCGCGTCCGGCaccccc                                                                  >  1:1380488/1‑150 (MQ=255)
                                                                                gCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTGGCCCGCGAAAGGGGCAGCAGTGTTACGAGCAACGCCCCAGGTAAAGCTGGTTCTCTCGCCAGGCCTCGTGCTGCaccacc                                                                  >  1:2974767/1‑150 (MQ=255)
                                                                                gCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGCAAGGGCCAGCACTCGTACGACCAACGCCCCAGGTAAAGCAGCGGCTCTCGCCAGGCCGCGTCCAGCaccacc                                                                  >  1:263549/1‑150 (MQ=255)
                                                                                   aGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGcc                                                               <  2:3892669/150‑1 (MQ=255)
                                                                                   aGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGcc                                                               <  2:3318704/150‑1 (MQ=255)
                                                                                   aGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGcc                                                               <  2:2846342/150‑1 (MQ=255)
                                                                                        tGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGc                                                                                                                >  2:2332921/1‑97 (MQ=255)
                                                                                        tGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGc                                                                                                                <  1:2332921/97‑1 (MQ=255)
                                                                                        tGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGaa                                                          <  2:3432265/150‑1 (MQ=255)
                                                                                        tGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGaa                                                          <  1:5124320/150‑1 (MQ=255)
                                                                                         ggaggaCGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCGGAAc                                                         >  2:5111134/1‑150 (MQ=255)
                                                                                                     gcgcTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTc                                             >  2:2912371/1‑150 (MQ=255)
                                                                                                       gcTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCACGTAAAGCAGCTGCTCTCCCCACGCCTCGTCCTGCACCCCCGCCCTGAACCCCACCACCgaacc                                           >  2:2773336/1‑146 (MQ=255)
                                                                                                            cTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCa                                      <  2:78644/150‑1 (MQ=255)
                                                                                                             tGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCAt                                     <  1:4925831/150‑1 (MQ=255)
                                                                                                               tCCAAGTGAGGTTGCCTGTACCGGCCCTTTGGCGGGCGCGCCCCCCCAACGGGCCCCACGCGTAACACCAAAGCCCCCGGGAAAACAGCGGCTCAGGCCACGCCCCCTCCAGCACCAACCCCCCGAAACCCCCCCCCCGCCCCCGCCagc                                    >  1:2279262/1‑148 (MQ=255)
                                                                                                                ccAAGTGAGGTTGCCTGTTCCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCa                                  <  1:4948987/150‑1 (MQ=255)
                                                                                                                ccAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCa                                  <  2:4796457/150‑1 (MQ=255)
                                                                                                                     tGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCCGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTgcg                             >  2:4763211/1‑150 (MQ=255)
                                                                                                                     tGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTgcg                             >  2:4763199/1‑150 (MQ=255)
                                                                                                                             ccTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTgc                     <  1:4939993/150‑1 (MQ=255)
                                                                                                                                 tACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTgcgcca                 <  2:3002954/150‑1 (MQ=255)
                                                                                                                                   ccGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTgcgc                            >  2:4386815/1‑137 (MQ=255)
                                                                                                                                   ccGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTgcgc                            <  1:4386815/137‑1 (MQ=255)
                                                                                                                                    tggCCCTTTCGCGGGCGAGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTgcgccagcg              >  2:1167204/2‑150 (MQ=255)
                                                                                                                                    cGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACCACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCCCCCAGAACCGCACCACCGTCCCCGGCATGCACAGCGCCAGTTgcgccagcg              >  2:2919354/1‑150 (MQ=255)
                                                                                                                                        ccTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGcgccagcgccag          <  1:3685093/150‑1 (MQ=255)
                                                                                                                                        ccTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGcgccagcgccag          >  1:3092986/1‑150 (MQ=255)
                                                                                                                                        ccTTTCGCGGGCGAGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGcgccagcgccag          >  1:5211032/1‑150 (MQ=255)
                                                                                                                                         cTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCgccagcgccagc         >  2:2694458/1‑150 (MQ=255)
                                                                                                                                         cTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCAGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCgccagcgccagc         >  2:5418313/1‑150 (MQ=255)
                                                                                                                                          tttCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGca                                                     >  1:1324366/1‑105 (MQ=255)
                                                                                                                                          tttCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGca                                                     <  2:1324366/105‑1 (MQ=255)
                                                                                                                                          tttCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGca                                                     >  1:4075501/1‑105 (MQ=255)
                                                                                                                                          tttCGCGGGCGAGCCCGCGAAAGGGCCGGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGccagcgccagcg        >  2:2356675/1‑150 (MQ=255)
                                                                                                                                          tttCGCGGGCGAGCCCGCGAAAGGGCCGGCAGTCTTACGACCAACGCCCCAGGTAAAGCCGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGccagcgccagcg        >  2:3502137/1‑150 (MQ=255)
                                                                                                                                          tttCGCGGGCGAGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGccagcgccagcg        >  2:1880768/1‑150 (MQ=255)
                                                                                                                                          tttCGCGGGCGAGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGccagcgccagcg        >  2:324386/1‑150 (MQ=255)
                                                                                                                                            tCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGCCAGCGCCAgcggc      <  1:1735754/150‑1 (MQ=255)
                                                                                                                                            tCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGCCAGCGCCAgcggc      <  1:38980/150‑1 (MQ=255)
                                                                                                                                            tCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGCCAGCGCCAgcggc      <  2:3458183/150‑1 (MQ=255)
                                                                                                                                              gcgGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCCCCGTCCCCGGCATGCACTGCGCCAGTTGCGCCAGCGCCAgcggcgg    >  1:1394293/1‑149 (MQ=255)
                                                                                                                                                gggCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGCCAGCGCCAGCGGCGTCa  <  2:4382594/150‑1 (MQ=255)
                                                                                                                                                gggCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCAGCTGATCTCGCCACGCCTCGACCTGCACCACCGCCCTGAAACGCACCACCGTCCCCGGCATGCACTGCGCCAGTAGCGCCAGCGCCAGCGGCGTCa  >  1:4510234/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CCCAAGCAGGCCTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGCCAGCGCCAGCGGCGTCA  >  NC_002947/5196709‑5197002

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: