Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,956,632 T→G 8.8% intergenic (‑170/‑111) 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,6320TG 8.8% 275.6 / 13.0 91intergenic (‑170/‑111)PP_2586/PP_2587putative amino acid transporter/LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (55/28);  new base G (5/3);  total (60/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.01e-01

GCAGGACACCGCGGCAGGCGAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCAGAACAATGGTTCGAGC  >  NC_002947/2956484‑2956779
                                                                                                                                                     |                                                                                                                                                    
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                                                                                                                      <  1:2008025/150‑1 (MQ=255)
                              cgACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGt                                                                                                                      <  2:3560714/150‑1 (MQ=255)
                              cgACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGt                                                                                                                      <  2:3560694/150‑1 (MQ=255)
                                   gTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTg                                                                                                                 >  1:4358576/1‑150 (MQ=255)
                                   gTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTg                                                                                                                 >  2:2162395/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:1572549/1‑150 (MQ=255)
                                                                                              tCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGaa                                                      >  2:626514/1‑150 (MQ=255)
                                                                                              tCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGaa                                                      >  2:626526/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:2954000/1‑150 (MQ=255)
                                                                                                                         aCCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACcaatcaa                           >  1:4933070/1‑150 (MQ=255)
                                                                                                                         aCCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACcaatcaa                           >  1:2511709/1‑150 (MQ=255)
                                                                                                                           ccTTTGGGGTGAGCCGGGCAAGGGCTGTGTCCGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAAtt                         >  1:2036814/1‑150 (MQ=25)
                                                                                                                           ccTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAAtt                         >  1:3134204/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)
                                                                                                                             tttGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGc                       >  2:3787139/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       >  1:4617350/1‑150 (MQ=255)
                                                                                                                                              acaGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCAGAACAATGGtt       >  2:4288643/1‑150 (MQ=255)
                                                                                                                                              acaGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCAGAACAATGGtt       >  2:4288670/1‑150 (MQ=255)
                                                                                                                                              acaGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCAGAACAATGGtt       >  2:3254374/1‑150 (MQ=255)
                                                                                                                                              acaGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCAGAACAATGGtt       >  1:2773353/1‑150 (MQ=255)
                                                                                                                                              acaGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCAGAACAATGGtt       >  2:645023/1‑150 (MQ=255)
                                                                                                                                                aGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCAGAACAATGGTTCg     >  2:3477011/1‑150 (MQ=255)
                                                                                                                                                aGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCAGAACAATGGTTCg     >  1:5122316/1‑150 (MQ=255)
                                                                                                                                                aGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCAGAACAATGGTTCg     >  2:904631/1‑150 (MQ=255)
                                                                                                                                                  cccTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCAGAACAATGGTTCGAg   <  2:4154221/150‑1 (MQ=255)
                                                                                                                                                   ccTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCAGAACAATGGTTCGAGc  >  2:1133144/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCAGGACACCGCGGCAGGCGAGCAGGCTCGCGACCGTCAGGTGTTGTGCAAGGAGGTCGAGCGGAAAGCTCGCAGCGCTGGGGCGGGTTGATGCTCCCACACCCTGCAGGCACCCGGAATCACCCTTTGGGGTGAGCGGACACAGCCCTTGCCCGGCTCACCCCAAAGGGTGATTCCGGTGCCTGCAGCGCCGTGGGAGTATCTGCTCAACGACGCACCGCACTGGCGGCGGTGCCACTGAGAAAACGGCAGGCTTACAGTGACCAATCAATTGCTTTCAGAACAATGGTTCGAGC  >  NC_002947/2956484‑2956779

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: