Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,945,616 C→G 5.8% R322P (CGG→CCG)  flhB ← flagellin export apparatus substrate specificity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,945,6160CG 5.8% 222.5 / 3.3 69R322P (CGG→CCG) flhBflagellin export apparatus substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base C (32/33);  new base G (2/2);  total (34/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 = 9.93e-01

GCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATG  >  NC_002947/4945470‑4945764
                                                                                                                                                   |                                                                                                                                                    
gCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCcgc                                                                                                                                                    >  7:122941/1‑149 (MQ=255)
  aGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCgcgc                                                                                                                                                  <  8:290669/149‑1 (MQ=255)
    gtcgtcCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCa                                                                                                                                                <  3:173649/149‑1 (MQ=255)
        tcCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTCGAGTAATAGATCGCCGGCGCTCgcgc                                                                                                                                            >  8:189580/1‑149 (MQ=255)
           ttCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCgg                                                                                                                                         >  4:113800/1‑149 (MQ=255)
           ttCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCgg                                                                                                                                         >  6:245705/1‑149 (MQ=255)
            tCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGc                                                                                                                                        >  5:180321/1‑149 (MQ=255)
              aGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGa                                                                                                                                      <  8:122941/149‑1 (MQ=255)
                tGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACTTACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTCGAGTAATAGATCGCCGGCGCTTGCGCGGGCGAtc                                                                                                                                    <  7:189580/149‑2 (MQ=255)
                  gCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTc                                                                                                                                   >  4:304108/1‑148 (MQ=255)
                      tGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGa                                                                                                                              >  8:172794/1‑149 (MQ=255)
                        ggCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGATTCGGTGGAGTAATAAATCGCCGGCGCTGGCGCGGGCGAtctattac                                                                                                                              <  1:71234/148‑8 (MQ=255)
                        ggCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATc                                                                                                                            >  6:120620/1‑149 (MQ=255)
                          cgcTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTg                                                                                                                          <  3:209328/149‑1 (MQ=255)
                               gcccttgccGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGAtg                                                                                                                      <  3:252987/148‑1 (MQ=255)
                               gcccttgccGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGAt                                                                                                                       <  3:29089/147‑1 (MQ=255)
                               gcccttgccGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTGCGATGGAGTAATAGATCGCCCGCGCCAGCGCCGACGATTCGAGGAGAGGGAtgt                                                                                                                     >  8:66213/1‑149 (MQ=255)
                                cccttgccGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGTATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGAtgtt                                                                                                                    <  4:118714/149‑1 (MQ=255)
                                cccttgacGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGAtgtt                                                                                                                    <  3:306393/149‑1 (MQ=255)
                                 ccttgccGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTATTAGATCGCCCGCGCCGGCGCCGGCGATGCGGGGAGCTGCAtgttg                                                                                                                   >  7:167166/1‑149 (MQ=255)
                                       cGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTc                                                                                                              >  5:54403/1‑148 (MQ=255)
                                              gTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCaa                                                                                                      >  4:127085/1‑149 (MQ=255)
                                                    gCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCAc                                                                                                <  6:54403/149‑1 (MQ=255)
                                                    gCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCAc                                                                                                <  3:304108/149‑1 (MQ=255)
                                                     cGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACg                                                                                               <  2:193248/149‑1 (MQ=255)
                                                      gAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACg                                                                                               >  5:176661/1‑148 (MQ=255)
                                                            ggAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCtt                                                                                         >  2:196140/1‑148 (MQ=255)
                                                             gAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGATCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCAACGCCAATTTCACGAATCTTc                                                                                        >  5:234507/1‑148 (MQ=255)
                                                                  caTACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGcc                                                                                  <  7:231035/149‑1 (MQ=255)
                                                                   aTACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCa                                                                                 <  5:245705/149‑1 (MQ=255)
                                                                    tACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCAt                                                                                <  5:120620/149‑1 (MQ=255)
                                                                    tACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCAt                                                                                >  3:35886/1‑149 (MQ=255)
                                                                      cGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATg                                                                               >  4:151109/1‑148 (MQ=255)
                                                                       gCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGa                                                                              >  2:121927/1‑148 (MQ=255)
                                                                         cAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAAt                                                                           >  8:299781/1‑149 (MQ=255)
                                                                          aGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATc                                                                          <  7:172794/149‑1 (MQ=255)
                                                                          aGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATc                                                                          <  3:151109/149‑1 (MQ=255)
                                                                                      aCCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGc                                                              <  1:121927/149‑1 (MQ=255)
                                                                                      aCCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGc                                                              >  4:306636/1‑149 (MQ=255)
                                                                                          ccAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCt                                                                       <  7:288275/136‑1 (MQ=255)
                                                                                          ccAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGAGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCt                                                                       >  8:288275/1‑136 (MQ=255)
                                                                                              gTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCagcagc                                                       >  4:63312/1‑148 (MQ=255)
                                                                                                  aGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACgg                                                  <  3:113800/149‑1 (MQ=255)
                                                                                                  aGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACgg                                                  >  3:277467/1‑149 (MQ=255)
                                                                                                   gCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTGCACGCCAATTTCACGAATCTTCAACGCCATGAAATAGCTACCCTTGGCCAGCAGCAACGgc                                                 >  3:277553/1‑149 (MQ=255)
                                                                                                    ccccGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGgcg                                                >  1:254956/1‑149 (MQ=255)
                                                                                                     cgctccgGGATCTCCTGCTCATGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGgcgc                                               <  6:234507/145‑1 (MQ=255)
                                                                                                      ccgccgGGATCTCATGCTCAAGTTCGGTGGAGTAATACATCGCCCGCGCCAGCGCCGGCGAATCGACGATCTGGATGTTGTGCTCCACTACAATTTCACGAATCTTCACAGCCATGAAAAAGCAACCCTTGGCCAGCAGCAACGgcgca                                              >  8:259551/1‑148 (MQ=255)
                                                                                                       cgccgGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGcgccg                                             <  3:63312/149‑1 (MQ=255)
                                                                                                         ccgGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTAc                                                                     <  8:29078/123‑1 (MQ=255)
                                                                                                         ccgGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTAc                                                                     >  7:29078/1‑123 (MQ=255)
                                                                                                          cgGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCAc                                          <  7:299781/149‑1 (MQ=255)
                                                                                                                  ccTGCGCCAGTTCCGTGGAGTAATCGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAACTCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTtc                                  <  4:35886/149‑1 (MQ=255)
                                                                                                                   cTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTtc                                  <  6:122910/148‑1 (MQ=255)
                                                                                                                   cTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGACCAGCAGCAACGGCGCCGCACCGCCtt                                    <  6:176661/146‑1 (MQ=255)
                                                                                                                    tGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGAtgt                                                                                                                     >  7:217682/1‑64 (MQ=255)
                                                                                                                    tGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGAtgt                                                                                                                     <  8:217682/64‑1 (MQ=255)
                                                                                                                      cTCAAGTTCGGTGGAGTAATAGATGGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAAt                                                                           >  4:172423/1‑104 (MQ=255)
                                                                                                                      cTCAAGTTCGGTGGAGTAATAGATGGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAAt                                                                           <  3:172423/104‑1 (MQ=255)
                                                                                                                       tCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCgg                             >  3:137061/1‑149 (MQ=255)
                                                                                                                         aaGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGt                           >  6:58511/1‑149 (MQ=255)
                                                                                                                                ggTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTAtt                     <  2:275323/149‑1 (MQ=255)
                                                                                                                                     aGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCaaa                <  3:127085/149‑1 (MQ=255)
                                                                                                                                         tcgaaTCGCCGGCGCTGGCGCGGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCa            >  6:197569/5‑149 (MQ=255)
                                                                                                                                         aTAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCa            >  1:104053/1‑149 (MQ=255)
                                                                                                                                          tAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCAc           >  3:202907/1‑149 (MQ=255)
                                                                                                                                           aGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCAcc          <  6:180321/149‑1 (MQ=255)
                                                                                                                                            gATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCAcc          <  1:196140/148‑1 (MQ=255)
                                                                                                                                              tCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTAtt                     >  1:33168/1‑135 (MQ=255)
                                                                                                                                              tCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTAtt                     <  2:33168/135‑1 (MQ=255)
                                                                                                                                                 cccGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCTGGGTCGTATTGCAAAGCCAcc          <  7:3027/143‑1 (MQ=255)
                                                                                                                                                 cccGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCTGGGTCGTATTGCAAAGCCAcc          >  8:3027/1‑143 (MQ=255)
                                                                                                                                                  ccGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAAt   <  2:254956/149‑1 (MQ=255)
                                                                                                                                                   cgcgCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATg  <  8:167166/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATG  >  NC_002947/4945470‑4945764

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: