Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,137,670 C→T 9.4% intergenic (‑69/‑144) PP_4522 ← / → PP_4523 LysR family transcriptional regulator/agmatinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,137,6700CT 9.4% 204.4 / 8.4 74intergenic (‑69/‑144)PP_4522/PP_4523LysR family transcriptional regulator/agmatinase
Reads supporting (aligned to +/- strand):  ref base C (26/41);  new base T (3/4);  total (29/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-02

GGGCGTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATCGACCGA  >  NC_002947/5137523‑5137810
                                                                                                                                                   |                                                                                                                                               
gggCGTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTcc                                                                                                                                                >  1:133586/1‑149 (MQ=255)
 ggCGTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTcc                                                                                                                                                <  1:121569/148‑1 (MQ=255)
 ggCGTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCt                                                                                                                                               <  8:363012/149‑1 (MQ=255)
 ggCGTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCt                                                                                                                                               <  2:140067/149‑1 (MQ=255)
 ggCGTGTGCAACGCCCTGGTGGCGCACCACGCTGATGTAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTc                                                                                                                                                 <  8:296914/147‑1 (MQ=255)
    gTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGaa                                                                                                                                            <  7:374327/149‑1 (MQ=255)
     ttGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAAc                                                                                                                                           >  3:409135/1‑149 (MQ=255)
       ggCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGGAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACt                                                                                                                                          <  5:321040/148‑1 (MQ=255)
        gCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTaa                                                                                                                                        <  8:166112/149‑1 (MQ=255)
         cAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGTCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTaa                                                                                                                                        <  3:208189/148‑1 (MQ=255)
          aaaGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGGAATTTCTGAACTaaag                                                                                                                                      <  7:112356/149‑3 (MQ=255)
            aGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGGAATTTCTGAACTaaaggc                                                                                                                                    <  3:29149/149‑5 (MQ=255)
              cccTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGGAATTTCTGAACTaaaggcgg                                                                                                                                  <  7:394061/149‑7 (MQ=255)
              cccTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTAttt                                                                                                                                  <  6:61896/149‑1 (MQ=255)
               ccTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGGAATTTCTGAACTaaaggcggg                                                                                                                                 <  8:335333/149‑8 (MQ=255)
                cTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGt                                                                                                                                <  5:416534/149‑1 (MQ=255)
                cTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGt                                                                                                                                <  8:214920/149‑1 (MQ=255)
                      gcgcACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGAGCGCGCAGTCTAACTGCTCTACCAGCCTTTAGTTCAGAAATTCATGAACTAAGTATTTGTCGGTAg                                                                                                                         >  7:119352/1‑149 (MQ=255)
                            cacGCTGATGAAGATGCTCAGCAGTTTCAGGTCGTGCAAGGTCGATGCCATGGTGCATCGGTTCCACAAGTGTATGTGCGCGCAGTCTAACTGCTGTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GAttc                                                                                                                 <  7:260190/149‑1 (MQ=255)
                               gatgaagaAGATGCGCAGCAGTTTCAGGTCGCGCAAGGAAGACGCCATGGTGCAACGGTGCCACATGGGTATGTGCGCGCAGTCAAACTGCTGTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GAttcttc                                                                                                              <  8:407614/147‑1 (MQ=255)
                               gCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GAttcttc                                                                                                              <  3:47087/149‑1 (MQ=255)
                                   atgaAGATGCGCAGCAGTTTCATGTCGGGCAAGGTCGATGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTcccc                                                                                                           <  4:285704/148‑1 (MQ=255)
                                         aTGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGAAGC‑‑GATTCTTCCCCGGcgctc                                                                                                    <  7:339354/149‑1 (MQ=255)
                                            cgcAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGcgctcgc                                                                                                  <  1:19679/148‑1 (MQ=255)
                                            cgcAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGcgctcgc                                                                                                  <  1:393018/148‑1 (MQ=255)
                                             gcagcaGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCtt                                                                                                >  4:219581/1‑149 (MQ=255)
                                               agcagTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑Ga                                                                                                                    <  6:73829/127‑1 (MQ=255)
                                               agcagTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑Ga                                                                                                                    >  5:73829/1‑127 (MQ=255)
                                                        aGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTg                                                                                     >  2:431457/1‑149 (MQ=255)
                                                          gTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTACCTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGGCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGc                                                                                   <  3:416735/149‑1 (MQ=255)
                                                              ggCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGa                                                                                >  5:340660/1‑148 (MQ=255)
                                                               gCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGaaa                                                                              <  8:272177/149‑1 (MQ=255)
                                                                  aGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGt                                                                           <  7:312791/149‑1 (MQ=255)
                                                                    gTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTcc                                                                         <  3:219581/149‑1 (MQ=255)
                                                                     tCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTccc                                                                        >  5:64773/1‑149 (MQ=255)
                                                                     tCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTccc                                                                        <  2:133586/149‑1 (MQ=255)
                                                                       gAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTg                                                                      <  4:337943/149‑1 (MQ=255)
                                                                         ggCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGcc                                                                    >  7:160944/1‑149 (MQ=255)
                                                                            cATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGTGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCccgc                                                                 >  8:213578/1‑149 (MQ=255)
                                                                            cATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCccgc                                                                 >  5:28113/1‑149 (MQ=255)
                                                                             aTGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCccgcc                                                                <  5:403802/149‑1 (MQ=255)
                                                                              tGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCccgccg                                                               <  6:340660/149‑1 (MQ=255)
                                                                                   cAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCaaa                                                          >  2:90471/1‑149 (MQ=255)
                                                                                    aaCGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAAc                                                         <  2:205095/149‑1 (MQ=255)
                                                                                     aCGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAAc                                                         >  4:308729/1‑148 (MQ=255)
                                                                                       ggTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACca                                                       >  1:416497/1‑148 (MQ=255)
                                                                                            cacaAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGt                                                  <  3:4791/148‑1 (MQ=255)
                                                                                             acaAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGAAANGCGGGATTAGAAACCGGCGCACGCTTGATCAGACTTGGCCGAAATTTCCCTGCCCGCCGGCAAACCACACg                                                   >  1:141041/1‑149 (MQ=255)
                                                                                                      tgtgCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGgcgcg                                        <  6:64773/148‑1 (MQ=255)
                                                                                                           gcgcAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGcggcgg                                  <  1:431457/149‑1 (MQ=255)
                                                                                                               aGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGca                                >  1:185676/1‑147 (MQ=255)
                                                                                                                 tCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAgg                            <  8:160944/149‑1 (MQ=255)
                                                                                                                   tAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGtt                          >  3:53912/1‑149 (MQ=255)
                                                                                                                     aCTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGtt                          <  4:409135/147‑1 (MQ=255)
                                                                                                                      cTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGa                       <  4:53912/149‑1 (MQ=255)
                                                                                                                       tGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGaa                      >  2:76735/1‑149 (MQ=255)
                                                                                                                        gCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGaa                      >  7:160286/1‑148 (MQ=255)
                                                                                                                          tctACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGgcg                                          <  8:244306/126‑1 (MQ=255)
                                                                                                                          tctACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGgcg                                          >  7:244306/1‑126 (MQ=255)
                                                                                                                             aCCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCGCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGaacaaaca                >  1:226604/1‑149 (MQ=255)
                                                                                                                             aCCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGaacaaaca                >  2:307725/1‑149 (MQ=255)
                                                                                                                             aCCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGaacaaaca                >  4:217836/1‑149 (MQ=255)
                                                                                                                              cccGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGaacaaa                  <  3:308729/146‑1 (MQ=255)
                                                                                                                              cccGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGaacaaa                  <  2:185676/146‑1 (MQ=255)
                                                                                                                               ccGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAacac               >  1:237054/1‑148 (MQ=255)
                                                                                                                               ccGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAAcaca              >  4:208356/1‑149 (MQ=255)
                                                                                                                               aaatacttAGTTCAGGAATTTCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAacac               >  7:317437/7‑148 (MQ=255)
                                                                                                                                    cttAGTTCAGGAATTTCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAAt         >  5:92591/2‑149 (MQ=255)
                                                                                                                                     ttAGTTCAGGAATTTCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATc        >  5:204624/1‑149 (MQ=255)
                                                                                                                                      tAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAAt         <  1:90471/147‑1 (MQ=255)
                                                                                                                                      tAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATc        <  7:106319/148‑1 (MQ=255)
                                                                                                                                        gTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCGCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATcga      <  2:226604/148‑1 (MQ=255)
                                                                                                                                         ttCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGtt                          >  3:309693/1‑127 (MQ=255)
                                                                                                                                         ttCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGtt                          <  4:309693/127‑1 (MQ=255)
                                                                                                                                         ttCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATcgacc    <  1:115320/149‑1 (MQ=255)
                                                                                                                                           cAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACCCAACACTcgaccga  >  7:215145/1‑149 (MQ=255)
                                                                                                                                            aGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATcgaccga  <  2:141041/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GGGCGTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGC‑‑GATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATCGACCGA  >  NC_002947/5137523‑5137810

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: