Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,382,269 A→T 10.0% E124D (GAA→GAT smpB → SsrA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,382,2690AT10.0% 265.3 / 15.6 97E124D (GAA→GATsmpBSsrA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (42/45);  new base T (5/5);  total (47/50)
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.88e-01

GCACCCACGTCATCGCCGACCCTACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGTTCCGGCCTCTTCGCGGGCAA  >  NC_002947/5382121‑5382416
                                                                                                                                                    |                                                                                                                                                   
gCACCCACGTCATCGCCGACCCTACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGaa                                                                                                                                                     <  1:144408/149‑1 (MQ=255)
      aCGTCATCGCCGACCCTACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCa                                                                                                                                               <  2:48582/149‑1 (MQ=255)
        gTCATCGCCGACCCTACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACt                                                                                                                                             <  6:63289/149‑1 (MQ=255)
        gTCATCGCCGACCCTACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACt                                                                                                                                             >  7:43049/1‑149 (MQ=255)
        gTCATCGCCGACCCTACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACAACTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACt                                                                                                                                             <  7:354573/149‑1 (MQ=255)
         tCATCGCCGACCCTACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTg                                                                                                                                            >  7:289085/1‑149 (MQ=255)
         tCATCGCCGACCCTACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTg                                                                                                                                            <  2:699/149‑1 (MQ=255)
          cATCGCCGACCCTACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTgg                                                                                                                                           >  2:139597/1‑149 (MQ=255)
           aTCGCCGACCCTACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTgg                                                                                                                                           <  1:139597/148‑1 (MQ=255)
              gCCGACCCTACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGATTTCGCACttgatc                                                                                                                                        <  3:329822/148‑6 (MQ=255)
               ccGACCCTACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGATTTCGCACttgatcag                                                                                                                                      <  4:364464/149‑8 (MQ=255)
               ccGACCCTACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTgggcaag                                                                                                                                      <  8:285521/149‑1 (MQ=255)
                 gACCCTACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTgggcaaggg                                                                                                                                    >  2:74634/1‑149 (MQ=255)
                      tACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGggcaagggcaag                                                                                                                                >  6:220252/1‑148 (MQ=255)
                        cacGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCaagaag                                                                                                                             <  1:5358/149‑1 (MQ=255)
                        cacGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCaagaag                                                                                                                             <  8:46394/149‑1 (MQ=255)
                        cacGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCaagaa                                                                                                                              <  4:57312/148‑1 (MQ=255)
                        cacGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCaagaa                                                                                                                              <  5:220252/148‑1 (MQ=255)
                        cacGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACATGATCAAGTGCGAAATCGCACTGGGCAAGGGCaagaag                                                                                                                             >  8:348922/1‑149 (MQ=255)
                              cccGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTc                                                                                                                       <  6:260573/149‑1 (MQ=255)
                                  cAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACa                                                                                                                   <  1:74634/149‑1 (MQ=255)
                                   aaGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACaa                                                                                                                  <  8:255543/149‑1 (MQ=255)
                                     gctgctgctGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAgc                                                                                                                >  6:249474/1‑149 (MQ=255)
                                         ctgctgAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGaa                                                                                                                          <  6:236915/135‑1 (MQ=255)
                                         ctgctgAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGa                                                                                                            <  7:348922/149‑1 (MQ=255)
                                              gAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCACGCACCTGATCAAGTGCGAAGTCGCGCTGGGCAGGGGGAGGAGGGAGTGCGGGAGGCGCGaggcgg                                                                                                       >  4:267525/1‑144 (MQ=255)
                                                aCAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCAt                                                                                                      <  1:322131/148‑1 (MQ=255)
                                                 cAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGTCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGTTCAAGAAGGAATTCGACAAGCGCGACACCATgc                                                                                                    <  6:201521/149‑1 (MQ=255)
                                                 cAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATgc                                                                                                    >  6:368868/1‑149 (MQ=255)
                                                 cAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATgc                                                                                                    <  1:136090/149‑1 (MQ=255)
                                                    gcgcGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGggcaagggcaag                                                                                                                                <  8:61623/118‑1 (MQ=255)
                                                    gcgcGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGggcaagggcaag                                                                                                                                >  7:61623/1‑118 (MQ=255)
                                                       cgAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGac                                                                                                           >  3:367339/1‑136 (MQ=255)
                                                       cgAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGac                                                                                                           <  4:367339/136‑1 (MQ=255)
                                                           cTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGa                                                                                          >  2:274350/1‑149 (MQ=255)
                                                              gAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGATATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTc                                                                                       >  6:198539/1‑149 (MQ=255)
                                                              gAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTc                                                                                       >  7:75038/1‑149 (MQ=255)
                                                                 cgcgTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCaa                                                                                    <  2:343887/149‑1 (MQ=255)
                                                                  gcgTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGATATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAAc                                                                                   <  5:198539/149‑1 (MQ=255)
                                                                    gTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACcg                                                                                 >  6:62611/1‑149 (MQ=255)
                                                                    gTGGAAGCGGACGTCGCGCAAAAAGGCTACACCTGTGTGGCTTTCGCGCTGTACTGGACCAACCACATGATCAACTCCGAAATCGCACTCGACAAGGGCAAGAACGAATTCGACAAGCGACACCCCCTGCGCGAACGCGATTCCAAcca                                                                                 >  8:180029/1‑148 (MQ=255)
                                                                       gAAGCGGCCGTGGCGCAAAAAGGCTACTCCTGTGTGTCTTTGGCGCGGTACTGCAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGCGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGa                                                                              <  8:123943/149‑1 (MQ=255)
                                                                       gAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACcgcg                                                                               <  8:43049/148‑1 (MQ=255)
                                                                         aGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCaa                                                                                    >  7:59386/1‑141 (MQ=255)
                                                                         aGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCaa                                                                                    <  8:59386/141‑1 (MQ=255)
                                                                         aGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGc                                                                            <  7:180029/149‑1 (MQ=255)
                                                                            ggCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGa                                                                                          <  8:135862/132‑1 (MQ=255)
                                                                            ggCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGa                                                                                          >  7:135862/1‑132 (MQ=255)
                                                                            ggCCGTGGCGCAAAAAGGCTACACCCGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTg                                                                          >  7:201210/1‑148 (MQ=255)
                                                                               cGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGc                                                                      >  3:299244/1‑149 (MQ=255)
                                                                                   gcgcAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCggg                                                                   <  8:75038/148‑1 (MQ=255)
                                                                                     gcAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCtg                                                                >  6:382416/1‑149 (MQ=255)
                                                                                      cAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCtgt                                                               >  6:305958/1‑149 (MQ=255)
                                                                                       aaaaaGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCggg                                                                   >  8:300732/1‑144 (MQ=255)
                                                                                       aaaaaGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCggg                                                                   <  7:300732/144‑1 (MQ=255)
                                                                                        aaaaGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGc                                                             >  5:6037/1‑149 (MQ=255)
                                                                                        aaaaGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGc                                                             >  6:331715/1‑149 (MQ=255)
                                                                                         aaaGGCTTCACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATGCCAACCGCGAGCTGCAGCGGGCTGTGCg                                                            >  6:271031/1‑149 (MQ=255)
                                                                                                 cacCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAgg                                                    <  4:268278/149‑1 (MQ=255)
                                                                                                 cacCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAgg                                                    >  6:216929/1‑149 (MQ=255)
                                                                                                    cTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGc                                                             >  6:19371/1‑137 (MQ=255)
                                                                                                    cTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGc                                                             <  5:19371/137‑1 (MQ=255)
                                                                                                    cTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGc                                                  >  3:267689/1‑148 (MQ=255)
                                                                                                    cTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGc                                                  >  4:264601/1‑148 (MQ=255)
                                                                                                       tgtgGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAg                                               >  4:267074/1‑148 (MQ=255)
                                                                                                          ggCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAgg                                                    <  1:274539/140‑1 (MQ=255)
                                                                                                          ggCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAgg                                                    >  2:274539/1‑140 (MQ=255)
                                                                                                           gCTTTGGCGCTGTACTGGAGCAAGCACCTTATCATGTGCGAAATGGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAaga                                          <  5:249474/149‑1 (MQ=255)
                                                                                                           gCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAaga                                          <  5:368868/149‑1 (MQ=255)
                                                                                                            cTTTGGCGCTGTACTGGAGCAAGCACCTTATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAagag                                         <  8:289085/149‑1 (MQ=255)
                                                                                                             tttGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAagag                                         >  3:26176/1‑148 (MQ=255)
                                                                                                                   gcTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAactac                                  >  7:118822/1‑149 (MQ=255)
                                                                                                                     atgtcTGGAGCAAGCACGTGATCAAGTGCGAAATCGCACTTGGCAAGGGAAAGAAGTAATTCGACAAGAGCGACACCTTGCGCGAACGCGATTCCAACGGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTActc                                <  8:271486/145‑1 (MQ=255)
                                                                                                                          tGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGAAAAGCGCGACACCATGCGCGAACGCGATTCCAACAGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTActc                                >  7:271486/1‑144 (MQ=255)
                                                                                                                            gAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAg                                                     <  8:346812/121‑1 (MQ=255)
                                                                                                                            gAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAg                                                     >  7:346812/1‑121 (MQ=255)
                                                                                                                            gAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAAcc                         >  2:335355/1‑149 (MQ=255)
                                                                                                                            gAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAAcc                         >  2:44553/1‑149 (MQ=255)
                                                                                                                            gAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAAc                          <  4:299244/148‑1 (MQ=255)
                                                                                                                               caagcaCCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGt                      >  3:182359/1‑149 (MQ=255)
                                                                                                                                  gcaCCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGTTcc                   <  8:148369/149‑1 (MQ=255)
                                                                                                                                  gcaCCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGTTcc                   <  3:267074/149‑1 (MQ=255)
                                                                                                                                   caCCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGTTCCg                  <  7:32981/149‑1 (MQ=255)
                                                                                                                                   caCCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGTTCCg                  <  1:274350/149‑1 (MQ=255)
                                                                                                                                    aCCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGTTCCgg                 >  4:5992/1‑149 (MQ=255)
                                                                                                                                     ccTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGTTCCGGc                <  8:118822/149‑1 (MQ=255)
                                                                                                                                      cttgcccaGTGCGATTTCGCACTTGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGTTCCGGcc               >  3:365786/8‑149 (MQ=255)
                                                                                                                                      cTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCg                                                                     <  2:374583/95‑1 (MQ=255)
                                                                                                                                      cTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCg                                                                     >  1:374583/1‑95 (MQ=255)
                                                                                                                                         aTCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGTTCCGGCct              <  5:305958/147‑1 (MQ=255)
                                                                                                                                          cccaGTGCGATTTCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGTTCCGGCCTCtt           >  3:57893/4‑149 (MQ=255)
                                                                                                                                           ccaGTGCGATTTCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGTTCCGGCCTCTTc          >  4:368667/3‑149 (MQ=255)
                                                                                                                                             aGTGCGATTTCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGTTCCGGc                <  8:159566/141‑1 (MQ=255)
                                                                                                                                             aGTGCGATTTCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGTTCCGGc                >  7:159566/1‑141 (MQ=255)
                                                                                                                                              gTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGTTCCGGCCTCTTcgcg       <  4:267689/149‑1 (MQ=255)
                                                                                                                                               tGCGATTTCGCACTTGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGTTCCGGCCTCTTCGCgg      <  4:365786/149‑1 (MQ=255)
                                                                                                                                                gCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTaa                                      >  6:40747/1‑116 (MQ=255)
                                                                                                                                                  gAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGTTCCGGCCTCTTCGCGGGc    >  6:96042/1‑148 (MQ=255)
                                                                                                                                                    aaTCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGTTCCGGCCTCTTCGCGGGCaa  <  4:26176/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCACCCACGTCATCGCCGACCCTACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGTTCCGGCCTCTTCGCGGGCAA  >  NC_002947/5382121‑5382416

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