Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,936,918 C→G 8.8% 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 8.8% 197.4 / 7.3 66intergenic (‑329/‑89)flhA/PP_4345flagellar export pore protein/GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (27/32);  new base G (3/3);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTACGGCAATGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATAAAACCG  >  NC_002947/4936773‑4937061
                                                                                                                                                  |                                                                                                                                               
ctACGGCAATGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTCGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAAcat                                                                                                                                               >  8:32955/1‑148 (MQ=255)
   cGGCAATGGCCAAAGGCTTAAGCTCTTTGAAAGTCAAAGATAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATg                                                                                                                                           <  7:277109/149‑1 (MQ=255)
   cGGCAATGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAAcatcat                                                                                                                                            <  7:32955/148‑1 (MQ=255)
    ggCAATGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGc                                                                                                                                          <  6:258713/149‑1 (MQ=255)
        aTGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGATGTTCATGccttt                                                                                                                                      <  8:710/149‑4 (MQ=255)
         ttgCCAAACGCATAATCCTTTTGAAATCCACACACATGGTTGCGGCAACAAAAAAAAAATCTTACCCAACTTATCATTCAAAACCTCACGGGTTCAAAGCCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAAc                                                                                                                                      <  1:285937/147‑1 (MQ=255)
         tGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACt                                                                                                                                     <  1:209877/149‑1 (MQ=255)
           gCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGATGTTCATGcctttga                                                                                                                                    <  5:138390/148‑6 (MQ=255)
            ccAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGATGTTCATGcctttga                                                                                                                                    >  2:193356/1‑142 (MQ=255)
            ccAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGATGTTCATGcctttga                                                                                                                                    <  1:193356/147‑6 (MQ=255)
             cAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGc                                                                                                                                          <  3:301077/140‑1 (MQ=255)
             cAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTc                                                                                                                                 <  8:117516/149‑1 (MQ=255)
                aGGCTTAAGCTCTTGGAAACCCACAGACAGGGTCGCGCCACGACAAAAAAAATCTGAGCGACTTGATCAAGCAAAAGCTCACGCTTTCTATCCGGATCTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGAt                                                                                                                              <  8:217398/149‑1 (MQ=255)
                 ggCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATg                                                                                                                             <  3:180241/149‑1 (MQ=255)
                    ttAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGCGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATgg                                                                                                                            >  5:226600/1‑147 (MQ=255)
                    ttAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATgg                                                                                                                            >  3:185486/1‑147 (MQ=255)
                    ttAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTa                                                                                                                          <  4:185486/149‑1 (MQ=255)
                        gCTCTTGGAAAGCCACAGAGAGTGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTGGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGctc                                                                                                                      <  8:154271/149‑1 (MQ=255)
                        gCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGctc                                                                                                                      <  4:213649/149‑1 (MQ=255)
                            ttGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCGAATGACTACTCAAAGGCATGAACATCATGCCAAGTGCTCGATGGTAGCTCTAGt                                                                                                                  >  5:192028/1‑149 (MQ=255)
                               gAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAg                                                                                                               >  1:73771/1‑149 (MQ=255)
                               gAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAg                                                                                                               >  4:265192/1‑149 (MQ=255)
                                 aaGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTc                                                                                                             >  3:188314/1‑149 (MQ=255)
                                   gCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCtg                                                                                                           <  6:226600/149‑1 (MQ=255)
                                              ggTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCg                                                                                                <  2:73771/149‑1 (MQ=255)
                                                         gACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGt                                                                                                                  <  3:276754/120‑1 (MQ=255)
                                                         gACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGt                                                                                                                  >  4:276754/1‑120 (MQ=255)
                                                         gACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCtttt                                                                                     <  3:265192/149‑1 (MQ=255)
                                                         gACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCtttt                                                                                     <  5:142596/149‑1 (MQ=255)
                                                           aaaaaaaaaTCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTATAGATCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTTAGTCTGTGATACCTTCGCCAGCGCTTTTGa                                                                                   <  5:374423/148‑1 (MQ=255)
                                                                        aGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGa                                                                      <  8:289993/149‑1 (MQ=255)
                                                                            aCTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCg                                                                  <  8:275525/149‑1 (MQ=255)
                                                                              ttGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGga                                                                <  5:137536/149‑1 (MQ=255)
                                                                                gATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGgag                                                               >  8:308968/1‑148 (MQ=255)
                                                                                 aTCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCtt                                                                                       <  8:75334/123‑1 (MQ=255)
                                                                                 aTCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCtt                                                                                       >  7:75334/1‑123 (MQ=255)
                                                                                    aTGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACAt                                                          >  3:190592/1‑149 (MQ=255)
                                                                                     tGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACAtg                                                         <  5:138548/149‑1 (MQ=255)
                                                                                     tGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACAtg                                                         <  4:188314/149‑1 (MQ=255)
                                                                                      gCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCgg                                                                 <  3:221626/140‑1 (MQ=255)
                                                                                      gCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCgg                                                                 >  4:221626/1‑140 (MQ=255)
                                                                                      gCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACAtgg                                                        >  6:279834/1‑149 (MQ=255)
                                                                                        aaaaGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACAtggtg                                                      >  5:224191/1‑149 (MQ=255)
                                                                                          acgCTCACGTCTTCTAGGCGGATGTAGAGACCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGGGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGc                                                    <  4:261555/147‑1 (MQ=255)
                                                                                                aCGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAg                                                                                                                   >  8:187728/1‑80 (MQ=255)
                                                                                                aCGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAg                                                                                                                   <  7:187728/80‑1 (MQ=255)
                                                                                                aCGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGa                                                                                   >  1:210692/1‑112 (MQ=255)
                                                                                                aCGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGa                                                                                   <  2:210692/112‑1 (MQ=255)
                                                                                                aCGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAActc                                              >  5:334534/1‑149 (MQ=255)
                                                                                                aCGGATTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAAct                                               <  8:309070/148‑1 (MQ=255)
                                                                                                  ggTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAActctc                                            <  7:308968/149‑1 (MQ=255)
                                                                                                   gTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAActctct                                           <  6:224191/149‑1 (MQ=255)
                                                                                                        tATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGa                                        >  5:265818/1‑147 (MQ=255)
                                                                                                           gCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCt                                   >  4:332278/1‑149 (MQ=255)
                                                                                                              ggATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCtt                                  >  5:372341/1‑148 (MQ=255)
                                                                                                               gATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTa                                >  6:86089/1‑149 (MQ=255)
                                                                                                                    agagagCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGtt                           >  2:331416/1‑149 (MQ=255)
                                                                                                                         gCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGa                      >  3:87956/1‑149 (MQ=255)
                                                                                                                          ccAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGAc                     <  5:279834/149‑1 (MQ=255)
                                                                                                                           cAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATAc                                                                                                     <  8:133744/68‑1 (MQ=255)
                                                                                                                           cAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATAc                                                                                                     >  7:133744/1‑68 (MQ=255)
                                                                                                                                     tCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCa          >  6:343178/1‑149 (MQ=255)
                                                                                                                                     tCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCa          >  5:381633/1‑149 (MQ=255)
                                                                                                                                     tCAAAGGCATGAAAATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTATCTCGATAGCTTTATCGTTCCTGACAAACATCGCCa          >  4:344887/1‑149 (MQ=255)
                                                                                                                                       aaaGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCAt         >  2:374590/1‑148 (MQ=255)
                                                                                                                                        ttggCATGATGTTCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATa        >  7:6085/3‑148 (MQ=255)
                                                                                                                                          ggCATGATGTTCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATaaaa     >  3:340639/1‑149 (MQ=255)
                                                                                                                                             aTGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCAATTGACTACATCGCTAGACACGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTACTGACAAACATAGCCATAAAAcag  >  1:311414/1‑147 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
CTACGGCAATGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATAAAACCG  >  NC_002947/4936773‑4937061

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