Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,102,008 A→G 7.5% intergenic (+124/+29) 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,102,0080AG 7.5% 278.1 / 9.0 94intergenic (+124/+29)PP_4488/phhRhypothetical protein/sigma‑54 dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (50/36);  new base G (5/2);  total (56/38)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

TGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCAGC  >  NC_002947/5101859‑5102157
                                                                                                                                                     |                                                                                                                                                     
tGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCAGCGGGCCCCTATCGCGGCACAAGGCCCCGCCCACACCTAAAGCGCTCAGCCAGCGCCACGCGGATTACCTGTGGGAGCGGCCGTGTGGCGCGAAAGGACTGCa                                                                                                                                                       >  2:2679273/1‑150 (MQ=255)
 ggggAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCGGCGGGCCCCTATCGCGACACAAGGCCGCTCCCCCAGATAAATCGCGCAGCCAGAGCACAGCGGAATAACAGGGGGAGGGGCCGTTTGTGGCGAAAAGAATGCaa                                                                                                                                                      >  1:2945726/1‑150 (MQ=255)
  gggAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACACATAAAGCGCTCAGCCACAGCCCAGCGGATTACCAGTGGGAGCGGCCATGTGTCGCGGAAGGACGGCAAc                                                                                                                                                     >  2:1457994/1‑150 (MQ=255)
    gAACGGCTGCTCCAACGGCTGATACCTGGGAATCGGGGGCCCGGCGGGCCCCGATCGCGACACCAGGCCCCCCCCCCCGCGAAGGCGCGCCGCCCGCGCCCCGCGGCTTACCGGTGGGGGCGGCCGTGTGGCGCGACAGGCCGGCCACGc                                                                                                                                                   >  2:4013915/1‑150 (MQ=255)
      aCGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAg                                                                                                                                                 <  2:1163320/150‑1 (MQ=255)
         gctgctCCAACGGCTGATACCTTGGAATCGGGGGCCCGGCGGGCCCCGATCGCGGCACAAGGCCGCGCCCACAGATAATGCGCTCAGCCAGAGCACAGCGGATTACCAGTGGGAGCGGCCGTGTGTCGCGAAAGGACTGCAACGCAGGcc                                                                                                                                              >  2:2446565/1‑150 (MQ=255)
          ctgctCCACGGGCTGATACCTTGGATTGGGGGGCCCGGCGGGCCCCTATCGCGACACACGGCCGCTCCCACAGATACTGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCGTTGCAGTcct                                                                                                                                             <  1:1538761/150‑2 (MQ=37)
              tCCAACGGCTGATACCTTGGAAGCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCACCCAAACATAAAGCGCTCAGCCCGAGCACAGCGTATTACCTGTGGGAGCGGCCATGTGTCGCGAAAGGACTGCACCGCAGGCCCAGGt                                                                                                                                         >  2:3280251/1‑150 (MQ=255)
               ccACCGGCTGATACCTTGGATTGGGGGGCCCGGCGGGCCCCTTTCGCCACACCCGGCCGCTCCCACAGATATGGCGCCCACCCCGAGCACTGCGTTTTACGGGTGGGCGCGGCCTTGGGCCGCGAAAGGACTGCAACGCAGTCCCAGGtt                                                                                                                                        <  1:3901753/150‑1 (MQ=255)
                   cGGCTGTTACCTTGGAATGGGGGGCCCTGCGGGCCCCTATCGCCACACACGGCCGCTCCCACAGATAATGCGCTCACCCAGAGCACAGCGTTTTACCTGTGGGCGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCttt                                                                                                                                    <  1:4853635/150‑1 (MQ=255)
                   cGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTTTTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCttt                                                                                                                                    <  1:2035133/150‑1 (MQ=255)
                   cGGCTGATACCTTGGAATCGGGGGCCCGGCGGGCCCCAATCGCGACACAAGGCCCCACCCACAGATAATGCGCTCAGCCAGAGCACAGCGTAATACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCttt                                                                                                                                    >  2:2956657/1‑150 (MQ=39)
                   cGGCTGATACCTTGGAATCGGGGGCCCAGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCCCAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCttt                                                                                                                                    >  2:2956701/1‑150 (MQ=255)
                                 gAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAAAGCGCTCAGCCAGAGCACAGCGGATTACCTGTGGGAGCGGCCATGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTTACAAACAACTaaa                                                                                                                      >  2:3355480/1‑149 (MQ=39)
                                   atggggggCCCGGGGGGCCCCTAACGCCACACACGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCTTATTACCTGTGGGCGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATcc                                                                                                                    <  1:3144810/147‑1 (MQ=255)
                                    tCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCGCACCCCGAGCACCGCGGATTACCGGTGGGAGCGGCCGTGGGTCGCGAAAGGCCTGCACCGCAGTCCCAGGTGCTGTGACCCACAACGAATccg                                                                                                                   >  2:2672563/1‑150 (MQ=255)
                                     taggggCCCGCAGGGCCCCTATCGCGACACAAGGCCCCTCCCACAGATACTGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACGGCAACGCAGTCCCAGGTTCTGTGGCCAACAACTAATccgc                                                                                                                  >  2:3920313/3‑150 (MQ=38)
                                        gggCCCGGGGGGCCCCTAGCGCGCCCCACGGCCGCTCCCACCTATAATGCGCTCAGCCAGAGCAAAGCTTTTTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGa                                                                                                               <  1:4928541/150‑1 (MQ=255)
                                         ggCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGAc                                                                                                              <  2:962707/150‑1 (MQ=255)
                                                 gggCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGAGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAg                                                                                                      >  1:4659610/1‑150 (MQ=255)
                                                     cccTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCCGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAg                                                                                                  <  1:1253793/150‑1 (MQ=255)
                                                     cccTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAg                                                                                                  <  2:226599/150‑1 (MQ=255)
                                                       cTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCttt                                                                                                                                    >  1:2910360/1‑114 (MQ=255)
                                                       cTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCttt                                                                                                                                    <  2:2910360/114‑1 (MQ=255)
                                                        tATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTaa                                                                                               >  2:2018755/1‑150 (MQ=255)
                                                        tATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTaa                                                                                               >  2:4487484/1‑150 (MQ=255)
                                                         agcccgACCCACGGCCGCTCCCACAGATCATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAAt                                                                                              <  2:3199482/146‑1 (MQ=255)
                                                         aTCGCGACCCCAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCGGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAAt                                                                                              <  2:2657028/150‑1 (MQ=255)
                                                         aTCGCGACACACGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAAt                                                                                              <  1:4482603/150‑1 (MQ=255)
                                                            gcgACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGtt                                                                                                                                        <  2:1601069/105‑1 (MQ=255)
                                                            gcgACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGtt                                                                                                                                        >  1:1601069/1‑105 (MQ=255)
                                                            gcgACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATccc                                                                                           >  2:185638/1‑150 (MQ=255)
                                                            gcgACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATccc                                                                                           >  2:4084265/1‑150 (MQ=255)
                                                                cacaAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCACCGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCaa                                                                                       >  2:1161750/1‑150 (MQ=255)
                                                                cacaAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCaa                                                                                       >  2:3755323/1‑150 (MQ=255)
                                                                   aaGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACtt                                                                                    >  1:3332043/1‑150 (MQ=255)
                                                                      gCCGCTCCCCCAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTATTCCCTCAACttgtt                                                                                 <  1:2446565/150‑1 (MQ=255)
                                                                      gCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACttgtt                                                                                 <  2:2456488/150‑1 (MQ=255)
                                                                         gCTCCCACAGGTAATACGCTGTGCTCTGGCTGAGCGCATTATCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGc                                                                              >  1:1860138/1‑150 (MQ=21)
                                                                             ccACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAAtt                                                                          <  2:4326653/150‑1 (MQ=255)
                                                                             ccACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAAtt                                                                          >  2:2869948/1‑150 (MQ=255)
                                                                              cacaGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTg                                                                         >  1:4698331/1‑150 (MQ=255)
                                                                              cacaGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTg                                                                         >  1:2211166/1‑150 (MQ=255)
                                                                               acaGGTAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGt                                                                        >  2:5047023/1‑150 (MQ=255)
                                                                               acaGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGTCCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGt                                                                        >  2:5081515/1‑150 (MQ=255)
                                                                               acaGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGt                                                                        >  1:1255945/1‑150 (MQ=255)
                                                                                cataTAATGCGCTCACCCAGAGCACAGCGTTTTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTc                                                                       <  1:1298912/147‑1 (MQ=255)
                                                                                  gATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCgt                                                                     <  2:3328024/150‑1 (MQ=255)
                                                                                  gATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCgt                                                                     <  2:2291924/150‑1 (MQ=255)
                                                                                   aTAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCgtg                                                                    >  2:2327097/1‑150 (MQ=255)
                                                                                      aTGCGCCCACCCAGAGCACAGCTTTTTACCGGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGc                                                                 <  1:2805912/150‑1 (MQ=255)
                                                                                                ccAGAGCACCGCTTTTTACGTGTGGGAGCGGCGTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCa                                                       <  1:3856624/150‑1 (MQ=255)
                                                                                                ccAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCGTTGCAGTcct                                                                                                                                             >  2:4797140/1‑63 (MQ=21)
                                                                                                ccAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCGTTGCAGTcct                                                                                                                                             <  1:4797140/64‑2 (MQ=21)
                                                                                                ccAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCa                                                       >  2:4017306/1‑150 (MQ=255)
                                                                                                 cAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAg                                                      >  1:893903/1‑150 (MQ=255)
                                                                                                         aGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGc                                              >  2:3465254/1‑150 (MQ=255)
                                                                                                             tATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCaa                                          >  2:495338/1‑150 (MQ=255)
                                                                                                              aTTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAg                                         <  2:1952112/150‑1 (MQ=255)
                                                                                                               ttACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAgg                                        <  2:2766349/150‑1 (MQ=255)
                                                                                                                  tcTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGccc                                     <  1:4606536/149‑1 (MQ=255)
                                                                                                                  ccTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGccc                                     <  1:3627618/150‑1 (MQ=255)
                                                                                                                   cTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCg                                    >  1:500841/1‑150 (MQ=255)
                                                                                                                    tgtgGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGa                                   >  2:939540/1‑150 (MQ=255)
                                                                                                                          aGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTg                             <  1:1640374/150‑1 (MQ=255)
                                                                                                                            cGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTggg                           <  1:1457994/150‑1 (MQ=255)
                                                                                                                              gCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGaa                         <  1:2988620/150‑1 (MQ=255)
                                                                                                                                    tgccgcgAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTccc                   <  1:1438480/147‑1 (MQ=255)
                                                                                                                                     gtCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTcccc                  <  2:2945726/150‑1 (MQ=255)
                                                                                                                                        gcgAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTCCTGGGAAACTCCCCCAt               >  2:4106313/1‑150 (MQ=255)
                                                                                                                                        gcgAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCAt               >  1:967034/1‑150 (MQ=255)
                                                                                                                                        gcgAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCAt               >  2:1144130/1‑150 (MQ=255)
                                                                                                                                          gAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCa             >  1:4991983/1‑150 (MQ=255)
                                                                                                                                          gAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCa             >  1:2891145/1‑150 (MQ=255)
                                                                                                                                           ctgggACTGCGTTGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCaa            >  1:2213229/4‑150 (MQ=255)
                                                                                                                                           ctgggACTGCGTTGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAAATCCCCCATCaa            >  1:2213009/4‑150 (MQ=255)
                                                                                                                                            tgggACTGCGTTGCAGTCCTAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCaaa           >  1:2127953/3‑150 (MQ=255)
                                                                                                                                            tgggACTGCGTTGCAGTCCTAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCaaa           >  1:2127825/3‑150 (MQ=255)
                                                                                                                                             aGGACTGCAACGCAGTCCCAGGTTCTTTGTCCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAAc          <  1:5081515/150‑1 (MQ=255)
                                                                                                                                             aGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAAc          >  1:50534/1‑150 (MQ=255)
                                                                                                                                             aGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACGTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCCTCAAAc          >  2:3841234/1‑150 (MQ=255)
                                                                                                                                                  tGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTACTCCCTCAACTTGTTGGCCATTGTCGTGTGCGACACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGACACTCCCCCATCAAACTTTcc     >  2:3636992/1‑150 (MQ=255)
                                                                                                                                                  tGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTcc     >  2:1179006/1‑150 (MQ=255)
                                                                                                                                                  tGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTcc     >  1:579516/1‑150 (MQ=255)
                                                                                                                                                  tGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTcc     >  1:550194/1‑150 (MQ=255)
                                                                                                                                                  tGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTcc     <  1:2869948/150‑1 (MQ=255)
                                                                                                                                                   gCAACGCAGTCCCAGGTTCTTTGACCAACAACTACTCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCa    <  1:3828122/150‑1 (MQ=255)
                                                                                                                                                   gCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCa    <  1:2679273/150‑1 (MQ=255)
                                                                                                                                                   gCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCa    >  2:3237198/1‑150 (MQ=255)
                                                                                                                                                   gCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCa    <  1:3280251/150‑1 (MQ=255)
                                                                                                                                                    cAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCCAACTTTCCAg   >  2:3360343/1‑150 (MQ=255)
                                                                                                                                                    cAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCAg   >  2:1343143/1‑150 (MQ=255)
                                                                                                                                                    cAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCAg   >  1:2109791/1‑150 (MQ=255)
                                                                                                                                                    cAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCAg   >  1:1704784/1‑150 (MQ=255)
                                                                                                                                                    cAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCAg   >  1:1705045/1‑150 (MQ=255)
                                                                                                                                                    cAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTCAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCAg   >  1:74793/1‑150 (MQ=255)
                                                                                                                                                     aaCGCAGTCCCTGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCAGc  >  2:1217921/1‑150 (MQ=255)
                                                                                                                                                     aaCGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCAGCAAACTTTCCAGc  >  2:2902595/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCAGC  >  NC_002947/5101859‑5102157

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: ‑