Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,120,818 A→G 27.5% L246P (CTC→CCC)  wcaJ ← colanic biosynthesis UDP‑glucose lipid carrier transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,120,8180AG27.5% 109.4 / 13.1 63L246P (CTC→CCC) wcaJcolanic biosynthesis UDP‑glucose lipid carrier transferase
Reads supporting (aligned to +/- strand):  ref base A (10/34);  new base G (15/3);  total (26/37)
Fisher's exact test for biased strand distribution p-value = 2.41e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACACTTTGATCGGCTTGCCATCCATGCCGTAGCGAGTCTGGCGGAAAATAACCGGCCCTGGTGAACTGAGTTTCACCGCCAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAGCTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCA  >  NC_000913/2120581‑2121053
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             
cgcatttggacgggCTTGCCATCCCTCCCGTAGGGGTGCGAGCGGAAAAACACCGGCCCTGGTGAACTGAGTTTCACCGCAAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGgag                                                                                                                                                                                                                                              <  1:510250/227‑1 (MQ=255)
 catcttgggacgtcTTGCCAGTCATGCCGTAGCGGTTCTGGCGGAAAATACCCGGCCCTGACGAACTGGGTTTCCCCGCCAGAGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGgaga                                                                                                                                                                                                                                             <  1:589306/226‑1 (MQ=255)
         aTCGGCTTGCCATCCATGCCGTAGCGAGTCTGGCGGAAAATAACCGGCCCTGGTGAACTGAGTTTCACCGCCAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGGAAGCGGCGTGTTGTAAAGCGGCACCACCGGGACGCCGTTCCTCTCTTCGAGGGGTGGATGGGGAATGTTgga                                                                                                                                                                                                                                     >  1:934668/1‑237 (MQ=255)
             gCTTGCCATCCATGCCGTAGCGAGTCTGGCGGAAAATAACCGGCCCTGGTGAACTGAGTTTCACCGCCAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAACGGCGTGTCGTAAAGCGGACCCCCCGGTACGCCGTTCATTTCTTCCAGGCGGGGATGGGGAATGTTGAAGGTa                                                                                                                                                                                                                                 >  1:947276/1‑239 (MQ=255)
               ttGCCATCCATGCCGTAGCGAGTCTGGCGGAAAATAACCGGCCCTGGTGAACTGAGTTTCACCGCCAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAACGGCGTGTCGTAAAGGGGCCCCACCGGTACGCCGTTCATCTCTTCGAGGCGGGAATGGGGAATGTTGAAAGGaaa                                                                                                                                                                                                                               >  1:275251/1‑239 (MQ=255)
                 gCCATCCATGCCGTAGCGAGTCTGGCGGAAAATAACCGGCCCTGGTGAACTGAGTTTCACCGCCAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGTCAGCACAATGTCTTCCGGACGTTTGAGCAGGCGGTTAAACCCGGAAAGCGGGGGAGCGTAAAGCGGCACCCCCGGTACGGCGCTCACCTCCTCGAGGGCGCGAAGGAGAAGGCTGCAGGGaaaaa                                                                                                                                                                                                                             >  2:661488/1‑237 (MQ=255)
                        atGCCGTAGCGAGTCTGGCGGAAAATAACCGGCCCTGGTGAACTGAGTTTCACCGCCAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGGGTCGTAAAGCGGCACCACCGGTACGCCGGTCATCTCTTCGAGGCGTGAATGGGGAATGTTGAAGGGAAAAACGTCggg                                                                                                                                                                                                                      >  1:247344/1‑239 (MQ=255)
                         tGCCGTAGCGGGTCGGGCGGAAAATACCCGGCCTTGGTGAACTGAGTTTCACCGCCAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCgggg                                                                                                                                                                                                                     <  1:1083219/239‑1 (MQ=255)
                         ggcggtggCGAGTCTGGCGGAAAATACCCGGCCCGGGTAAACTGAGTTTCACCGCCAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCgggg                                                                                                                                                                                                                     <  1:863602/232‑1 (MQ=255)
                            cGTAGCGAGTCTGGCGGAAAATAACCGGCCCTGGTGAACTGAGTTTCACCGCCAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAACGGCGTGTCGTAAAGCGGCCCCACCGGTACGCCGTTCATCCCTTCCAGGGGGGAGTGGGGAATGTTTAAGGGAAAGACCTCGGGGGTc                                                                                                                                                                                                                  >  1:407794/1‑239 (MQ=255)
                              tAGCGAGTCTGGCGGAAAATAACCGGCCCTGGTGAACTGAGTTTCACCGCCAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAAACCCGGAAAGCGGCGTGTCGGAAAACGGCCCCACCGGTACGCCGTTCATCTCTTCCAGGCGTGAATGGGGAATTGTTAAGGGAAAGACGGCGGGGGTcag                                                                                                                                                                                                                >  1:1087918/1‑239 (MQ=255)
                              gaggcaGTCTGGGCGAAAATAACCGGCCCTGGTGAACTGAGTTTCACCGCCAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATcag                                                                                                                                                                                                                <  1:399014/234‑1 (MQ=255)
                                  gtgTCGGGGGGAAAATACCCCGCCCGGGGGAACTGGTTTTACCCGCCAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCagcagc                                                                                                                                                                                                            <  2:247344/237‑1 (MQ=255)
                                    ttcaggCGGAAATAAACCGGCCCTGGTGAATGGAGTTTCCCCGCCAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAgcagca                                                                                                                                                                                                           <  1:757851/234‑1 (MQ=255)
                                     tCTGGCGGAAAATAACCGGCCCTGGTGAACTGAGTTTCACCGCCAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGGGTCGGAAAGGTGCACCACCGGTACGCCGTTCATCTCTTCGCGGCGGGAATGGGGAATGTTGAAGGTAAAGACGTCGGGGGTCaacagccac                                                                                                                                                                                                         >  2:855197/1‑236 (MQ=255)
                                     tCTGGCGAAAAATACCCGGCCCTGGTGACCTGAGTTTCACCGCCAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCAc                                                                                                                                                                                                          <  1:165824/238‑1 (MQ=255)
                                        ggCGAAAAATACCCGGCCCTGGTAAACTGGGTTTCACCGCCAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGaa                                                                                                                                                                                                      <  1:715521/239‑1 (MQ=255)
                                                             gTGAACTGAGTTTCACCGCCAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGGATGGGGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGGGGGGTCCGCCCGTTGa                                                                                                                                                                                >  1:573972/1‑239 (MQ=255)
                                                                 accggttttCCCCGCCAGCGCAATCCAGCAACGCCCCGGGGAGATCAGGAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGa                                                                                                                                                                            <  2:1087918/233‑1 (MQ=255)
                                                                  cTGAGTTCCACCGCCAGCGCAATACAGCACAGCCCCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGAc                                                                                                                                                                           <  1:429265/239‑1 (MQ=255)
                                                                                   gccAATACCGACCAACACCGGGGGGGGCCGGAGCAGATAAGGCGCCGCCAGCAAAATGGCATCCTCGCGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGGGTGCCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGGCCAGTTTTTTCACTcgcg                                                                                                                                                          <  2:395733/237‑1 (MQ=255)
                                                                                                agcacCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGGAAAGGGGGGGGCCGCCCGTTGTTGGGCCAGGTTTTTCCCCCGCGCGCCGCCGCCCAt                                                                                                                                             >  1:941945/1‑239 (MQ=255)
                                                                                                  cacCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGGAAGCGGCGTGTCGTAAAGCGGCACCACCGGGACGCCGTTCATCCCTTCGAGGCGTGAATGGCGAATGTTTACGGGAAAAAGGTCGGGGGTCAGAGGGACCGGACAGGGGGGTTTCGCCAGTTGGTTGGCCAGTTTTTTTACTcccggccccggccccactt                                                                                                                                           >  2:68273/1‑224 (MQ=255)
                                                                                                  cacCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGgag                                                                                                                                                                                                                                              >  1:701813/1‑141 (MQ=255)
                                                                                                  cacCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGgag                                                                                                                                                                                                                                              <  2:701813/141‑1 (MQ=255)
                                                                                                  cacCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGGGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAAAAGGGGGGGTCCGCCAGTTGATGGACCAGtttttttccccgccccccgccggccattt                                                                                                                                           >  1:955612/1‑216 (MQ=255)
                                                                                                  cacCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGGGTCCGCCCGTTGATGGACCAGGTTTTTTCCTCGCGCGCGGTGGCACAttt                                                                                                                                           >  1:588820/1‑239 (MQ=255)
                                                                                                  cacCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGCGTCGTAAAGCGGCACCACCGGTAAGCCGGTCATCTCTTCGAGGCGTGAATGGGGAATGTTGAAGGGACAGACGTCGGGGATCAGCACAACCGAAAAGAGGGTGTCCGCCAGTGGGGGGACAAGTTTTTTCACCAGCGCGCCGGCGcacaatg                                                                                                                                           >  2:177755/1‑236 (MQ=255)
                                                                                                   cccGGGGAGATCAGCCGCAAAATAAGCGGCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGGAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTg                                                                                                                                          <  2:559389/238‑1 (MQ=255)
                                                                                                    ccGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTg                                                                                                                                          <  1:855197/238‑1 (MQ=255)
                                                                                                      gggggcgtCAGCAGCAGAAAAAGGCGCGCCAGGCCAAGTTCTTCCGCACGTTTGGGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCAt                                                                                                                                       <  2:1072072/232‑1 (MQ=255)
                                                                                                       gggAGATCAGCAGCAGCATAAGCGTCGCCAGCACAATATCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGGAAGCGGCGTGGCGTAAAGGGGCACCACCGGTACGGCGGTCCACTCTTCGTGGCGGGAAAGGGAGATGTTTACGGGAAAGACCGCCGGGGTCCACAGCCCCCAACCAGGGGGTTCCGCCAGTGGGCGGGCCCGTGTTTTCCCTCCCCCGcccgcggcaaatttg                                                                                                                                         >  2:60296/1‑223 (MQ=255)
                                                                                                       gggAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGGGTCGTAAAGCGGCACCACCGGTACGCCCTTCCTCTCTTCGAGGCGTGAATTGGGAAAGGTGCAGGGAAAGACGTCGGGGGTAAGCGGCCCCGAACAGGGGGGGTCCGCCACTTTGTTGAGAAGGGTTTTTACTCGCccgccgccgaccaagtgcgcc                                                                                                                                      >  2:507305/1‑227 (MQ=255)
                                                                                                        ggtgttccccaggcagatgagggTCGCCCGCACAATGTTTTCCGCCCGTTTGCGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAACGCGGCTCCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATcg                                                                                                                                     <  2:455973/217‑1 (MQ=255)
                                                                                                        ggAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATc                                                                                                                                      <  1:292315/238‑1 (MQ=255)
                                                                                                         gagaTCCGCAGCAAAATAAGCCTCCCCAGCACAATGTCTTCCGCACGTTTGAGCAGGGGGTTAACCCCGGAAAGGGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATcg                                                                                                                                     <  2:955612/238‑1 (MQ=255)
                                                                                                         gagaTCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATcg                                                                                                                                     <  1:316474/238‑1 (MQ=255)
                                                                                                            aTCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCAc                                                                                                                                                                                                          >  1:410012/1‑167 (MQ=255)
                                                                                                            aTCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCAc                                                                                                                                                                                                          <  2:410012/167‑1 (MQ=255)
                                                                                                                 ccacaaaaTAAGCTTCCCCAGCACAATGTTTTCCGCACGTTTGGGCTGGCGGTGAACCCCGGAAAGGGGCGTGTCGTAAAGCGGCCCCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGGCGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGa                                                                                                                            <  2:931371/233‑1 (MQ=255)
                                                                                                                   gcagAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGAGGCGTGTCGTAAAGCGGGCCCACCGGTACGTCCTTCCTCTCTTTGAGGCGTGAATGGGGAATGTTGAAGGGCAAGACGCCGGGGATCCGCCGCGCCCAACAGGGGGGGTTCGCCACCTGAGTGACCAGATTTTTTTCCCTCGCGCCCGCGCACCCTTGCATccccctaaaaaag                                                                                                                          >  2:517542/1‑227 (MQ=255)
                                                                                                                     ggaaacagcgtccccgcgacaATGGCTTCCGCAGGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTGCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGAc                                                                                                                           <  2:588820/218‑1 (MQ=255)
                                                                                                                               tCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGGGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGGTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGGTTTAGACGTTTTGaag                                                                                                                 >  1:931371/1‑234 (MQ=255)
                                                                                                                                 gCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCTCGCCGTCGCACATTTGCATCGCGATATAGCCGTTATGAATCtttc                                                                                                            >  1:198788/1‑236 (MQ=255)
                                                                                                                                              tttCCGCCCGTTTGAGCAGGGCGTTAACCCCGGAAAGGGGCGGTTCGTAAAGCGGCCCCTCCGGTACGCCGTTCATCTGTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTc                                                                                               <  2:275251/237‑1 (MQ=255)
                                                                                                                                                    cccGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCCCCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGAcc                                                                                        <  1:68273/237‑1 (MQ=255)
                                                                                                                                                               cggtggTTAACCCGGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCCCCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAGCTGTTGCa                                                                              <  1:517542/234‑1 (MQ=255)
                                                                                                                                                                 gCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGGACGCCGTTCATCTCTTCGAGGCGTGAATGGGGAATGTTGAAGGTAAAGACGGCGGGGATCAGCAGCACCGAACAGGGGGGGTCCGCCCGTTGATGGGCCAGTTTTTTCACCCGCGCGCCGTCCCCACTTTGCATCGCGCAATAAACGCTTTTAAATTTTCCCCCCTTCcgctcccccccccccctggtgtggg                                                                            >  2:788789/1‑217 (MQ=255)
                                                                                                                                                                         ccctgagagccgcgggtctttAAAGGGGCACGCCCGTGATGCAGTTCACTCCTTCCGGGCGGGACTGGAGAATTGTGAAGGTAAAGGCGTCGGGGATCAGCAGCACCGAACAGGTGTTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACACTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAGCTGTTGCAGGTTACCCGcc                                                                   <  2:410019/219‑1 (MQ=255)
                                                                                                                                                                         ccccGGAAAGCGGGGTTTCGTAAAGCGGCACCACGGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAGCTGTTGCAGGTTACCCGcc                                                                   <  2:407794/239‑1 (MQ=255)
                                                                                                                                                                         ccccGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTTAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAGCTGTTGCAGGTTACCCGcc                                                                   <  1:788789/239‑1 (MQ=255)
                                                                                                                                                                              gAAAGCGGCGTGTCGTAAAGCGTCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGCATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGCACAGGTGGGTTCCGCCAGTTGGTGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCCCTTTCGCGTCCTCCACCAACTGTTCCAAGTTAACCGCCCAGTc                                                              >  2:518093/1‑239 (MQ=255)
                                                                                                                                                                                             ttaaggggcccccgcGGTTCGCCGTTCATCTTTTGGGGGCGTGAATGGGGAATGGTGAAGGTAAAGACGTCGGGGGTCAGCAGCCCCGAACAGGGGGTGTCCGCCAGGTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAGCTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACgccgcc                                               <  2:573972/225‑1 (MQ=255)
                                                                                                                                                                                              aaaGCGGCACCCCCGGTACGCCGTTCATCTCTTCGAGCCGTGAATGGAAAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAGCTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACgccgcc                                               <  1:644095/238‑1 (MQ=255)
                                                                                                                                                                                              aaaGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAGCTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACgccgcc                                               <  1:887180/238‑1 (MQ=255)
                                                                                                                                                                                                           cGGTACTCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAGCTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGt                                 <  1:726797/239‑1 (MQ=255)
                                                                                                                                                                                                                  ccGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCCCGTCCTCGACCACCTGTTGCAGGTTACCCCCCCAGTCGTTCGAAACCCCCCCCGCTTTCGGGCGGTGGTaac                           >  2:254285/1‑237 (MQ=255)
                                                                                                                                                                                                                       tctcttTTCGGGGCGTTATGGGGAAATGTTTAAGGGAAAGGCGTCGGGGATCAGCAGCCCCGAACGGGTGGGGCCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAGCTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCcca                     <  2:934668/237‑1 (MQ=255)
                                                                                                                                                                                                                        cccttttCGCGGCGGGAATGGAGAATGTTGAAGGGAAAGCCGTCGGGGTTCCGGACCCCCGAACAGGGGGGTGCCGCCAGTTGATGGACCCTTTTTTTCCCTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAGCTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCccac                    <  2:941945/234‑1 (MQ=255)
                                                                                                                                                                                                                         actcgTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAGCTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCccacc                   <  1:254285/238‑1 (MQ=255)
                                                                                                                                                                                                                             ttCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCTTCCTCGACCAGCTGTTGCAGGTTACCCGCCCAGCCGTTCGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACtt               >  1:337331/1‑239 (MQ=255)
                                                                                                                                                                                                                               cGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGCTATAGACGTTATGAATCTTGCCCGCTTTCCGGTCCTCGACCCGCGGTTGCAGGTTACCCGCCCAGTCGTTAGAAACCCCGCCCCGTTTCGGGTGGAGGTAAACGCCACCCACTCCa             >  2:449307/1‑239 (MQ=255)
                                                                                                                                                                                                                               cGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGCCCAGCTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCCCCCGGCTTCGGGTGGGGGTAAAccccccccccccca             >  2:356716/1‑233 (MQ=255)
                                                                                                                                                                                                                                gAGGGGTGAATGGAGAATGGTGAAGGTAAAGACGGCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAGCTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCaa            <  1:356716/239‑1 (MQ=255)
                                                                                                                                                                                                                                 ggggCGTGAATGGATAATGTTGAAGGTAAAGGCGGCGGGGATCAGCAGCCCCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTCTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAGCTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCaa            <  1:556639/237‑1 (MQ=255)
                                                                                                                                                                                                                                           tGGGGAATGGTTAAGGTAAAGGCGGCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCTGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAGCTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCa  <  1:456335/239‑1 (MQ=255)
                                                                                                                                                                                                                                           tGGAGATGGTTGAAGGTAAAGCCGTCGGGGATCAGCAGCCCCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAGCTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCa  <  1:147043/239‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             
CACACTTTGATCGGCTTGCCATCCATGCCGTAGCGAGTCTGGCGGAAAATAACCGGCCCTGGTGAACTGAGTTTCACCGCCAGCGCAATACAGCACAGCACCGGGGAGATCAGCAGCAGAATAAGCGTCGCCAGCACAATGTCTTCCGCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAGCTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCA  >  NC_000913/2120581‑2121053

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: