Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,609,000 G→C 15.2% P200A (CCT→GCT)  waaA ← 3‑deoxy‑D‑manno‑octulosonic acid transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,609,0000GC15.2% 168.9 / 21.2 66P200A (CCT→GCT) waaA3‑deoxy‑D‑manno‑octulosonic acid transferase
Reads supporting (aligned to +/- strand):  ref base G (26/30);  new base C (5/5);  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 = 8.54e-01

CACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTGGCCAG  >  NC_002947/5608859‑5609144
                                                                                                                                              |                                                                                                                                                
cACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGc                                                                                                                                            <  8:181780/148‑1 (MQ=255)
   cgtcgtGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGATCTTCATCGAGCCGGTGACTTTCACGCGCTCAGGGCGTGCAcc                                                                                                                                         >  6:321342/1‑148 (MQ=255)
        tGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTc                                                                                                                                   <  4:154455/149‑1 (MQ=255)
              atgtcgtcAATCCACACATGGCGCTGGCCGCCCCCCAACTGTTCAAGCAATGCCCTCTCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑TGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGaa                                                                                                                             <  8:228559/142‑1 (MQ=255)
               tGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAg                                                                                                                            >  2:19995/1‑149 (MQ=255)
                   ggcAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGAATTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGaa                                                                                                                             >  7:228559/1‑144 (MQ=255)
                     cAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTgg                                                                                                                      <  6:11974/149‑1 (MQ=255)
                         ccACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGGACGCGCTCAGGGCGTGCGCCCAGGTCACGGACGCGCGGGGCCt                                                                                                                  >  4:294913/1‑149 (MQ=255)
                           acacaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCCGAACTTGATCGAGCCGGTGACTTGCACGCCCTGAGCGCGTGc                                                                                                                                            <  1:68798/122‑1 (MQ=255)
                           acacaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCCCTGAGCGCGTGc                                                                                                                                            >  2:68798/1‑121 (MQ=255)
                            cacaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCCCTGAGCGCGTGCaa                                                                                                                                          >  5:63267/1‑121 (MQ=255)
                            cacaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCCCTGAGCGCGTGCaa                                                                                                                                          <  6:63267/122‑2 (MQ=255)
                                 ggCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCg                                                                                                          <  5:22521/149‑1 (MQ=255)
                                 ggCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCg                                                                                                          <  6:224874/149‑1 (MQ=255)
                                 ggCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCg                                                                                                          <  4:329958/149‑1 (MQ=255)
                                 ggCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCg                                                                                                          <  1:279896/149‑1 (MQ=255)
                                   cgcTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGt                                                                                                        >  6:27957/1‑149 (MQ=255)
                                   cgcTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGt                                                                                                        <  8:34080/149‑1 (MQ=255)
                                    gcTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTc                                                                                                       <  6:332616/149‑1 (MQ=255)
                                     cTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCt                                                                                                      <  4:166696/149‑1 (MQ=255)
                                          cgccCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACg                                                                                                 >  2:362172/1‑149 (MQ=255)
                                                ccAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGAt                                                                                            >  4:181480/1‑148 (MQ=255)
                                                ccAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGAt                                                                                            >  1:75581/1‑148 (MQ=255)
                                                ccAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGAt                                                                                            >  8:149106/1‑148 (MQ=255)
                                                ccAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCCTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGAt                                                                                            >  6:339797/1‑148 (MQ=255)
                                                 cAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCa                                                                                          <  4:49114/149‑1 (MQ=255)
                                                   aCTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAgg                                                                                        <  6:59594/149‑1 (MQ=255)
                                                    cTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGc                                                                                       <  1:362172/149‑1 (MQ=255)
                                                       ttCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTc                                                                                     >  6:27717/1‑148 (MQ=255)
                                                          aCGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGACTCAtct                                                                                 >  6:241382/1‑149 (MQ=255)
                                                             cAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGAGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCGCATCTCgg                                                                              >  8:88145/1‑149 (MQ=255)
                                                               aTGCCCTCGCGCGTGGCAGCAGCTGGTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGcc                                                                            <  5:321342/149‑1 (MQ=255)
                                                                     tCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATc                                                                      >  8:174138/1‑149 (MQ=255)
                                                                       gcgcgTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCg                                                                     <  6:39742/148‑1 (MQ=255)
                                                                              gcagcaGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGt                                                              <  2:75581/148‑1 (MQ=255)
                                                                                           gtcgATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAg                                                 <  5:290371/148‑1 (MQ=255)
                                                                                           gtcgATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAg                                                 >  6:290371/1‑148 (MQ=255)
                                                                                            tcgATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGt                                                                                                        >  5:236583/1‑92 (MQ=255)
                                                                                            tcgATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGt                                                                                                        <  6:236583/92‑1 (MQ=255)
                                                                                              gATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAgcgc                                              <  5:27717/148‑1 (MQ=255)
                                                                                              gATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGcc                                             <  7:174138/149‑1 (MQ=255)
                                                                                                   tCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAg                                        <  3:181480/149‑1 (MQ=255)
                                                                                                       gT‑CGAACTTGATCTAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCa                                           >  1:350849/1‑142 (MQ=255)
                                                                                                       gT‑CGAACTTGATCTAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCa                                           <  2:350849/142‑1 (MQ=255)
                                                                                                       gT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCAc                                    >  2:156525/1‑149 (MQ=255)
                                                                                                             aCTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTc                                                                                                           <  8:170733/73‑1 (MQ=255)
                                                                                                               ttGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGcc                               >  2:50181/1‑147 (MQ=255)
                                                                                                                 gATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGa                             <  1:19995/147‑1 (MQ=255)
                                                                                                                  aTCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGTATCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCa                                                    <  6:185206/123‑1 (MQ=255)
                                                                                                                  aTCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGTATCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCa                                                    >  5:185206/1‑123 (MQ=255)
                                                                                                                           gggACTTGCACGCGCTCAGGGCGTGCACCCAGGGCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGAc                   <  7:313826/145‑1 (MQ=255)
                                                                                                                           gTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACgg                                  <  7:160437/132‑1 (MQ=255)
                                                                                                                           gTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACgg                                  >  8:160437/1‑132 (MQ=255)
                                                                                                                           gTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCAGAGCCCCGGGCCGAGCGCTCCGACaa                 >  8:313826/1‑149 (MQ=255)
                                                                                                                             tgggtGCACGCCCTGAGCGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGcc                               >  1:324922/5‑133 (MQ=255)
                                                                                                                             tgggtGCACGCCCTGAGCGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGcc                               <  2:324922/129‑1 (MQ=255)
                                                                                                                              gggtGCACGCCCTGAGCGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGa                    >  3:210820/4‑143 (MQ=255)
                                                                                                                              gggtGCACGCCCTGAGCGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGa                    <  4:210820/140‑1 (MQ=255)
                                                                                                                               ggtGCACGCCCTGAGCGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTc                       <  2:114941/137‑1 (MQ=255)
                                                                                                                               ggtGCACGCCCTGAGCGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTc                       >  1:114941/3‑139 (MQ=255)
                                                                                                                               cTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAgcgc             <  7:149106/149‑1 (MQ=255)
                                                                                                                                ttGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAgcgcg            >  7:245251/1‑149 (MQ=255)
                                                                                                                                ttGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAgcgcg            >  6:203755/1‑149 (MQ=255)
                                                                                                                                    aCGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCAtt        >  3:37101/1‑149 (MQ=255)
                                                                                                                                         cTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTGGCCa   <  8:245251/149‑1 (MQ=255)
                                                                                                                                          tCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTGGCCAg  >  2:63759/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                
CACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGT‑CGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTGGCCAG  >  NC_002947/5608859‑5609144

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