Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,981,916 T→C 6.3% A12A (GCA→GCG flgE ← flagellar hook protein FlgE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,981,9160TC 6.3% 179.7 / 2.3 63A12A (GCA→GCGflgEflagellar hook protein FlgE
Reads supporting (aligned to +/- strand):  ref base T (29/30);  new base C (2/2);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

TTACCCCGGTGCCGATCACGGTCTTGCCCGCGCTGGTGCCGCGAATGGAGTTGGAGTACTGGTCGGCGAACTCGGCACGCGACGACTTGAAGCCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATGCTTTGGATTTTGGAAAGCGCAACACTGCCCAGCCCGGCCAGGTTGAGCGTCATCTCCGAACCGTTGCTACC  >  NC_002947/4981768‑4982063
                                                                                                                                                    |                                                                                                                                                      
ttACCCCGGTGCCGATCACGGTCTTGCCCGCGCTGGTGCCGCGAATGGAGTTGGAGTACTGGTCGGCGAACTCGGCACGCGACGACTTGAAGCCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTtgct                                                                                                                                                        >  3:109905/1‑149 (MQ=255)
           ccGATCACGGTCTTGCCCGCGCTGGTGCCGCGAATGGAGTTGGAGTACTGGTCGGCGAACTCGGCACGCGACGACTTGAAGCCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGcc                                                                                                                                             >  3:2025/1‑149 (MQ=255)
                     tCTTGCCCGCGCTGGTTCCGCGAATGCAGTTGGAGTACTGGGCGGCGAACTCGGCACGCGACGACTTGAAGCCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCg                                                                                                                                   <  2:371473/149‑1 (MQ=255)
                         gCCCGCGCTGGTGCCGCGAATGGAGTTGGAGTACTGGTCGGCGAACTCGGCACGCGACGACTTGAAGCCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGatat                                                                                                                              >  2:370019/1‑149 (MQ=255)
                         gCCCGCGCTGGTGCCGCGAATGGAGTTGGAGTACTGGTCGGCGAACTCGGCACGCGACGACTTGAAGCCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGat                                                                                                                                 <  1:239910/147‑1 (MQ=255)
                         gCCCGCGCTGGTGCCGCGAATGGAGTTGGAGTACTGGTCGGCGAACTCGGCACGCGACGACTTGAAGCCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGa                                                                                                                                  <  3:225854/146‑1 (MQ=255)
                           ccGCGCTGGTGCCGCGAATGGAGTTGGAGTACTGGTCGGCGAACTCGGCACGCGACGACTTGAAGCCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATAtt                                                                                                                             >  6:259174/1‑148 (MQ=255)
                           ccGCGCTGGTGCCGCGAATGGAGTTGGAGTACTGGTCGGCGAACTCGGCACGCGACGACTTGAAGCCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATAtt                                                                                                                             >  4:366900/1‑148 (MQ=255)
                           ccGCGCTGGTGCCGCGAATGGAGTTGGAGTACTGGTCGGCGAACTCGGCACGCGACGACTTGAAGCCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATAtt                                                                                                                             >  7:108397/1‑148 (MQ=255)
                                   gTGCCGCGAATGGAGTTGGAGTACTGGTCGGCGAACTCGGCACGCGACGACTTGAAGCCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACAt                                                                                                                   <  5:81398/149‑1 (MQ=255)
                                      ccGCGAATGGAGTTGGAGTACTGGTCGGCGAACTCGGCACGCGACGACTTGAAGCCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACAtgtg                                                                                                                >  1:140581/1‑149 (MQ=255)
                                      ccGCGAATGGAGTTGGAGTACTGGTCGGCGAACTCGGCACGCGACGACTTGAAGCCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACAtgt                                                                                                                 <  6:45991/148‑1 (MQ=255)
                                        gcgAATGGAGTTGGAGTACTGGTCGGCGAACTCGGCACGCGACGACTTGAAGCCTGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCTATATTGAAAGACATgtgt                                                                                                               >  6:317431/1‑148 (MQ=255)
                                         cgAATGGAGTTGGAGTACTGGTCGGCGAACTCGGCACGCGACGACTTGAAGGCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTaa                                                                                                             <  8:39358/149‑1 (MQ=255)
                                            aTGGAGTTGGAGTACTGGTCGGCGAACTCGGCACGCGACGACTTTAAGCCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTc                                                                                                          <  5:269349/149‑1 (MQ=255)
                                                       gtACTGGTCGGCGAACTCGGCACGCGACGACTTGAAGCCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTg                                                                                               <  5:259174/149‑1 (MQ=255)
                                                           tGGTCGGCGAACTCGGCACGCGACGACTTGAAGCCGGTGGTGGCAACGTTCGCAATGATGTTGCCGGAAACGTTCAGGGCCTGGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGCCATGTGTTACTCCTGTAGTTTTGCt                                                                                             >  7:291006/1‑147 (MQ=255)
                                                             gTCGGCGAACTCGGCACGCGACGACTTGAAGCCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAgg                                                                                         >  4:317617/1‑149 (MQ=255)
                                                             gTCGGCGAACTCGGCACGCGACGACTTGAAGCCGGAGGTGGCAACGTTGGCAATGTTGATGCCTGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAGAGACATGTGTAACACCTGTCGTTTTGCTCAgg                                                                                         >  2:253687/1‑149 (MQ=255)
                                                                       tCGGCACGCGACGACTTGAAGCCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGc                                                                               >  2:280728/1‑149 (MQ=255)
                                                                       tCGGCACGCGACGACTTGAAGCCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGc                                                                               >  2:111759/1‑149 (MQ=255)
                                                                         ggCACGCGACGACTTGAAGCCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCg                                                                             <  4:2025/149‑1 (MQ=255)
                                                                          gCACGCGACGACTTGAAGCCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGc                                                                                                                                              <  3:385547/85‑1 (MQ=255)
                                                                          gCACGCGACGACTTGAAGCCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGc                                                                                                                                              >  4:385547/1‑85 (MQ=255)
                                                                           cACGCGACGACTTGAAGCCGGTTGTGGCAACGTTGGCAATGTTGTTGACGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCAGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGAGCTAGATGCCGAt                                                                           >  8:155561/1‑149 (MQ=255)
                                                                            aCGCGACGACTTGAAGCCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATg                                                                          >  8:115691/1‑149 (MQ=255)
                                                                            aCGCGACGACTTGAAGCCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATg                                                                          <  3:366900/149‑1 (MQ=255)
                                                                            aCGCGACGACTTGAAGCCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCCGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATg                                                                          >  2:316047/1‑149 (MQ=255)
                                                                               cgacgaCTTGAAGCCGGTGGTGGCAGCGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATGCtt                                                                       <  6:219337/149‑1 (MQ=255)
                                                                               cgacgaCCTGAAGCCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATGCtt                                                                       >  5:355656/1‑149 (MQ=255)
                                                                                      ttGAAGCCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGaaa                                                                                                                         <  7:157335/93‑1 (MQ=255)
                                                                                      ttGAAGCCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGaaa                                                                                                                         >  8:157335/1‑93 (MQ=255)
                                                                                        gAAGACGGTGGTGGCAACGTTGGCAATGTTGTTGCCGTTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAATGACATGTGTAACTCTTGTCGTTTTGCTCAGGCCCTAGATGCCGATGCTTTGGAtttt                                                               <  1:75013/148‑1 (MQ=255)
                                                                                         aaGCCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATGCTTTGGATTTTg                                                              <  8:108397/148‑1 (MQ=255)
                                                                                         aaGCCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATGCTTTGGATTTTg                                                              <  6:355656/148‑1 (MQ=255)
                                                                                          tgCCGGTGGTGGCAACGTTGGCAATGTTGTTGCGGGTAACGTTCTGGGCCTTGTTTGATGCATTTAGGCGGCTAAGGCCGATATTTTTTAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATGCTTTGGATTTTg                                                              <  1:111759/147‑1 (MQ=255)
                                                                                          agtcGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATGCTTTGGATTTTGGa                                                            <  5:144005/146‑1 (MQ=255)
                                                                                               gtggtgGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATGCTTTGGATTTTGGAAAgcg                                                       >  8:311678/1‑149 (MQ=255)
                                                                                                 ggtggCAACCTTGGCAATGATGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTCTAACTACTGTCGTTTTGCTCAGGCACTAGATGCCTATGCTTTGGATTTTGGAAAgcgc                                                      >  5:298016/1‑148 (MQ=255)
                                                                                                  gtggCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATGCTTTGGATTTTGGAAAGCGCaa                                                    <  1:370019/149‑1 (MQ=255)
                                                                                                        aCGTTTGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTTTTTGCTGGATAGACGCCGCTAATGCCGATATTGAAAGACATTTGTAACTCCTGTCGTTTTTCTCATGCCCTAGATGCCGATGCTTTGGATTTTGGAAAGCGCAACACTGc                                             <  7:155561/149‑1 (MQ=255)
                                                                                                          gTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATGCTTTGGATTTTGGAAAGCGCAACACTGcc                                            >  2:228116/1‑148 (MQ=255)
                                                                                                             ggCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATGCTTTGGATTTTGGAAAGCGCAACACTGCCCAg                                         <  2:140581/148‑1 (MQ=255)
                                                                                                             ggCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATGCTTTGGATTTTGGAAAGCGCAACACTGCCCAGc                                        >  5:49729/1‑149 (MQ=255)
                                                                                                                aaTGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCCGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATGCTTTGGATTTTGGAAAGCGCAACACTGCCCAGCCCg                                     <  1:316047/149‑1 (MQ=255)
                                                                                                                       ttgCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCa                                                                                           >  8:251262/1‑89 (MQ=255)
                                                                                                                       ttgCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCa                                                                                           <  7:251262/89‑1 (MQ=255)
                                                                                                                       ttgCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATGCTTTGGATTTTGGAAAGCGCAACACTGCCCAGCCCGGCCAGGt                              >  5:66922/1‑149 (MQ=255)
                                                                                                                       ttgCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATGCTTTGGATTTTGGAAAGCGCAACACTGCCCAGCCCGGCCAGGt                              >  2:189409/1‑149 (MQ=255)
                                                                                                                       ttgCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCCCTAAGTCCGAAATTGAAAGACATGTGTAACTCATCTCGTTTTGCTCAGGCCCTACATTCCCATGCTTTGGATTTTGGAAAGCGCAACACTGCCCAGCCCGCCCACGt                              >  7:65364/1‑149 (MQ=255)
                                                                                                                        ggCCGGTAACCGTCAGCGCCTTGTTTGCTGCATAGAGGCCGCTCAGCCCGATATTGAAAGACATGTGTAACTCCTGTCGTCTTGCTCAGGCCCTAGATGCCGATGCTTTGGATTTTGGAAAGCGCAACACTGCCCAGCCCGGCCAGGtt                             <  4:109905/148‑1 (MQ=255)
                                                                                                                         gCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGAGATGTGTAACTCCTGTAGTTTTGCTCAGGCCCTAGATGCCGATGCTTTGGATTTTGGAAAGCGCAACACTGCCCAGCCCGGCCAGGtt                             <  1:253687/148‑1 (MQ=255)
                                                                                                                          ccGGTATAGTTCAGGGCCTTGTTTGCTGCATAGAGGCCTCTAAGGCCGATATTGAAAGTCATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATTCCGATGCTTTGGATTTTGGAAAGCGCAACACTGCCCAGCCCGGCCAGGTTGa                           <  3:90377/149‑1 (MQ=255)
                                                                                                                          ccGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATGCTTTGGATTTTGGAAAGCGCAACACTGCCCAGCCCGGCCAGGTTGa                           <  4:291398/149‑1 (MQ=255)
                                                                                                                                actttCAGGGCCTTGTTTGCTGCATAGACGCCGCTAAGTCCGATATTGATAGACATGTGTAACTCCTGTCGTGTTGCTCTGGCCCTAGATGCCGATGCTTTGGATTTTGGAAAGCGCAACACTGCCCAGCCCGGCCAGGTTGAGCGTc                      <  6:298016/145‑1 (MQ=255)
                                                                                                                                aCGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATGCTTTGGATTTTGGAAAGCGCAACACTGCCCAGCCCGGCCAGGTTGAGCGTc                      <  3:317617/148‑1 (MQ=255)
                                                                                                                                  gTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGAGCATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATGCTTTGGATTTTGGAAAGCGCAACACTGCCCAGCCCGGCCAGGTTGAGCGTCAt                    >  1:96976/1‑149 (MQ=255)
                                                                                                                                  gTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATGCTTTGGATTTTGGAAAGCGCAACACTGCCCAGCCCGGCCAGGTTGAGCGTCAtc                   >  8:68073/1‑149 (MQ=255)
                                                                                                                                   ttCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATGCTTTGGATTTTGGAAAGCGCAACACTGCCCAGCCCGGCCAGGTTGAGCGTCAtct                  >  2:27800/1‑149 (MQ=255)
                                                                                                                                        ggCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATGCTTTGGATTTTGGAAAGCGCAACACTGCCCAGCCCGGCCAGGTTGAGCGTCATCTCCGaa             >  6:19390/1‑149 (MQ=255)
                                                                                                                                         gCCTTGTTTGCCGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATGCTTTGGATTTTGGAAAGCGCAACACTGCCCAGCCCGGCCAGGTTGAGCGTCATCTCCGAAc            >  4:82699/1‑149 (MQ=255)
                                                                                                                                          ccTTGTTTGCCGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATGCTTTGGATTTTGGAAAGCGCAACACTGCCCAGCCCGGCCAGGTTGAGCGTCATCTCCGAAcc           <  3:82699/149‑1 (MQ=255)
                                                                                                                                              gTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATGCTTTGGATTTTGGAAAGCGCAACACTGCCCAGCCCGGCCAGGTTGAGCGTCATCTCCGAACCGtt        <  7:115691/148‑1 (MQ=255)
                                                                                                                                               tttGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATGCTTTGGATTTTGGAAAGCGCAACACTGCCCAGCCCGGCCAGGTTGAGCGTCATCTCCGAACCGTTg       <  6:49729/148‑1 (MQ=255)
                                                                                                                                                   ctgcATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATGCTTTGGATTTTGGAAAGCGCAACACTGCCCAGCCCGGCCAGGTTGAGCGTCATCTCCGAACCGTTGCTAcc  >  6:133840/1‑149 (MQ=255)
                                                                                                                                                    tgcATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATGCTTTGGATTTTGGAAAGCTAGGTAGCATACCAACCCCGCCAGGTTGATCGTAATATACGAACCGTTgatta   >  1:324638/1‑144 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
TTACCCCGGTGCCGATCACGGTCTTGCCCGCGCTGGTGCCGCGAATGGAGTTGGAGTACTGGTCGGCGAACTCGGCACGCGACGACTTGAAGCCGGTGGTGGCAACGTTGGCAATGTTGTTGCCGGTAACGTTCAGGGCCTTGTTTGCTGCATAGAGGCCGCTAAGGCCGATATTGAAAGACATGTGTAACTCCTGTCGTTTTGCTCAGGCCCTAGATGCCGATGCTTTGGATTTTGGAAAGCGCAACACTGCCCAGCCCGGCCAGGTTGAGCGTCATCTCCGAACCGTTGCTACC  >  NC_002947/4981768‑4982063

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: