Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,137,665 A→G 6.6% intergenic (‑64/‑149) 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,6650AG 6.6% 216.2 / 4.8 75intergenic (‑64/‑149)PP_4522/PP_4523LysR family transcriptional regulator/agmatinase
Reads supporting (aligned to +/- strand):  ref base A (31/39);  new base G (2/3);  total (33/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.77e-01

GCGCTGGGCGTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATCGACC  >  NC_002947/5137518‑5137808
                                                                                                                                                    |                                                                                                                                                 
gcgcTGGGCGTTGGCAAAGCCCTTGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGa                                                                                                                                                   <  1:119229/148‑1 (MQ=255)
gcgcTGGGCGTTGGCAAAGCCCTGGTGGCGCACCTCGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGa                                                                                                                                                   <  5:319284/148‑1 (MQ=255)
gcgcTGGGCGTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGa                                                                                                                                                   >  5:212622/1‑148 (MQ=255)
        cGTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTg                                                                                                                                          >  6:376278/1‑148 (MQ=255)
        cGTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGa                                                                                                                                         <  2:333187/149‑1 (MQ=255)
             gCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTaa                                                                                                                                    <  1:269582/149‑1 (MQ=255)
             gCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTaa                                                                                                                                    <  3:140232/149‑1 (MQ=255)
             gCAAAGCCATGGTGGCGCACCACGCTGATGAAGATGCTCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGTTTACACAAGGGTATGTGCGCGCTGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTaa                                                                                                                                    <  1:326117/149‑1 (MQ=255)
               aaaGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGGAATTTCTGAACTaaag                                                                                                                                  <  4:65426/149‑3 (MQ=255)
               aaaGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGt                                                                                                                                 <  5:376278/149‑1 (MQ=255)
                 aGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTGAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTAt                                                                                                                               <  2:401689/149‑1 (MQ=255)
                 aGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCCTAACTGCTCTACCCGCCTTTAGTTCAGGAATTTCTGAACTaaagg                                                                                                                                 <  1:293511/149‑4 (MQ=255)
                 aGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGGAATTTCTGAACTaaaggc                                                                                                                                <  7:275880/149‑5 (MQ=255)
                 aGCCCTGGTGGAGAACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGAAAGGTCGATGCCATGGTGCAACGGTTCCACAATGGTATGTGAGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTAt                                                                                                                               <  7:44866/149‑1 (MQ=255)
                   cccTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTAttt                                                                                                                             >  7:93145/1‑149 (MQ=255)
                   cccTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTAttt                                                                                                                             <  4:252973/149‑1 (MQ=255)
                   cccTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTAttt                                                                                                                             >  2:324876/1‑149 (MQ=255)
                        gtggCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCg                                                                                                                         >  2:318838/1‑148 (MQ=255)
                        gtggCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCg                                                                                                                         <  5:35191/148‑1 (MQ=255)
                        gtggCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCg                                                                                                                         <  1:318838/148‑1 (MQ=255)
                        gtggCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCg                                                                                                                         <  3:119060/148‑1 (MQ=255)
                           gcgcACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAg                                                                                                                     >  4:92348/1‑149 (MQ=255)
                                 cacGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGAttc                                                                                                               <  2:59650/149‑1 (MQ=255)
                                 caagCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGGCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGAttc                                                                                                               <  4:241097/146‑1 (MQ=255)
                                   cGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGAttctt                                                                                                             <  2:40428/149‑1 (MQ=255)
                                        attaaGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTccc                                                                                                          <  4:216150/144‑1 (MQ=255)
                                        atgaAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCg                                                                                                        >  1:368665/1‑149 (MQ=255)
                                          gaAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGAttcttc                                                                                                            <  4:21560/143‑1 (MQ=255)
                                          gaAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGAttcttc                                                                                                            >  3:21560/1‑143 (MQ=255)
                                          gaAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGgc                                                                                                      >  2:302353/1‑149 (MQ=255)
                                            aGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTTTGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGATATTCCTGAAGTATGTTTTTTTCGGTAGGGATTCTTCACCGgcgg                                                                                                    >  3:327048/1‑148 (MQ=255)
                                                     gcagTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAg                                                                                           <  3:92348/149‑1 (MQ=255)
                                                              ggTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTgt                                                                                  >  8:90693/1‑148 (MQ=255)
                                                               gTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTc                                                                                                                          >  8:300441/1‑108 (MQ=255)
                                                               gTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTc                                                                                                                          <  7:300441/108‑1 (MQ=255)
                                                                tagggCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGc                                                                                 <  1:94077/147‑1 (MQ=255)
                                                                tCGGGCAATGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTTCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGACAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGc                                                                                 <  1:324876/149‑1 (MQ=255)
                                                                   ggCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGa                                                                              >  3:260044/1‑148 (MQ=255)
                                                                     cAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAAt                                                                           <  2:368665/149‑1 (MQ=255)
                                                                          tCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTccc                                                                      >  1:197742/1‑149 (MQ=255)
                                                                          tCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTccc                                                                      >  4:55832/1‑149 (MQ=255)
                                                                               gCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGccc                                                                 <  3:262273/149‑1 (MQ=255)
                                                                                 cATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCccgc                                                               >  8:173236/1‑149 (MQ=255)
                                                                                  aTGGGGCAACGGGTCCACAAGCGAATGTGCGCGCAGACTAACAGCGCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGGCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCccgcc                                                              <  4:210365/149‑1 (MQ=255)
                                                                                     gTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCgg                                                            <  6:282419/148‑1 (MQ=255)
                                                                                        cAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCaaa                                                        >  5:54181/1‑149 (MQ=255)
                                                                                          aCGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAAc                                                       >  6:350442/1‑148 (MQ=255)
                                                                                                 cagaaGGGTATGTGCGCGCAGTCTAACTGCTCTACCCTCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCGCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACAGGt                                                <  6:45354/145‑1 (MQ=255)
                                                                                                  acaAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTgg                                                                                  >  5:180308/1‑113 (MQ=255)
                                                                                                  acaAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTgg                                                                                  <  6:180308/113‑1 (MQ=255)
                                                                                                  acaAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTg                                               >  2:246221/1‑148 (MQ=255)
                                                                                                    aaGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAAt                                            >  1:112755/1‑149 (MQ=255)
                                                                                                     aGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATg                                           <  8:93145/149‑1 (MQ=255)
                                                                                                     aGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATg                                           <  6:54181/149‑1 (MQ=255)
                                                                                                                  gcAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGca                              <  1:246221/149‑1 (MQ=255)
                                                                                                                      tCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCgcggcg                                 >  2:22687/1‑142 (MQ=255)
                                                                                                                      tCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCgcggcg                                 <  1:22687/142‑1 (MQ=255)
                                                                                                                      tCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGGACAgg                          <  6:212622/149‑1 (MQ=255)
                                                                                                                           aCTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTc                       >  8:375978/1‑148 (MQ=255)
                                                                                                                             tGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGaa                    >  6:229758/1‑149 (MQ=255)
                                                                                                                               ctctACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGaa                    <  4:260044/147‑1 (MQ=255)
                                                                                                                                 ctACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGt                                                <  1:218066/117‑1 (MQ=255)
                                                                                                                                 ctACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGt                                                >  2:218066/1‑117 (MQ=255)
                                                                                                                                  tACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGaacaaa                >  7:72294/1‑148 (MQ=255)
                                                                                                                                    cccGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGaacaaa                <  2:197742/146‑1 (MQ=255)
                                                                                                                                    cccGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGCTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAacac             <  5:350442/149‑1 (MQ=255)
                                                                                                                                      cGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGgc                                         >  1:336025/1‑119 (MQ=255)
                                                                                                                                      cGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGgc                                         <  2:336025/119‑1 (MQ=255)
                                                                                                                                           ttAGTTCAGGAATTTCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATc      >  5:110577/1‑149 (MQ=255)
                                                                                                                                           ttAGTTCAGGAATTTCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATc      >  1:25027/1‑149 (MQ=255)
                                                                                                                                           ttAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGTGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCg                      <  7:303836/133‑1 (MQ=255)
                                                                                                                                           ttAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGTGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCg                      >  8:303836/1‑133 (MQ=255)
                                                                                                                                           ttAGTTCACAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCACGGAGCACGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGAGGCACAGGTTCGAACAAACACAACAATc      >  2:4887/1‑149 (MQ=255)
                                                                                                                                              gTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATcgac   <  7:90693/149‑1 (MQ=255)
                                                                                                                                              gTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATcga    <  5:229758/148‑1 (MQ=255)
                                                                                                                                               ttCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACACTcgacc  >  6:50709/1‑149 (MQ=255)
                                                                                                                                               ttCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATcgacc  >  3:58918/1‑149 (MQ=255)
                                                                                                                                                tCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCTCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATcgacc  >  2:50134/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GCGCTGGGCGTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATCGACC  >  NC_002947/5137518‑5137808

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: