Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,936,919 A→T 9.2% 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.2% 180.0 / 8.0 65intergenic (‑330/‑88)flhA/PP_4345flagellar export pore protein/GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (28/31);  new base T (4/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.05e-01

CAATGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATAAAACCG  >  NC_002947/4936779‑4937061
                                                                                                                                            |                                                                                                                                                   
cAATGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGcc                                                                                                                                              >  6:29543/1‑148 (MQ=255)
 aaTGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCGATGACTACTCAAAGGCATGAACATCATGCCa                                                                                                                                             >  4:19937/1‑148 (MQ=255)
  aTGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGATGTTCATGccttt                                                                                                                                           <  2:294432/149‑4 (MQ=255)
  aTGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAAc                                                                                                                                           <  6:90851/149‑1 (MQ=255)
  aTGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAAc                                                                                                                                           >  4:307406/1‑149 (MQ=255)
      ccAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCt                                                                                                                                       <  5:279994/149‑1 (MQ=255)
       gatatgCTTACGCTGTTGGATATCCATAGAGTGGGTGGCTCCAGGACAAAAAAAATCTGAGCGACTTTTTCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTc                                                                                                                                      <  8:313920/144‑1 (MQ=255)
       cAAAGTCTTAAGGTCTTGGATAGCCACAGAGAGGGTGTCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTGAGGGTTTCTATGCGGATGTAGAGAGGCAATGACTACTCAAAGGCATGAACATCATGc                                                                                                                                               <  8:42505/140‑1 (MQ=255)
       cAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTc                                                                                                                                      >  7:163587/1‑149 (MQ=255)
       cAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGCGTGGCGCCACGACAAAAAAAATCTGAGCGACTTCATCATGCAAAAGAACACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAACGCATGAACAACATGc                                                                                                                                               >  7:42505/1‑140 (MQ=255)
          aGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGAt                                                                                                                                   <  8:35275/149‑1 (MQ=255)
           ggCTTAAGCTCTTGGAAAGCCACAGATAGGGTGGCGCCACGACAAAAATAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATg                                                                                                                                  <  3:34604/149‑1 (MQ=255)
                aaGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTa                                                                                                                               <  3:307406/147‑1 (MQ=255)
                aaGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGc                                                                                                                             >  7:361069/1‑149 (MQ=255)
                  gCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGctc                                                                                                                           >  7:229209/1‑149 (MQ=255)
                      ttGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTa                                                                                                                         <  2:270457/147‑1 (MQ=255)
                        ggAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGCTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑Gt                                                                                                                     <  1:119054/147‑1 (MQ=255)
                              ccACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCtgt                                                                                                           >  5:179888/1‑149 (MQ=255)
                               cacaGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCGATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCtgtg                                                                                                          <  3:19937/149‑1 (MQ=255)
                                  agagagGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATa                                                                                                       >  1:275843/1‑149 (MQ=255)
                                              gcCACGACAAAAAAAATCTGAGCTACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGa                                                                                                         <  1:241390/135‑1 (MQ=255)
                                              gcCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGa                                                                                                         >  2:241390/1‑135 (MQ=255)
                                                       aaaaaaaTCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACt                                                                                 <  3:36609/149‑1 (MQ=255)
                                                        aaaaaaTCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTa                                                                                >  1:105557/1‑149 (MQ=255)
                                                                tGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGc                                                                          >  1:83380/1‑147 (MQ=255)
                                                                    cGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGa                                                                                   <  5:64066/134‑1 (MQ=255)
                                                                    cGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGa                                                                                   >  6:64066/1‑134 (MQ=255)
                                                                       cTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCtt                                                                                       <  2:325634/127‑1 (MQ=255)
                                                                       cTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCtt                                                                                       >  1:325634/1‑127 (MQ=255)
                                                                       cTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCg                                                                  >  5:384845/1‑148 (MQ=255)
                                                                          gATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGgag                                                               >  1:317812/1‑148 (MQ=255)
                                                                               tGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACAtg                                                         <  8:163587/149‑1 (MQ=255)
                                                                               tGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACAtg                                                         >  6:104985/1‑149 (MQ=255)
                                                                                 cAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACAtggt                                                       >  4:106079/1‑149 (MQ=255)
                                                                                 cAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACAtgg                                                        <  2:83380/148‑1 (MQ=255)
                                                                                  aaaaGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGTCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACAtggtg                                                      <  4:118481/149‑1 (MQ=255)
                                                                                  aaaaGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACAtggtg                                                      <  5:29543/149‑1 (MQ=255)
                                                                                  aaaaGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACAtggtg                                                      >  3:118481/1‑149 (MQ=255)
                                                                                   aaaGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACAtggtgg                                                     <  1:251377/149‑1 (MQ=255)
                                                                                     aGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCg                                                   <  2:105557/149‑1 (MQ=255)
                                                                                        tCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAAc                                                >  1:21176/1‑149 (MQ=255)
                                                                                            ggTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAActctc                                            <  5:104985/149‑1 (MQ=255)
                                                                                              tttCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGATCCGGAGAACATGGTGGCGAACtctctc                                          <  8:361069/149‑1 (MQ=255)
                                                                                              tttCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACtctctc                                          <  2:275843/149‑1 (MQ=255)
                                                                                                  tATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGGGTAGGAGTCTGTGATATGATACCCCAGCGATTTTGACAACATAGCTAGCCCAGTAGAACATGGTGGCGaaatatc                                            >  1:117644/1‑143 (MQ=255)
                                                                                                   aTGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGAccc                                                                   >  7:5994/1‑119 (MQ=255)
                                                                                                   aTGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGAccc                                                                   <  8:5994/119‑1 (MQ=255)
                                                                                                    tGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGAt                                       <  2:21176/146‑1 (MQ=255)
                                                                                                           gTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGgag                                                               >  6:190950/1‑115 (MQ=255)
                                                                                                           gTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGgag                                                               <  5:190950/115‑1 (MQ=255)
                                                                                                              gagagCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTc                          <  8:229209/149‑1 (MQ=255)
                                                                                                               agagCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTcc                         <  6:384845/149‑1 (MQ=255)
                                                                                                                      aTGACTACTCAAGGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAAGATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACaa                   >  1:315790/1‑148 (MQ=255)
                                                                                                                          ctactCAAAGGCATGAACACCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATc              >  3:144005/1‑149 (MQ=255)
                                                                                                                                gttggCATGATGTTCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCAt         >  8:177669/4‑148 (MQ=255)
                                                                                                                                aaaGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCAt         >  1:183014/1‑148 (MQ=255)
                                                                                                                                 ttggCATGATGTTCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGc                                                                          <  2:203007/80‑1 (MQ=255)
                                                                                                                                 ttggCATGATGTTCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGc                                                                          >  1:203007/3‑82 (MQ=255)
                                                                                                                                 ttggCATGATGTTCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATa        >  7:266967/3‑148 (MQ=255)
                                                                                                                                 ttggCATGATGTTCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATa        >  3:151568/3‑148 (MQ=255)
                                                                                                                                   ggCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATaaaa     >  7:64299/1‑149 (MQ=255)
                                                                                                                                   ggCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATaaaa     <  2:183014/149‑1 (MQ=255)
                                                                                                                                      aTGAACACCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATAAAACcg  <  4:144005/149‑1 (MQ=255)
                                                                                                                                          aCATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGa                                        <  7:277266/107‑1 (MQ=255)
                                                                                                                                          aCATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGa                                        >  8:277266/1‑107 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
CAATGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTA‑‑GT‑‑GAGTCTGTGATA‑CCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATAAAACCG  >  NC_002947/4936779‑4937061

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: ‑