Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,944,743 A→G 6.5% intergenic (+103/+695) PP_5673 → / ← flhB D‑alanine‑‑D‑alanine ligase/flagellin export apparatus substrate specificity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,944,7430AG 6.5% 273.2 / 6.2 92intergenic (+103/+695)PP_5673/flhBD‑alanine‑‑D‑alanine ligase/flagellin export apparatus substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base A (39/47);  new base G (3/3);  total (42/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 = 8.04e-01

TTATCTTAAACAACTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATACACAGACCAGC  >  NC_002947/4944603‑4944881
                                                                                                                                             |                                                                                                                                          
ttATCTTAAACAACTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTc                                                                                                                                    >  1:337950/1‑149 (MQ=255)
        aacaacTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGTATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGa                                                                                                                            <  6:41414/149‑1 (MQ=255)
        aacaacTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTg                                                                                                                             <  1:110854/148‑1 (MQ=255)
        aacaacTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGa                                                                                                                            <  6:184462/149‑1 (MQ=255)
        aacaacTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGa                                                                                                                            <  3:117770/149‑1 (MQ=255)
          caacTTCGATCCTACCCAATACTGCTGACCGAACAACAACTCTAATTCCTACCATCAGGCTATTGGATGTCGTAGCGGTTTTTCGGATTGTAATTTATTGATCACCCAATTGAGATTAACCGCCGGTTGAGCGGCGCTTAATCtcaatt                                                                                                                          <  8:288802/149‑6 (MQ=255)
           aacTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTGATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAg                                                                                                                           >  1:294235/1‑148 (MQ=255)
                  aTCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAg                                                                                                                  <  1:25163/149‑1 (MQ=255)
                      tAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAAc                                                                                                              >  4:44459/1‑149 (MQ=255)
                       aGCCCATACTGCTGACCCAAGAACAAGTATAATTCGTACCATCAGGCTAATGGATGTCGTTCGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACg                                                                                                             <  8:241954/149‑1 (MQ=255)
                        gCCAATACTTCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGATATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAAc                                                                                                              <  2:313113/147‑1 (MQ=255)
                        gCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGTTTGAGCGGCGCttaa                                                                                                                                  <  4:153897/127‑4 (MQ=255)
                        gCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGTTTGAGCGGCGCttaa                                                                                                                                  >  3:153897/1‑124 (MQ=255)
                        gCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACgg                                                                                                            >  8:221236/1‑149 (MQ=255)
                        gCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACAGCGCTCGATCTTGAGCAATCCCAGAAACgg                                                                                                            <  4:79678/149‑1 (MQ=255)
                         ccAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGt                                                                                                           <  6:314858/149‑1 (MQ=255)
                            aTACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGTTTGAGCGGCGCTTAATCtcaa                                                                                                                            <  2:178176/129‑4 (MQ=255)
                            aTACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGTTTGAGCGGCGCTTAATCtcaa                                                                                                                            >  1:178176/1‑126 (MQ=255)
                              aCTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATTGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGa                                                                                                      <  3:280934/149‑1 (MQ=255)
                                tgctgACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGaa                                                                                                     >  1:48930/1‑148 (MQ=255)
                                     aCCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAg                                                                                                                           >  6:304321/1‑121 (MQ=255)
                                     aCCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACttt                                                                                                <  8:114034/148‑1 (MQ=255)
                                      ccGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTc                                                                                               >  1:203410/1‑148 (MQ=255)
                                       cGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTc                                                                                                          <  6:110979/136‑1 (MQ=255)
                                       cGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTc                                                                                                          >  5:110979/1‑136 (MQ=255)
                                            aacaaGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGt                                                                                         >  7:116370/1‑148 (MQ=255)
                                             acaaGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCa                                                                                       <  3:225703/149‑1 (MQ=255)
                                               aaGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCAcg                                                                                     >  3:332136/1‑149 (MQ=255)
                                               aaGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCAcg                                                                                     >  5:198442/1‑149 (MQ=255)
                                                 gTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCAcgcc                                                                                   <  2:203410/149‑1 (MQ=255)
                                                      aTTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAAc                                                                                                   <  2:361151/128‑1 (MQ=255)
                                                      aTTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAAc                                                                                                   >  1:361151/1‑128 (MQ=255)
                                                       ttCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTg                                                                                                                             >  2:250220/1‑101 (MQ=255)
                                                       ttCGTACAATCATGCTATTGGATGTGGTTGCTGTTTTTCGCATGTTAATTTTTTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAACCGGTCAAGAAACTTTCAACTGTCACGCCGCAAAt                                                                             <  2:360909/149‑1 (MQ=255)
                                                        tCGTACAATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGAGCTCGATCTTGAGCAATCACAGAAACGGTCAAGAAACTTTCAAATTTCAAGCCGCAAATc                                                                            >  5:25087/1‑149 (MQ=255)
                                                           tACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCcgg                                                                         <  6:361572/149‑1 (MQ=255)
                                                           tACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCACTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCcgg                                                                         >  5:145672/1‑149 (MQ=255)
                                                                 gaGGCTTTTTGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTTATCACACAATTGAGATTAAGCGCATCTCAAACGGCGCTGGATCTTGAGCAATCCCAGAAACGTTCAAGAAACTTTCTACTGTCAAGCCGCAAATCCGGGGGaa                                                                    <  2:48930/147‑1 (MQ=255)
                                                                 cAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAAt                                                                   >  3:353380/1‑149 (MQ=255)
                                                                  aGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATc                                                                  <  2:337950/149‑1 (MQ=255)
                                                                   ggCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCg                                                                 >  6:328352/1‑149 (MQ=255)
                                                                     cTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGa                                                               <  5:328352/149‑1 (MQ=255)
                                                                      tATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGa                                                               >  3:133119/1‑148 (MQ=255)
                                                                        aTTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGa                                                               >  5:217049/1‑147 (MQ=255)
                                                                         ttGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAg                                                            >  2:14258/1‑149 (MQ=255)
                                                                           ggATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGc                                                           <  7:179594/148‑1 (MQ=255)
                                                                            gATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGa                                                                                                                            >  7:130314/1‑82 (MQ=255)
                                                                            gATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGa                                                                                                                            <  8:130314/82‑1 (MQ=255)
                                                                            gATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAg                                                         >  3:328318/1‑149 (MQ=255)
                                                                                 cGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATgg                                                    >  8:157352/1‑149 (MQ=255)
                                                                                   ttGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTc                                                  <  2:55805/149‑1 (MQ=255)
                                                                                    tGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGAAAAGCAGGATGGTCg                                                 >  1:301692/1‑149 (MQ=255)
                                                                                      gtgtTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCggg                                               <  8:116370/149‑1 (MQ=255)
                                                                                      gtgtTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCgg                                                <  4:133119/148‑1 (MQ=255)
                                                                                      gtgtTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCgg                                                <  4:353380/148‑1 (MQ=255)
                                                                                      gtgtTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCgg                                                >  8:29186/1‑148 (MQ=255)
                                                                                      gtgtTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCgg                                                <  1:146680/148‑1 (MQ=255)
                                                                                               ggATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGc                                       >  3:86417/1‑148 (MQ=255)
                                                                                               ggATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGc                                       >  3:91622/1‑148 (MQ=255)
                                                                                                    gTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACt                                                                                           <  1:301398/91‑1 (MQ=255)
                                                                                                    gTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACt                                                                                           >  2:301398/1‑91 (MQ=255)
                                                                                                     tAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTc                                 >  1:374135/1‑148 (MQ=255)
                                                                                                       aTTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGa                               <  4:86417/148‑1 (MQ=255)
                                                                                                       aTTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGa                               <  4:91622/148‑1 (MQ=255)
                                                                                                        tttATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGAAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGTATGGTCGGGAACGCGGCAAAGACGAGc                             >  7:334499/1‑149 (MQ=255)
                                                                                                            ttGATCACCCAAGTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGgc                         <  8:295777/149‑1 (MQ=255)
                                                                                                               aTCACCCTATTGATATTAAGCGCCGCTCTAACGGCGCTCGATCTTGATCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACTCCTCAAATCCGGCGGAATCGGAAAGCAGGATGGGCGGGAACGCGGCAAAGGCGAGCGGgcgc                       <  8:334499/148‑1 (MQ=255)
                                                                                                                tCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCaa                     <  7:29186/149‑1 (MQ=255)
                                                                                                                 cACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAg                    >  1:285318/1‑149 (MQ=255)
                                                                                                                   cccAATTGAGATTAAGCGCCTCTCAAACGGCGCTCGATCTTGAGCAATCCCATAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGTAACGCGGCAAAGTCGAGCGGGCGCAAGGc                  <  6:25087/149‑1 (MQ=255)
                                                                                                                    ccAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCaaa                                    <  3:251077/130‑1 (MQ=255)
                                                                                                                    ccAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCaaa                                    >  8:288916/1‑130 (MQ=255)
                                                                                                                    ccAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCaaa                                    >  4:251077/1‑130 (MQ=255)
                                                                                                                    ccAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCaaa                                    <  7:288916/130‑1 (MQ=255)
                                                                                                                       aTTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGta               >  7:300086/1‑148 (MQ=255)
                                                                                                                       aTTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGAAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCCGCAAAGTCGAGCGGGCGCAAGGCGta               >  4:343220/1‑148 (MQ=255)
                                                                                                                        ttGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGtata             >  6:30045/1‑149 (MQ=255)
                                                                                                                          aagaTCAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATaca           >  2:314111/2‑149 (MQ=255)
                                                                                                                           ccaTTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGGCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATaca           <  8:59262/146‑1 (MQ=255)
                                                                                                                           agaTTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATacac          <  7:221236/149‑1 (MQ=255)
                                                                                                                           agaTTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATaca           <  2:285318/148‑1 (MQ=255)
                                                                                                                            gatcgaGCGCCGTTTGAGCGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATaca           >  5:138500/6‑147 (MQ=38)
                                                                                                                            gaTTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATc                                                                            <  7:275738/82‑1 (MQ=255)
                                                                                                                            gaTTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATc                                                                            >  8:275738/1‑82 (MQ=255)
                                                                                                                            gaTTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATAcaca         <  4:328318/149‑1 (MQ=255)
                                                                                                                                aaGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATACACAGAc      <  4:332136/148‑1 (MQ=255)
                                                                                                                                aaGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATACACAGAc      <  2:374135/148‑1 (MQ=255)
                                                                                                                                  gcgcCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATACACAGAcc     >  4:240790/1‑147 (MQ=255)
                                                                                                                                  gcgcCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATACACAGAcc     <  3:240790/147‑1 (MQ=255)
                                                                                                                                  gcgcCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATACACAGACcag   <  1:229267/149‑1 (MQ=255)
                                                                                                                                   cgccgcTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATACACAGACcagc  <  6:198442/149‑1 (MQ=255)
                                                                                                                                   cgccgcTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATACACAGACcagc  <  5:155893/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                          
TTATCTTAAACAACTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATACACAGACCAGC  >  NC_002947/4944603‑4944881

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: