Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,156,328 A→G 7.6% S263G (AGC→GGC)  PP_4539 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,156,3280AG 7.6% 181.4 / 5.1 66S263G (AGC→GGC) PP_4539LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (30/31);  new base G (2/3);  total (32/34)
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.29e-01

CGCCGAGCAGCGCTGGACCGTCGGCAGCCTGGGCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACCGCGCCACTGGACGTGCCCTTCGCAGCCCCCCCACAAGCCTGAGAGGA  >  NC_002947/5156180‑5156470
                                                                                                                                                     |                                                                                                                                               
cgcCGAGCAGCGCTGGACCGGCGGCAGCCTGGGCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCa                                                                                                                                                 <  8:177461/149‑1 (MQ=255)
    gAGCAGCGCTGGACCGTCGGCAGCCTGGGCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAgccg                                                                                                                                             <  8:111169/149‑1 (MQ=255)
      gcagcGCTGGACCGTCGGCAGCCTGGGCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCGgctgcc                                                                                                                                           <  2:295602/149‑1 (MQ=255)
          cgcTGGACCGTCGGCAGCCTGGGCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGccacc                                                                                                                                        >  2:185556/1‑148 (MQ=255)
           gcTGGACCGTCGGCAGCCTGGGCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGccaccca                                                                                                                                      <  7:36632/149‑1 (MQ=255)
             tGGACCGTCGGCAGCCTGGGCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACt                                                                                                                                    >  2:194925/1‑149 (MQ=255)
               gACCGTCGGCTGCCTGGGCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACCTGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGGCATT‑GCCGCTGGATCAGGGTGGCGGCTGccaccctgatc                                                                                                                                  <  8:148395/149‑6 (MQ=255)
               gACCGTCGGCAGCCTGGGCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGt                                                                                                                                  >  4:196032/1‑149 (MQ=255)
                aCCGTCGGCAGCCTGGGCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGtt                                                                                                                                 <  1:185556/149‑1 (MQ=255)
                 ccGTCGGCAGCCTGGGCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTc                                                                                                                                <  2:210508/149‑1 (MQ=255)
                 ccGTCGGCAGCCTGGGCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTc                                                                                                                                <  2:317919/149‑1 (MQ=255)
                 ccGTCGGCAGCCTGGGCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCCGGATCAGGGTGGCAGCCGCCACCCACTGTTc                                                                                                                                <  6:123388/149‑1 (MQ=255)
                       gCAGCCTGGGCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCtt                                                                                                                          >  1:36068/1‑149 (MQ=255)
                           ccTGGGCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAt‑t                                                                                                                     <  2:201445/149‑1 (MQ=255)
                             tGGGCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCGGCTGccaccctga                                                                                                                                    >  6:52864/1‑130 (MQ=255)
                             tGGGCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCGGCTGccaccctga                                                                                                                                    <  5:52864/133‑4 (MQ=255)
                             tGGGCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCg                                                                                                                   >  4:23218/1‑149 (MQ=255)
                              gggCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGATGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGa                                                                                                                  <  1:254180/149‑1 (MQ=255)
                               ggCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATTAGCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGa                                                                                                                  <  1:268346/149‑1 (MQ=255)
                               ggCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGa                                                                                                                  <  3:196032/148‑1 (MQ=255)
                                             cTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGAccc                                                                                                   >  8:124954/1‑149 (MQ=255)
                                              ttCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCt                                                                                                  >  7:73895/1‑149 (MQ=255)
                                                cGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCTg                                                                                                 <  8:73895/148‑1 (MQ=255)
                                                   cagcagCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCtt                                                                                                                          >  6:192366/1‑121 (MQ=255)
                                                   cagcagCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCtt                                                                                                                          <  5:192366/121‑1 (MQ=255)
                                                     gcagcGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAaga                                                                                                           <  3:187523/133‑1 (MQ=255)
                                                     gcagcGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAaga                                                                                                           >  4:187523/1‑133 (MQ=255)
                                                         cGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCTGGGCCCGGCCa                                                                                       >  2:84648/1‑149 (MQ=255)
                                                          ggCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCTGGGCCCGGCCAc                                                                                      <  7:124954/149‑1 (MQ=255)
                                                          ggCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCTGGGCCCGGCCAc                                                                                      <  4:296974/149‑1 (MQ=255)
                                                           gCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGccaccca                                                                                                                                      <  1:270831/101‑1 (MQ=255)
                                                           gCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGccaccca                                                                                                                                      >  2:270831/1‑101 (MQ=255)
                                                            ccTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGc                                                                                    <  2:36068/149‑1 (MQ=255)
                                                             cTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCATGACATA‑GCCGCTCGATCAGGGTGGCAGCCGCCGCCGACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCTGGGCCCGGCCACg                                                                                     >  5:80485/1‑147 (MQ=255)
                                                                ggCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCCCTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGAt                                                                                >  7:327809/1‑149 (MQ=255)
                                                                 gCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGAtt                                                                               <  1:194925/149‑1 (MQ=255)
                                                                   tttGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATtct                                                                             >  3:129563/1‑149 (MQ=255)
                                                                      gCCTGGCTGCCCCGTCACATGATCGCGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCCGGGGGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCAt                                                                          >  5:250245/1‑149 (MQ=255)
                                                                        cTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGAcc                                                                                                    <  6:225222/121‑1 (MQ=255)
                                                                        cTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGAcc                                                                                                    >  5:225222/1‑121 (MQ=255)
                                                                        cTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGc                                                                                    >  7:121793/1‑137 (MQ=255)
                                                                        cTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGc                                                                                    <  8:121793/137‑1 (MQ=255)
                                                                           gCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTATCTCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGAc                                                                      <  1:84648/149‑1 (MQ=255)
                                                                                cccGTCACATGATCGAGCGCGAGCTGCGCGAAGCCCTCCTCAAGCCATT‑GCCGCTGGATCAGCGAGGCAGCCGCCACCCACTGTTCTACCTTAAT‑TCCAGCAAAGAGAAGACCCTGGGCCAGGCCACGCAGATTCTCAACGACctgctc                                                                >  6:23956/1‑148 (MQ=255)
                                                                                    tCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCt                                                                                                  <  4:127550/111‑1 (MQ=255)
                                                                                    tCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCt                                                                                                  >  3:127550/1‑111 (MQ=255)
                                                                                    tCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGACCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCa                                                            <  3:334031/149‑1 (MQ=255)
                                                                                           aTCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGa                                                     <  5:23956/149‑1 (MQ=255)
                                                                                              gAGCGCGAGCTGTGAGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGacac                                                  <  1:301350/149‑1 (MQ=255)
                                                                                                  gcgAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACcgcg                                              >  5:162373/1‑149 (MQ=255)
                                                                                                   cgAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACCgcgc                                             <  6:68969/149‑1 (MQ=255)
                                                                                                             gAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACCGCGCCACTGGACGt                                   >  5:13599/1‑149 (MQ=255)
                                                                                                                    tGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACCGCGCCACTGGACGTGCCCTTc                            >  3:118251/1‑149 (MQ=255)
                                                                                                                      cTCAACCCATT‑GCCGCTGGATCAGGGTCGCAGCCCCCACCCACTGATCTACCTTTAT‑TAGAGCAAAGAGAACACCCAGGGCCCCGCCACGCAGAGTCTCATCGACCTGCTGCGCAATTTCGACACCGCCCACCAGGACGAGCCCTTCg                           >  4:127883/1‑148 (MQ=255)
                                                                                                                       tCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCTGGGCCCGGCCACTCAGATTCTCATCGACCTGCTGCGCAATTTCGACACAGCGCCACTGGACGAGCCCATCg                           >  4:33928/1‑147 (MQ=255)
                                                                                                                         aaGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACCGCGCCACTGGACGTGCCCTTCGCAg                        <  6:162373/148‑1 (MQ=255)
                                                                                                                            caact‑GCCTCTGGATCACTGTGGCAGCCGCCACCCACTGAGCTCCCTTTAT‑TCGAGCAAAGAGAAGACCCTGGGCCCGTCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACCGCGCCACTGGACGTGCCCTTCGCAGcccc                    <  6:250245/145‑1 (MQ=255)
                                                                                                                             cATT‑GCCGCTGTAGCAGGGTTGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGATCACAGAGACCCCCCTGGGCCCGGCCACGCAGATTCTGATCGACCTGCTGCGCAATTTCGACACCGCGCCACTGTTCGTTCCCTTCGCAGccccc                   <  6:80485/149‑1 (MQ=255)
                                                                                                                              aTT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCTGGGCCCTGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACCGCGCCACTGGACGTGCCCTTCGCAGcccccc                  >  6:114530/1‑149 (MQ=255)
                                                                                                                                   ccGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACCGCGCCACTGGACGTGCCCTTCGCAGCCCCCCcac               >  2:22614/1‑148 (MQ=255)
                                                                                                                                   ccGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACCGCGCCACTGGACGTGCCCTTCGCAGCCCCCCcac               >  5:110258/1‑148 (MQ=255)
                                                                                                                                           tCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACCGCGCCACTGGACGTGCCCTTCGCAGCCCCCCCACAAGCCt         >  3:254538/1‑146 (MQ=255)
                                                                                                                                           tCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACCGCGCCACTGGACGTGCCCTTCGCAGCCCCCCCACAAGCCt         >  4:9094/1‑146 (MQ=255)
                                                                                                                                           agtgggTGGCGGCTGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACCGCGCCACTGGACGTGCCCTTCGCAGCCCCCCCACAAGCCt         >  4:219201/4‑146 (MQ=255)
                                                                                                                                            gaGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCTTGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACCGCGCCACTGGTCGTGCCCTTCGCAGCCCCCCCACAAGCCTgag      <  3:127883/147‑1 (MQ=255)
                                                                                                                                             aGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACCGCGCCACTGGACGTGCCCTTCGCAGCCCCCCCACAAGCCTgaga     <  3:33928/148‑1 (MQ=255)
                                                                                                                                               ggTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACCGCGCCACTGGACGTGCCCTTCGCAGCCCCCCCACAAGCCTGAGAGGa  >  5:251771/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
CGCCGAGCAGCGCTGGACCGTCGGCAGCCTGGGCACCGCCGCCACCTTCGTCAGCAGCGGCCTGGGCTTTGCCTGGCTGCCCCGTCACATGATCGAGCGCGAGCTGCGCGAAGGCGTGCTCAAGCCATT‑GCCGCTGGATCAGGGTGGCAGCCGCCACCCACTGTTCTACCTTTAT‑TCGAGCAAAGAGAAGACCCTGGGCCCGGCCACGCAGATTCTCATCGACCTGCTGCGCAATTTCGACACCGCGCCACTGGACGTGCCCTTCGCAGCCCCCCCACAAGCCTGAGAGGA  >  NC_002947/5156180‑5156470

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