Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,975,208 C→T 6.7% L132L (CTG→CTA)  flgK ← flagellar hook‑associated protein FlgK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,975,2080CT 6.7% 173.8 / 2.3 60L132L (CTG→CTA) flgKflagellar hook‑associated protein FlgK
Reads supporting (aligned to +/- strand):  ref base C (27/28);  new base T (2/2);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

TCGGCAATGTTGGCTGTGAGCTTGTTCACCTGGCCGGCCATGGTGTCGAGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGAACGTCAC  >  NC_002947/4975071‑4975348
                                                                                                                                         |                                                                                                                                            
tCGGCAATGTTGGCTGTGAGCTTGTTCACCTGGCCGGCCATGGTGTCGAGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGgcgc                                                                                                                                     >  1:61634/1‑147 (MQ=255)
  ggCAATGTTGGCTGTGAGCTTGTTCACCTGGCCGGCCATGGTGTCGAGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGa                                                                                                                                   <  7:45470/147‑1 (MQ=255)
  ggCAATGTTGGCTGTGAGCTTGTTCACCTGGCCGGCCATGGTGTCGAGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGa                                                                                                                                   <  2:320855/147‑1 (MQ=255)
  ggCAATGTTGGCTGTGAGCTTGTTCACCTGGCCGGCCATGGTGTCGAGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAgg                                                                                                                                 >  7:83752/1‑149 (MQ=255)
  ggCAATGTTGGCTGTGAGCTTGTTCACCTGGCCGGCCATGGTGTCGAGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAg                                                                                                                                  <  6:373539/148‑1 (MQ=255)
               gtgAGCTTGTTCACCTGGCCGGCCATGGTGTCGAGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCgg                                                                                                                    >  7:386133/1‑149 (MQ=255)
                     ttGTTCACCTGGCCGGCCATGGTGTCGAGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTcgc                                                                                                              >  7:55783/1‑149 (MQ=255)
                            ccTGGCCGGCCATGGTGTCGAGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCa                                                                                                       <  2:61634/149‑1 (MQ=255)
                             cTGGCCGGCCATGGTGTCGAGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAg                                                                                                      >  1:42850/1‑149 (MQ=255)
                                 ccggccATGGTGTCGAGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGGCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTg                                                                                                  <  8:386133/149‑1 (MQ=255)
                                   ggccATGGTGTCGAGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCtt                                                                                                                 >  1:167771/1‑132 (MQ=255)
                                   ggccATGGTGTCGAGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCtt                                                                                                                 <  2:167771/132‑1 (MQ=255)
                                      cATGGTGTCGAGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGTAGCTGGCGAGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAgcg                                                                                             >  1:119854/1‑149 (MQ=255)
                                      cATGGTGTCGAGCTGGGAGTTGATCGTGGCGTTCTCCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGTAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAgcg                                                                                             >  1:115015/1‑149 (MQ=255)
                                         ggTGTCGAGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGcc                                                                                           >  5:388701/1‑148 (MQ=255)
                                            gtCGAGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCgag                                                                                        >  7:204766/1‑148 (MQ=255)
                                            gtCGAGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCgag                                                                                        >  7:210578/1‑148 (MQ=255)
                                             tCGAGCTTGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGagaa                                                                                      <  8:210578/149‑1 (MQ=255)
                                             tCGAGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGagaa                                                                                      <  8:55783/149‑1 (MQ=255)
                                             tCGAGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGagaa                                                                                      <  8:204766/149‑1 (MQ=255)
                                               gAGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGcc                                                                                           <  6:166322/142‑1 (MQ=255)
                                               gAGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGcc                                                                                           >  5:166322/1‑142 (MQ=255)
                                                aGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAgaagaa                                                                                   <  6:388701/149‑1 (MQ=255)
                                                 gCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGcc                                                                                           <  1:260408/140‑1 (MQ=255)
                                                 gCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGcc                                                                                           >  2:260408/1‑140 (MQ=255)
                                                               gTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAgg                                                                                 <  8:354527/136‑1 (MQ=255)
                                                               gTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAgg                                                                                 >  7:354527/1‑136 (MQ=255)
                                                                tGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGACGAGAAGAAGGCGGTAAGGACCGAg                                                                   >  5:238444/1‑149 (MQ=255)
                                                                tGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCAGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAg                                                                   >  5:244316/1‑149 (MQ=255)
                                                                     tgctgctgGGTCAGCTGGGGGCTCACCGCGGTGAAACGGTTGCTCCGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGAt                                                              <  8:340802/147‑1 (MQ=255)
                                                                             gctCAGCTCGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGGCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGgtgt                                                      <  4:4909/147‑1 (MQ=255)
                                                                               tCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCa                                                    >  5:135499/1‑149 (MQ=255)
                                                                                  gatGGGTGCTCACCGCGTTGAAACGGTTTCTCAGGGTTTGCGCCTGGGTCAGCAGAAGCTGGCGCGAGGAAACGTCACCCGGCTTCGCCCCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACAGAGCTGATGACGGTGTCACg                                                  >  1:344280/3‑148 (MQ=255)
                                                                                  gCTGTGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGTGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGTACCGAGCTGATGCCGGTGTCACg                                                  <  2:344280/148‑1 (MQ=255)
                                                                                  gCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACgg                                                 >  7:315442/1‑149 (MQ=255)
                                                                                      ggTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGCTAAGAACCGAACTGATGCCGGTGTCACGGTCgg                                             >  8:367360/1‑149 (MQ=255)
                                                                                          ctctCCGCGTTTAAACGGTTGCTCAGGGTTTGCGCATGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCTCCCGGCTTCGGCGCCGCAGTCTGCAGCGCCGAGTAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCAGGGTCGGCca                                          <  7:367360/148‑1 (MQ=255)
                                                                                          cTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTcgccg                                                                                                            >  1:339367/1‑82 (MQ=255)
                                                                                          cTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTcgccg                                                                                                            <  2:339367/82‑1 (MQ=255)
                                                                                          cTCACCGCGTGGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCca                                          <  8:83752/148‑1 (MQ=255)
                                                                                             aCCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGTAGCTGGGGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGGCGTGATGAAGGCGGTAAGGACCGAGCTGATGGCGGTGTCAGGGTCGGCcagc                                        <  2:115015/147‑1 (MQ=255)
                                                                                             aCCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGTAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCcagc                                        <  2:119854/147‑1 (MQ=255)
                                                                                             aCCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCcagcag                                      >  6:56714/1‑149 (MQ=255)
                                                                                               cgcgTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGgtgt                                                      <  7:197772/131‑1 (MQ=255)
                                                                                               cgcgTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGgtgt                                                      >  8:197772/1‑131 (MQ=255)
                                                                                                 cgTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGc                                            <  6:374479/139‑1 (MQ=255)
                                                                                                 cgTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGc                                            >  5:374479/1‑139 (MQ=255)
                                                                                                                tCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGTAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGtgatc                   >  3:269576/1‑149 (MQ=255)
                                                                                                                  aGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACg                                                  <  3:407251/116‑1 (MQ=255)
                                                                                                                  aGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACg                                                  >  4:407251/1‑116 (MQ=255)
                                                                                                                  aGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGtgatc                   <  2:42850/147‑1 (MQ=255)
                                                                                                                   gggTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGtgatctga                <  8:315442/149‑1 (MQ=255)
                                                                                                                    ggTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGtgatctgat               <  1:29081/149‑1 (MQ=255)
                                                                                                                      tttGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTgg                        <  7:234971/138‑1 (MQ=255)
                                                                                                                      tttGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTgg                        >  8:234971/1‑138 (MQ=255)
                                                                                                                           gcCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGaa        >  6:279082/1‑149 (MQ=255)
                                                                                                                               gggTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGAACGt     >  2:364544/1‑148 (MQ=255)
                                                                                                                                 gTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTTATGCCGTTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGAACGTCAc  <  6:135499/149‑1 (MQ=255)
                                                                                                                                  tCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGa                                                               <  2:59/87‑1 (MQ=255)
                                                                                                                                  tCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGa                                                               >  1:59/1‑87 (MQ=255)
                                                                                                                                         |                                                                                                                                            
TCGGCAATGTTGGCTGTGAGCTTGTTCACCTGGCCGGCCATGGTGTCGAGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGAACGTCAC  >  NC_002947/4975071‑4975348

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