Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,944,740 C→T 6.4% intergenic (+100/+698) 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,7400CT 6.4% 276.8 / 5.2 93intergenic (+100/+698)PP_5673/flhBD‑alanine‑‑D‑alanine ligase/flagellin export apparatus substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base C (40/47);  new base T (3/3);  total (43/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.48e-01

ACTGGTACTTATCTTAAACAACTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATACACAGACCAGC  >  NC_002947/4944595‑4944881
                                                                                                                                                  |                                                                                                                                             
aCTGGTACTTATCTTAAACAACTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATTATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCCCCCAATTGAGATTAAGCGCCGCTCaa                                                                                                                                             >  8:179594/1‑148 (MQ=255)
        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                                                                                                                                  >  3:153897/1‑124 (MQ=255)
                                gCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGTTTGAGCGGCGCttaa                                                                                                                                  <  4:153897/127‑4 (MQ=255)
                                gCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACgg                                                                                                            >  8:221236/1‑149 (MQ=255)
                                gCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACAGCGCTCGATCTTGAGCAATCCCAGAAACgg                                                                                                            <  4:79678/149‑1 (MQ=255)
                                 ccAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGt                                                                                                           <  6:314858/149‑1 (MQ=255)
                                    aTACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGTTTGAGCGGCGCTTAATCtcaa                                                                                                                            >  1:178176/1‑126 (MQ=255)
                                    aTACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGTTTGAGCGGCGCTTAATCtcaa                                                                                                                            <  2:178176/129‑4 (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                                    >  8:288916/1‑130 (MQ=255)
                                                                                                                            ccAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCaaa                                    <  3:251077/130‑1 (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)
                                                                                                                                                  |                                                                                                                                             
ACTGGTACTTATCTTAAACAACTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATACACAGACCAGC  >  NC_002947/4944595‑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: