Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,956,626 C→A 10.4% intergenic (‑164/‑117) PP_2586 ← / → PP_2587 putative amino acid transporter/LuxR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,956,6260CA10.4% 260.9 / 15.4 84intergenic (‑164/‑117)PP_2586/PP_2587putative amino acid transporter/LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (47/28);  new base A (6/3);  total (53/31)
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.92e-01

AAGGGATGCAGGACACCGCGGCAGGCGAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCAGAACAATGGTT  >  NC_002947/2956477‑2956774
                                                                                                                                                      |                                                                                                                                                     
aaGGGATGCAGGACACCGCGGCAGGCGAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGacac                                                                                                                                                       <  1:1638696/150‑1 (MQ=255)
aaGGGATGCAGGACACCGCGGCAGGCGAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGacac                                                                                                                                                       <  1:2186586/150‑1 (MQ=255)
    gATGCAGGACACCGCGGCAGGCGAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGcc                                                                                                                                                   >  2:4331655/1‑150 (MQ=255)
    gATGCAGGACACCGCGGCAGGCGAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGcc                                                                                                                                                   >  2:595125/1‑150 (MQ=255)
       gCAGGACACCGCGGCAGGCGAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCtt                                                                                                                                                <  2:4133115/150‑1 (MQ=255)
       gCAGGACACCGCGGCAGGCGAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCAAACCATGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACCAAGCCAtt                                                                                                                                                >  2:2097160/1‑150 (MQ=255)
         aGGACACCGCGGCAGGCGAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGc                                                                                                                                              >  1:1526671/1‑150 (MQ=255)
           gACACCGCGGCAGGCGAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAACGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTGTGGGGTGAGCGGACACAGCCCTGGccc                                                                                                                                            >  2:1506243/1‑150 (MQ=255)
             cacCGCGGCAGGCGAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCgg                                                                                                                                          <  1:3496236/150‑1 (MQ=255)
                  cgGCAGGCGAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGTAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCAc                                                                                                                                     <  1:159702/150‑1 (MQ=255)
                  cgGCAGGCGAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCAc                                                                                                                                     <  1:2638609/150‑1 (MQ=255)
                  cgGCAGGCGAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCAc                                                                                                                                     <  1:2638798/150‑1 (MQ=255)
                   ggCAGGCGAGCGGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCAcc                                                                                                                                    <  2:5097239/150‑1 (MQ=255)
                    ggtcgCGAGCCTGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCAccc                                                                                                                                   >  1:4386468/5‑150 (MQ=255)
                     gtcgCGAGCCTGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCAcccc                                                                                                                                  >  1:3285405/4‑150 (MQ=255)
                      aGGCGAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCa                                                                                                                                 >  1:842894/1‑150 (MQ=255)
                      aGGCGAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCa                                                                                                                                 >  1:3414573/1‑150 (MQ=255)
                        gCGAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCaaa                                                                                                                               <  2:2947297/150‑1 (MQ=255)
                          gAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAgg                                                                                                                             <  1:4990989/150‑1 (MQ=255)
                          gAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAgg                                                                                                                             <  1:4990680/150‑1 (MQ=255)
                           aGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCACAggg                                                                                                                            >  1:4029367/1‑150 (MQ=255)
                           aGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAggg                                                                                                                            >  1:4029277/1‑150 (MQ=255)
                            gCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTTGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCCGGGCAAGGGCTGTGTCCGCTCACCCCAAAGGGt                                                                                                                           <  2:4002407/150‑1 (MQ=21)
                            gCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCCGGGCAAGGGCTGTGTCCGCTCACCCCAAAGGGt                                                                                                                           <  2:4002424/150‑1 (MQ=21)
                            gCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCGGAATCACCCTTTGGGGTGAGCCGGGCAAGGGCTGTGTCCGCTCACCCCAAAGGGTg                                                                                                                          <  1:4171164/150‑1 (MQ=18)
                               ggCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGAt                                                                                                                        >  1:1017678/1‑150 (MQ=255)
                                  tCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTcc                                                                                                                     >  1:4300603/1‑150 (MQ=255)
                                    gcgACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCCGCAGGCACCGGAATCACCCTTTGGGGTGAGCCGGGCAAGGGCTGTGTCCGCTCACCCCAAAGGGTGATTCCggg                                                                                                                 >  2:4220582/1‑149 (MQ=18)
                                     cgACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGt                                                                                                                 <  2:3560714/150‑1 (MQ=255)
                                     cgACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGt                                                                                                                 <  2:3560694/150‑1 (MQ=255)
                                     cgACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGt                                                                                                                 <  1:2008025/150‑1 (MQ=255)
                                          gTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTg                                                                                                            >  2:2162395/1‑150 (MQ=255)
                                          gTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTg                                                                                                            >  1:4358576/1‑150 (MQ=255)
                                                   tgtgCACGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAAGCACCCTTTGGGGTGAGCGGACACAGCCCTGGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGt                                                                                                   <  2:4386157/150‑1 (MQ=255)
                                                   tgtgCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTTAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGt                                                                                                   <  2:686659/150‑1 (MQ=255)
                                                      gCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTggg                                                                                                >  2:3124457/1‑150 (MQ=255)
                                                        aaGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAg                                                                                              >  2:1949007/1‑150 (MQ=255)
                                                        aaGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAg                                                                                              >  2:1336220/1‑150 (MQ=255)
                                                                      gAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAacga                                                                                >  2:2141167/1‑150 (MQ=255)
                                                                       aaaGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGGCCTTGCCCGGCTCACCCGAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGATCAacgac                                                                               >  2:3115247/1‑150 (MQ=255)
                                                                           cTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGAcgcac                                                                           <  2:1616018/150‑1 (MQ=255)
                                                                            tCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGAcgcacc                                                                          >  1:805645/1‑150 (MQ=255)
                                                                               cAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGAcgcaccgca                                                                       >  2:452595/1‑150 (MQ=255)
                                                                                    cTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTggc                                                                  >  1:1846382/1‑150 (MQ=255)
                                                                                    cTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTggc                                                                  >  1:2513605/1‑150 (MQ=255)
                                                                                         gCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCGGGCTCACCCCAAAGGGTTATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGcggcgg                                                             <  2:1386910/150‑1 (MQ=255)
                                                                                            ggTTGATGCCCCCACACCCTGCAGGCACCCGCAATCACCCTTTGGGGTGAGCGGACACCGCCCTTGCCCGGAACCCCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGc                                                          <  1:4538648/150‑1 (MQ=255)
                                                                                                 aTGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTg                                                     <  2:3290334/150‑1 (MQ=255)
                                                                                                  tGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTga                                                    <  2:1779450/150‑1 (MQ=255)
                                                                                                    cTCCAACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGAACCGCACTGGCGGCGGTGCCACTgaga                                                  >  2:2214767/1‑150 (MQ=255)
                                                                                                     tCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGGCGCACCGCACGGGCGGCGGTGCCCCTGAGaa                                                 >  1:4389622/1‑150 (MQ=255)
                                                                                                     tCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGaa                                                 >  2:626526/1‑150 (MQ=255)
                                                                                                     tCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGaa                                                 >  2:1572549/1‑150 (MQ=255)
                                                                                                     tCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGaa                                                 >  2:626514/1‑150 (MQ=255)
                                                                                                      cccACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGCGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGaaa                                                >  2:1722741/1‑150 (MQ=255)
                                                                                                      cccACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGaaa                                                >  1:5087146/1‑150 (MQ=255)
                                                                                                        cacaCCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACCCAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGCGAAAAc                                              >  1:3548143/1‑150 (MQ=255)
                                                                                                           acCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGTCACTGAGAAAACGGc                                           <  2:4797795/150‑1 (MQ=255)
                                                                                                           acCCGGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGc                                           <  2:1365066/150‑1 (MQ=255)
                                                                                                                     cACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAg                                 >  2:3445503/1‑150 (MQ=255)
                                                                                                                      aCCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGGGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGt                                >  1:4706893/1‑150 (MQ=255)
                                                                                                                       accGGAATCACCCTTTGGGGTGAGCCGGGCAAGGGCTGTGTCCGCTCACCCCAAAGGGTGATTCCGGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGt                                >  1:4955550/2‑150 (MQ=17)
                                                                                                                        ccGGAATCACCCTTGGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAATGGGTTATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGa                              <  1:2072781/150‑1 (MQ=255)
                                                                                                                         cGGAATCACCCTTTGGGGTGCGCGGACACAGCCCTTGCCCGGCACACCCCAAAGGGTTATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGAc                             <  1:5052218/150‑1 (MQ=255)
                                                                                                                              tCACCCTTTGGGGTGAGCCGGGCAAGGGCTGTGTCCGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACcaatc                        >  1:3109950/1‑150 (MQ=21)
                                                                                                                              tCACCCTTTGGGGTGAGCCGGGCAAGGGCTGTGTCCGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACcaatc                        >  1:554681/1‑150 (MQ=21)
                                                                                                                                aCCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACcaatcaa                      >  1:4933008/1‑150 (MQ=255)
                                                                                                                                aCCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACcaatcaa                      >  1:4933070/1‑150 (MQ=255)
                                                                                                                                aCCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACcaatcaa                      >  1:2511709/1‑150 (MQ=255)
                                                                                                                                aCCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACcaatcaa                      >  1:2954000/1‑150 (MQ=255)
                                                                                                                                  ccTTTGGGGTGAGCCGGGCAAGGGCTGTGTCCGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAAtt                    >  1:2036814/1‑150 (MQ=25)
                                                                                                                                  ccTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAAtt                    >  1:3134204/1‑150 (MQ=255)
                                                                                                                                    tttGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGc                  >  2:3787139/1‑150 (MQ=255)
                                                                                                                                    tttGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGc                  >  2:3017956/1‑150 (MQ=255)
                                                                                                                                    tttGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGc                  >  1:3939637/1‑150 (MQ=255)
                                                                                                                                    tttGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGc                  >  1:3939730/1‑150 (MQ=255)
                                                                                                                                       ggggTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCACTTGCtt                >  1:4870187/1‑149 (MQ=255)
                                                                                                                                       ggggTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCtt                <  2:4870187/149‑1 (MQ=255)
                                                                                                                                         ggTGAGCGGACACAGCCCTTGCCCGGCTCCCCCAAATGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCa             <  1:2097160/150‑1 (MQ=255)
                                                                                                                                         ggTGACCGCACACCGCCCTTCCCCCCCTCACCCCAAAGGCTGATACCGGTCCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCa             <  2:2601385/150‑1 (MQ=255)
                                                                                                                                          gTGCGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCAg            <  1:5011556/150‑1 (MQ=255)
                                                                                                                                              gCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCAGAACa        >  1:587214/1‑150 (MQ=255)
                                                                                                                                                cGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCAGAACaa       >  2:1640842/1‑150 (MQ=255)
                                                                                                                                                 ggACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCAGAACAAt      <  1:4193306/150‑1 (MQ=255)
                                                                                                                                                 ggACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCAGAACAAt      >  2:1076490/1‑150 (MQ=255)
                                                                                                                                                 ggACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCAGAACAAt      >  2:1214700/1‑150 (MQ=255)
                                                                                                                                                     acaGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCAGAACAATGGtt  >  2:4288670/1‑150 (MQ=255)
                                                                                                                                                     acaGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCAGAACAATGGtt  >  2:4288643/1‑150 (MQ=255)
                                                                                                                                                     acaGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCAGAACAATGGtt  >  1:2773353/1‑150 (MQ=255)
                                                                                                                                                     acaGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCAGAACAATGGtt  >  1:4617350/1‑150 (MQ=255)
                                                                                                                                                     acaGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCAGAACAATGGtt  >  2:645023/1‑150 (MQ=255)
                                                                                                                                                     acaGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCAGAACAATGGtt  >  2:3254374/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AAGGGATGCAGGACACCGCGGCAGGCGAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCAGAACAATGGTT  >  NC_002947/2956477‑2956774

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: