Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,102,009 A→T 7.4% intergenic (+125/+28) 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,0090AT 7.4% 295.1 / 11.4 94intergenic (+125/+28)PP_4488/phhRhypothetical protein/sigma‑54 dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (48/36);  new base T (5/2);  total (56/38)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCAGCA  >  NC_002947/5101860‑5102158
                                                                                                                                                     |                                                                                                                                                     
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                                                                                                                                     <  2:2910360/114‑1 (MQ=255)
                                                      cTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCttt                                                                                                                                     >  1:2910360/1‑114 (MQ=255)
                                                       tATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTaa                                                                                                >  2:4487484/1‑150 (MQ=255)
                                                       tATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTaa                                                                                                >  2:2018755/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:4084265/1‑150 (MQ=255)
                                                           gcgACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATccc                                                                                            >  2:185638/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                                                                                                                                              <  1:4797140/64‑2 (MQ=21)
                                                                                               ccAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCGTTGCAGTcct                                                                                                                                              >  2:4797140/1‑63 (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:2891145/1‑150 (MQ=255)
                                                                                                                                         gAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCa              >  1:4991983/1‑150 (MQ=255)
                                                                                                                                          ctgggACTGCGTTGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCaa             >  1:2213229/4‑150 (MQ=255)
                                                                                                                                          ctgggACTGCGTTGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAAATCCCCCATCaa             >  1:2213009/4‑150 (MQ=255)
                                                                                                                                           tgggACTGCGTTGCAGTCCTAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCaaa            >  1:2127825/3‑150 (MQ=255)
                                                                                                                                           tgggACTGCGTTGCAGTCCTAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCaaa            >  1:2127953/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      >  1:579516/1‑150 (MQ=255)
                                                                                                                                                 tGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTcc      >  1:550194/1‑150 (MQ=255)
                                                                                                                                                 tGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTcc      >  2:1179006/1‑150 (MQ=255)
                                                                                                                                                 tGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTcc      <  1:2869948/150‑1 (MQ=255)
                                                                                                                                                  gCAACGCAGTCCCAGGTTCTTTGACCAACAACTACTCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCa     <  1:3828122/150‑1 (MQ=255)
                                                                                                                                                  gCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCa     <  1:3280251/150‑1 (MQ=255)
                                                                                                                                                  gCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCa     >  2:3237198/1‑150 (MQ=255)
                                                                                                                                                  gCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCa     <  1:2679273/150‑1 (MQ=255)
                                                                                                                                                   cAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCCAACTTTCCAg    >  2:3360343/1‑150 (MQ=255)
                                                                                                                                                   cAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCAg    >  1:2109791/1‑150 (MQ=255)
                                                                                                                                                   cAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCAg    >  1:1705045/1‑150 (MQ=255)
                                                                                                                                                   cAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCAg    >  1:1704784/1‑150 (MQ=255)
                                                                                                                                                   cAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCAg    >  2:1343143/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)
                                                                                                                                                     aCGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCAGCa  <  1:4084265/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCAGTCTTTTGCCCAAGGCCCGACTACTGGGAAACTCCCCCATCAAACTTTCCAGCA  >  NC_002947/5101860‑5102158

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