Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,936,919 A→T 9.1% intergenic (‑330/‑88) 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,9190AT 9.1% 183.7 / 7.7 66intergenic (‑330/‑88)flhA/PP_4345flagellar export pore protein/GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (28/32);  new base T (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 = 7.82e-01

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: