Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,936,918 C→G 5.1% intergenic (‑329/‑89) flhA ← / → PP_4345 flagellar export pore protein/GntR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,936,9180CG 5.1% 257.4 / 3.7 78intergenic (‑329/‑89)flhA/PP_4345flagellar export pore protein/GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (38/36);  new base G (2/2);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.21e-01

CTCTACGGCAATGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATAAAACCGCGATG  >  NC_002947/4936771‑4937066
                                                                                                                                                    |                                                                                                                                                    
ctctACGGCAATGGCATAAGGCTTAAGCTCTTGTAAAGCCACTCTGTGGGTGGCGCAACGACTAAATAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAAca                                                                                                                                                     <  7:157078/149‑1 (MQ=255)
    aCGGCAATGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAAcatcat                                                                                                                                                 >  7:96082/1‑149 (MQ=255)
     cGGCAATGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAAcatcat                                                                                                                                                 <  8:96082/148‑1 (MQ=255)
     cGGCAATGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAAcatca                                                                                                                                                  <  4:220850/147‑1 (MQ=255)
        cAATGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGCCTACTCAAAGGCCTGAACATCATGca                                                                                                                                              >  8:400392/1‑147 (MQ=255)
        cAATGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGcc                                                                                                                                              >  5:148849/1‑148 (MQ=255)
          aTGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGATGTTCATGccttt                                                                                                                                           <  3:364243/149‑4 (MQ=255)
             gCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGc                                                                                                                                        <  4:347726/149‑1 (MQ=255)
               cAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTc                                                                                                                                      <  8:53055/149‑1 (MQ=255)
                aaaGGCTTAAGCTCTTGGAAAGCCACAGAGATGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTc                                                                                                                                      >  7:218995/1‑148 (MQ=255)
                aaaGGCTTAAGCTCTTGGAAAGCCACAGAGATGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTc                                                                                                                                      >  7:212901/1‑148 (MQ=255)
                aaaGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTc                                                                                                                                      <  6:209412/148‑1 (MQ=255)
                aaaGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTc                                                                                                                                      >  5:209412/1‑148 (MQ=255)
                  aGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGAt                                                                                                                                   <  8:61987/149‑1 (MQ=255)
                    gCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATgg                                                                                                                                 >  7:401853/1‑149 (MQ=255)
                     cTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATgg                                                                                                                                 <  1:125505/149‑1 (MQ=255)
                      ttAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATgg                                                                                                                                 >  7:251851/1‑147 (MQ=255)
                          gCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGctc                                                                                                                           >  4:244142/1‑149 (MQ=255)
                          gCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAActc                                                                                                                           <  4:254341/149‑1 (MQ=255)
                             ctTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAg                                                                                                                        >  7:102736/1‑149 (MQ=255)
                              ttGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGt                                                                                                                       >  7:346657/1‑149 (MQ=255)
                                   aaGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCTACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTc                                                                                                                  <  7:187872/149‑1 (MQ=255)
                                     gCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCtg                                                                                                                <  8:401853/149‑1 (MQ=255)
                                      ccACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCtgt                                                                                                               <  8:5199/149‑1 (MQ=255)
                                       cacaGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCtgtg                                                                                                              <  8:102736/149‑1 (MQ=255)
                                       cacaGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTAACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCtgt                                                                                                               >  1:304834/1‑148 (MQ=255)
                                           gagagGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATAc                                                                                                          <  2:242296/149‑1 (MQ=255)
                                               gggTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTc                                                                                                      <  2:6246/149‑1 (MQ=255)
                                                 gTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGc                                                                                                    <  3:244142/149‑1 (MQ=255)
                                                      gcCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAgcg                                                                                               >  1:267229/1‑149 (MQ=255)
                                                        cACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTa                                                                                                                               <  1:310875/115‑1 (MQ=255)
                                                        cACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTa                                                                                                                               >  2:310875/1‑115 (MQ=255)
                                                        cACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCt                                                                                             <  6:148849/149‑1 (MQ=255)
                                                         aCGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCtt                                                                                            <  1:106278/149‑1 (MQ=255)
                                                          cGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGt                                                                                                                       >  4:38202/1‑121 (MQ=255)
                                                          cGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGt                                                                                                                       <  3:38202/121‑1 (MQ=255)
                                                          cGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCttt                                                                                           <  4:363122/149‑1 (MQ=255)
                                                           gACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCtttt                                                                                          >  5:238575/1‑149 (MQ=255)
                                                             cAAAAAAAATCTGAGCGAAGTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGCGAGCCAATGACTACTCAAAGGCATTAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCtttt                                                                                          <  7:400392/147‑1 (MQ=255)
                                                                aaaaaaTCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTa                                                                                     >  1:181421/1‑149 (MQ=255)
                                                                aaaaaaTCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTa                                                                                     <  6:120539/149‑1 (MQ=255)
                                                                 aaaaaTCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTa                                                                                     >  6:111707/1‑148 (MQ=255)
                                                                   aaaTCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACAt                                                                                  >  4:281202/1‑149 (MQ=255)
                                                                   aaaTCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACAt                                                                                  <  8:346657/149‑1 (MQ=255)
                                                                    aaTCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATc                                                                                 <  1:38659/149‑1 (MQ=255)
                                                                        agAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTa                                                                             >  6:36044/2‑149 (MQ=255)
                                                                          aGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTa                                                                             <  5:111707/147‑1 (MQ=255)
                                                                             ggcTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGAccc                                                                        >  2:114614/3‑149 (MQ=255)
                                                                               cTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCgg                                                                      <  5:199470/149‑1 (MQ=255)
                                                                                ttGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCg                                                                       >  5:404100/1‑147 (MQ=255)
                                                                                  gATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGga                                                                     <  8:149185/147‑1 (MQ=255)
                                                                                  gATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGga                                                                     <  8:251851/147‑1 (MQ=255)
                                                                                        gCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACAtgg                                                             >  7:29984/1‑149 (MQ=255)
                                                                                        gCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACAtgg                                                             >  3:369797/1‑149 (MQ=255)
                                                                                          aaaaGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACAtggtg                                                           <  5:263264/149‑1 (MQ=255)
                                                                                                tCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAAc                                                     >  5:202427/1‑149 (MQ=255)
                                                                                                     gTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAActctct                                                <  8:29936/149‑1 (MQ=255)
                                                                                                         cTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGa                                             <  2:181421/148‑1 (MQ=255)
                                                                                                          tATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGa                                             >  5:115365/1‑147 (MQ=255)
                                                                                                               cGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGTTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCtt                                       <  4:313141/149‑1 (MQ=255)
                                                                                                                   tGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATc                                   <  2:304834/149‑1 (MQ=255)
                                                                                                                        agagCCAATGACTACTCAAAGGCATGAACATCATGACAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTcc                              >  7:295384/1‑149 (MQ=255)
                                                                                                                             cAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGAc                          >  5:41103/1‑148 (MQ=255)
                                                                                                                             cAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGAc                          >  2:370824/1‑148 (MQ=255)
                                                                                                                              aaTGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACa                         >  5:294530/1‑148 (MQ=255)
                                                                                                                                      ctCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCg                  <  3:281202/147‑1 (MQ=255)
                                                                                                                                       tCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCtgt                                                                                                               <  4:339690/53‑1 (MQ=255)
                                                                                                                                       tCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCtgt                                                                                                               >  3:339690/1‑53 (MQ=255)
                                                                                                                                         aaaGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCAt              >  7:50126/1‑148 (MQ=255)
                                                                                                                                         aaaGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCAt              >  4:369314/1‑148 (MQ=255)
                                                                                                                                         aaaGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATCGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCAt              >  5:3328/1‑148 (MQ=255)
                                                                                                                                         aaaGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATAGCTAGACAAGGAGAACATGGTGGCGAACTCTCTCTATAGCTTTATAGTTCCAGACAAACATCGCAAt              >  5:27153/1‑148 (MQ=255)
                                                                                                                                          ttggCATGATGTTCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATa             >  4:11865/3‑148 (MQ=255)
                                                                                                                                           tggCATGATGTTCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGaa                                                      <  1:64153/105‑1 (MQ=255)
                                                                                                                                           tggCATGATGTTCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGaa                                                      >  2:64153/2‑106 (MQ=255)
                                                                                                                                           aGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATaaa           <  6:202427/149‑1 (MQ=255)
                                                                                                                                             gCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATAAAAc         >  5:188707/1‑149 (MQ=255)
                                                                                                                                             gCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATAAAAc         >  2:301731/1‑149 (MQ=255)
                                                                                                                                                tGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATAAAACcgc      <  6:404100/149‑1 (MQ=255)
                                                                                                                                                 gAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATAAAACcgc      >  1:25248/1‑148 (MQ=255)
                                                                                                                                                 gAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCACTAGACCCGGAGAACATGTTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCTCCATAAAACcgc      >  4:92550/1‑148 (MQ=255)
                                                                                                                                                 catcatcatGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGACAACATGCTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATAAAACcgcg     <  2:106731/146‑1 (MQ=255)
                                                                                                                                                  aaCATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATAAAACcgcg     <  6:238575/148‑1 (MQ=255)
                                                                                                                                                    catcatGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATAAAACCGCGATg  >  4:198933/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTCTACGGCAATGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATAAAACCGCGATG  >  NC_002947/4936771‑4937066

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: