Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,101,905 T→G 14.6% intergenic (+21/+132) PP_4488 → / ← phhR hypothetical protein/sigma‑54 dependent transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,101,9050TG14.6% 186.9 / 6.8 82intergenic (+21/+132)PP_4488/phhRhypothetical protein/sigma‑54 dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (14/53);  new base G (10/4);  total (25/57)
Fisher's exact test for biased strand distribution p-value = 4.67e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTATTTCCATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCC  >  NC_002947/5101756‑5102046
                                                                                                                                                     |                                                                                                                                              
gTTATTTCCATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCt                                                                                                                                                <  2:4780068/150‑1 (MQ=255)
gTTATTTCCATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCt                                                                                                                                                <  2:4779989/150‑1 (MQ=255)
     ttCCATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCggg                                                                                                                                           <  2:3415532/150‑1 (MQ=255)
        cATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGccc                                                                                                                                        <  2:3359235/150‑1 (MQ=255)
        cATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGccc                                                                                                                                        >  1:4801417/1‑150 (MQ=255)
        cATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGccc                                                                                                                                        >  2:3773655/1‑150 (MQ=255)
        cATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGccc                                                                                                                                        <  1:3908732/150‑1 (MQ=255)
        cATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGccc                                                                                                                                        <  2:4127454/150‑1 (MQ=255)
         aTGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGcccc                                                                                                                                       >  2:187113/1‑150 (MQ=255)
          tGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCGGCGGGCCCCt                                                                                                                                      >  2:4160889/1‑150 (MQ=255)
          tGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCGCAGGGccccc                                                                                                                                      <  2:2381230/150‑2 (MQ=255)
            aGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTAt                                                                                                                                    >  2:5233421/1‑150 (MQ=255)
            aGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTAt                                                                                                                                    <  1:5233421/150‑1 (MQ=255)
                aaGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCGGCGGGCCCCTATcgcg                                                                                                                                >  2:2744865/1‑150 (MQ=255)
                  gTGCCCCCTGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGACTACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGac                                                                                                                              <  2:4261307/150‑1 (MQ=255)
                  gTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGACTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGac                                                                                                                              <  1:1021909/150‑1 (MQ=255)
                   tGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGGCCGGCGGGCCCCAATcgcgccc                                                                                                                             >  1:1363270/1‑147 (MQ=255)
                      ccAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACaa                                                                                                                          <  1:5261470/150‑1 (MQ=255)
                      ccAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGACCCTGCGGGCCCCTATCGCGACACaa                                                                                                                          <  2:2786296/150‑1 (MQ=255)
                       cAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATAGCTTGGAATCGGGGGCCCTGCGGGGCCCTATCGGGAgacatg                                                                                                                         >  2:3273705/1‑148 (MQ=255)
                         aTGTGTGCGAAGTCGTGCCGGCCCTCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGTCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGc                                                                                                                       <  1:252502/150‑1 (MQ=255)
                         aTGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGc                                                                                                                       <  2:4750571/150‑1 (MQ=255)
                              tgCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGGAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTc                                                                                                                 <  2:4183753/150‑1 (MQ=255)
                              tgCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTc                                                                                                                 <  1:4717602/150‑1 (MQ=255)
                              tgCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTc                                                                                                                 <  2:4653608/150‑1 (MQ=255)
                                cGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTccc                                                                                                               <  1:4395497/150‑1 (MQ=255)
                                 gAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCAGCGGGCCCCAATAGCGACACAACGCCGCACCca                                                                                                              >  2:3646098/1‑150 (MQ=255)
                                      cGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGAt                                                                                                         <  2:1336991/150‑1 (MQ=255)
                                         gCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAAt                                                                                                      <  1:4501796/150‑1 (MQ=255)
                                         gCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAAt                                                                                                      <  1:4389820/150‑1 (MQ=255)
                                            ggCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATgct                                                                                                   <  2:3107335/150‑2 (MQ=255)
                                             gCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATgcgc                                                                                                  <  2:2438640/150‑1 (MQ=255)
                                               ccGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCCCAGATAATGCGCTc                                                                                                >  2:1754370/1‑150 (MQ=255)
                                               ccGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCGGCGGGCCCCAATCGCGACACAAGGCCGCTCCCAAAGATAATGCGCTc                                                                                                >  2:2571354/1‑150 (MQ=255)
                                                cGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCa                                                                                               <  2:2311545/150‑1 (MQ=255)
                                                cGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCa                                                                                               <  2:1946945/150‑1 (MQ=255)
                                                  cATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGc                                                                                             <  1:660513/150‑1 (MQ=255)
                                                     gACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCag                                                                                          <  2:1127048/150‑1 (MQ=255)
                                                     gACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCag                                                                                          <  2:2604005/150‑1 (MQ=255)
                                                      aCCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCaga                                                                                         <  1:3434457/150‑1 (MQ=255)
                                                      aCCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCaga                                                                                         <  2:2303397/150‑1 (MQ=255)
                                                      aCCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCaga                                                                                         <  1:508728/150‑1 (MQ=255)
                                                       ccTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCagag                                                                                        <  1:2940017/150‑1 (MQ=255)
                                                        cTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCGGCGGGCCCCGATCGCGACACAAGGCCCCACCCACACATAAATCGCTCAGCCAGAGc                                                                                       >  2:156928/1‑150 (MQ=255)
                                                         tACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGca                                                                                      <  1:4368441/150‑1 (MQ=255)
                                                           cGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCTGGCTGAGcgca                                                                                    <  1:1058864/150‑1 (MQ=255)
                                                              ggACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCg                                                                                 <  2:3082223/150‑1 (MQ=255)
                                                                 cTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTAt                                                                              <  1:908231/150‑1 (MQ=255)
                                                                 cTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTAt                                                                              <  1:3185599/150‑1 (MQ=255)
                                                                  tCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTccc                                                                                                               >  2:898186/1‑116 (MQ=255)
                                                                  tCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTccc                                                                                                               <  1:898186/116‑1 (MQ=255)
                                                                  tCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCGGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACCCATAATGCGCTCAGCCACAGCACAGCGTAtt                                                                             >  2:687515/1‑150 (MQ=255)
                                                                    ccAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTAc                                                                           >  2:5409261/1‑150 (MQ=255)
                                                                    ccAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAACACAGCGTATTAc                                                                           <  2:1329935/150‑1 (MQ=255)
                                                                     cAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCGGCGGGCCCCGATCGCGACACAACGCCGCGCCCCCAGATAAAGCGCGCAGCCAGAGCCCAGCGGATTAcc                                                                          >  1:2190812/1‑150 (MQ=255)
                                                                          ccAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCtgtg                                                                      <  2:402928/150‑1 (MQ=255)
                                                                          ccAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCtgtg                                                                      <  2:3131865/150‑1 (MQ=255)
                                                                            aGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGa                                                                   <  1:2907634/150‑1 (MQ=255)
                                                                            aGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGa                                                                   <  1:735290/150‑1 (MQ=255)
                                                                                gCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCAATCGCGCCACAAGGCCGCTCCCACAGATAAAGCGCTCAGCCAGAGCCCAGCGGAGTACCTGTGGGAGCgg                                                               >  1:4084170/1‑150 (MQ=255)
                                                                                    aTTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCtt                                                           <  1:3646098/150‑1 (MQ=255)
                                                                                     ttCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCAATCGCGACACAAGGCCGCACCCACACAAAAATCGCTCACCCCCAGCACAACGGATTACCAGTGGGAGCCGCCGtg                                                          >  1:3617735/1‑150 (MQ=255)
                                                                                      tctcGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTtgt                                                         <  2:2483914/150‑1 (MQ=255)
                                                                                        tcGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGCCACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTgtgt                                                       <  1:2738875/150‑1 (MQ=255)
                                                                                        tcGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTCTTACCTGTGGGAGCGGCCTTgtgt                                                       >  2:3076612/1‑150 (MQ=255)
                                                                                        tcGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTgtgt                                                       >  2:4489245/1‑150 (MQ=255)
                                                                                         cGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTc                                                      <  2:1629254/150‑1 (MQ=255)
                                                                                               ccGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGGTAATACGCTGTGCTCTGGCTGAGCGCATTATCTGTGGGAGCGGCCTTGTGTCGCGata                                                <  1:3061484/150‑3 (MQ=17)
                                                                                                   gTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATAGCCGCACAAGGCCCCACCCACAGATAATTCGCACAGCCCGAGCACAGCGGATTAACAGTGGGGGCGGCCATGTGTCGCGTAAgggc                                            >  2:1414952/1‑148 (MQ=255)
                                                                                                        ggggAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCGGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAAATCGCACAGCCAGAGCACACCGTATTAACGGTGGGAGCGGCCTTGTGTAGCGGAAGGACAGCaa                                       >  2:618416/1‑150 (MQ=255)
                                                                                                                   gctCCAACGGCTGATACCTTGGAATGGGGGGCCCTGCGGGCCCCTATCGCGACCCACGGCCGCTCCCACATATAATGCGCTCAGCCAGAGAACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAg                            <  1:5278664/150‑1 (MQ=39)
                                                                                                                   cccccACCGGCTGATACTTTGGACTGGGGGCCCGGGCGGGCCCCTATCGCGCCACACGGCCGCCCCCACAGATATGGCGCTCCCCCAGCGAACAGCTTATTACCTGGGGGAGGGGCCTTGTCTCGCGAAAGGACTGCAACGCAGTCCCAg                            <  1:2744865/147‑1 (MQ=255)
                                                                                                                    ctCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCGTTGCAGTcctttc                           <  2:295571/150‑5 (MQ=255)
                                                                                                                    ctCCAACGGCTGATACCTTGGAATCGGGGGCCCGGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGCAAGGACTGCAACGCAGTCCCAgg                           >  2:3098216/1‑150 (MQ=255)
                                                                                                                     tCCAACGGCTGATACCGTGGAATCGGGGGCCCGGCGGGCCCCGATCGCGACACAAGGCCGCTCCCACAGCAAATGCGCACAGCCCCAGCACAGCGGATTACCGGGGGGAGCGGCCGGGTGTCGCGAAAGGACTGCAACGCCGTCCCCggg                          >  2:825410/1‑149 (MQ=255)
                                                                                                                            gCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGa                   <  2:5390807/150‑1 (MQ=255)
                                                                                                                                aTACCTTGTAAGGGGGGGCCCGGCGGGCCCCTACCGCGCCCCCCGCCCCCGCCCAAATATAAGGCGCTCCGCCAGAGCACAGCTTTTTACGTGTGGGAGGGGCGTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACcaa               <  1:987564/150‑1 (MQ=255)
                                                                                                                                     ttGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCACAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACt          <  2:311736/150‑1 (MQ=255)
                                                                                                                                     ttGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTACCACAGATAATTCGCTCAGCCAGAGCACAGCGTATTAACTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACcaacaaat          >  2:5393049/1‑148 (MQ=255)
                                                                                                                                       gcattggggggCCCGGGGGGCCCCTTTCGCGCCACACGGCCGCTCCCACATATAATGCGCTCAGCCTGCGCACAGGTTTTTACCTGTGGGTGGGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTaa        <  1:528585/144‑1 (MQ=255)
                                                                                                                                          atggggggCCCTGCGGGCCCCTATCGCGACACCAGGCCGATCCCACATATAATGCGCACCGCCAGAGCACAGCGTATTACGTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATcc     <  1:4160889/147‑1 (MQ=255)
                                                                                                                                          aTCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATcc     <  2:3128046/150‑1 (MQ=255)
                                                                                                                                             gggggCCCTGCGGGCCCCTATCGCGACACACGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACCGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATccgcc  <  2:262494/150‑1 (MQ=255)
                                                                                                                                             gggggCCCTGCGGGCCCCTATCGCGACACACGGCCGCTCCAACAGAAAATGCGCTCAGCCAGAGCACAGCTTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATccgcc  <  1:1231632/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
GTTATTTCCATGAGCGAAGTGCCCAATGTGTGCGAAGTCGTGCCGGCCCGCATGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCC  >  NC_002947/5101756‑5102046

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: