Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,620,227 C→G 5.6% intergenic (+20/+11) marB → / ← eamA periplasmic mar operon regulator/cysteine and O‑acetyl‑L‑serine efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,620,2270CG 5.6% 372.2 / 10.6 107intergenic (+20/+11)marB/eamAperiplasmic mar operon regulator/cysteine and O‑acetyl‑L‑serine efflux system
Reads supporting (aligned to +/- strand):  ref base C (51/50);  new base G (4/2);  total (55/52)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.45e-01

Only 100 of 111 total aligned reads displayed.

CGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCCAGTCCTACTACGGGCACCAGTAACGATAACGGTGCAACCCGCCAGGTTTCATAGCGTCCCAGTAACGTCCCCCAGATCCCATAACCAACAATTGTCGCCACAAACGCCAGATACATCAGAGACAAG  >  NC_000913/1619973‑1620496
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                             
aggcACTTAGGCAATCATGAAACCACTTTCACCCGCAATAGCAGCTGCGCTTATTCTCGTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCATCCAGTTGTTACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGTGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTccc                                                                                                                                                                                                                                                           <  1:83915/272‑1 (MQ=255)
                        attttCATCCACAAGAGCAGTGTCGCTTGTTCTCCTTTCCGCTCAGGGCGGTGCAGAACAAACCGCTCAGCCAGTGGTTACTTATTGCTCCAATGTCGTGGTTGTTCCCACATCGCAGGAAAGCCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTCTTCGGGCCCCGACAGCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTGACGCCACCGc                                                                                                                                                                                                                                   <  4:91625/273‑1 (MQ=255)
                             cATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACAtc                                                                                                                                                                                                                                                                               >  1:12430/1‑226 (MQ=255)
                                                 cTTATTCTCGTTCCCCGGAAGGGCGTTGCGGAACAAACCCAGAGACCAGTTGTTACTTCTTGTGCCAATGTCGTGGGTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAACTCACATGGATACTGGCAGTGATAAGTCGTATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGcc                                                                                                                                                                                                          <  4:141975/275‑1 (MQ=255)
                                                  ttATTCTCTTTTCCGCGCAGTGCGTTGCGGAACAAACCCCGCAGCCAGTTGTTACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGcc                                                                                                                                                                                                          <  1:57113/274‑1 (MQ=255)
                                                                gcgcAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATCCATTGCTATACAGCCCGGTCATAATGAgca                                                                                                                                                                                           >  2:344522/1‑275 (MQ=255)
                                                                  gcAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCaaa                                                                                                                                                                                                                           >  2:250792/1‑241 (MQ=255)
                                                                   cAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccg                                                                                                                                                                                        <  2:440057/275‑1 (MQ=255)
                                                                      ggCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACa                                                                                                                                                                                                                       <  4:250788/241‑1 (MQ=255)
                                                                                                                 ccAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACAtctc                                                                                                                                                                                                                                                                             <  5:355777/144‑1 (MQ=255)
                                                                                                                 ccAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACAtctc                                                                                                                                                                                                                                                                             >  6:355777/1‑144 (MQ=255)
                                                                                                                   aaTGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCGGGGcca                                                                                                                                                                                                                                                                      <  5:409872/149‑3 (MQ=255)
                                                                                                                   aaTGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCtt                                                                                                                                                                                                                                                                      >  7:198444/1‑149 (MQ=255)
                                                                                                                    aTGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCtta                                                                                                                                                                                                                                                                     <  8:198444/149‑1 (MQ=255)
                                                                                                                            gttgttCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTc                                                                                                                                                                                                                                                             >  5:231587/1‑149 (MQ=255)
                                                                                                                             ttgttCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTcc                                                                                                                                                                                                                                                            >  5:450227/1‑149 (MQ=255)
                                                                                                                             ttgttCCCCCATCGCAGGAACACCCACCATTTGATTTAAATCACTTCGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTcc                                                                                                                                                                                                                                                            >  8:293309/1‑149 (MQ=255)
                                                                                                                               gttCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgc                                                                                                                                                                                            >  2:55577/1‑211 (MQ=255)
                                                                                                                                ttCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCAc                                                                                                                                                                                                                                                         >  8:216055/1‑149 (MQ=255)
                                                                                                                                    cccATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCggg                                                                                                                                                                                                                                                                          >  5:130283/1‑128 (MQ=255)
                                                                                                                                    cccATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCggg                                                                                                                                                                                                                                                                          <  6:130283/128‑1 (MQ=255)
                                                                                                                                    cccATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCttt                                                                                                                                                                                                                                                     >  6:440110/1‑149 (MQ=255)
                                                                                                                                     ccATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTc                                                                                                                                                                                                     <  4:425714/196‑1 (MQ=255)
                                                                                                                                     ccATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGGAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTa                                                                                                                                                                                                                                                    >  8:267803/1‑149 (MQ=255)
                                                                                                                                      cATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAc                                                                                                                                                                                                                                                   <  7:38385/149‑1 (MQ=255)
                                                                                                                                      cATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAc                                                                                                                                                                                                                                                   <  6:318835/149‑1 (MQ=255)
                                                                                                                                        tCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCg                                                                                                                                                                                                                                                 <  6:270106/149‑1 (MQ=255)
                                                                                                                                          gCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGAAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTc                                                                                                                                                                                                     <  1:2570/191‑1 (MQ=255)
                                                                                                                                              gAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGGTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTAc                                                                                                                                                                                                                                           >  6:396187/1‑149 (MQ=255)
                                                                                                                                                   cccaccGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCAc                                                                                                                                                                                                                                      >  7:450748/1‑149 (MQ=255)
                                                                                                                                                   cccaccGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAc                                                                                                                                                                                                              >  1:35721/1‑173 (MQ=255)
                                                                                                                                                    ccaccGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTAc                                                                                                                                                                                                                                           >  7:95900/1‑143 (MQ=255)
                                                                                                                                                        cGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCaa                                                                                                                                                                                                                                 <  7:216055/149‑1 (MQ=255)
                                                                                                                                                        cGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCaa                                                                                                                                                                                                                                 <  7:223317/149‑1 (MQ=255)
                                                                                                                                                          tttGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGc                                                                                                                                                                                                                               >  7:430464/1‑149 (MQ=255)
                                                                                                                                                           ttGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCttt                                                                                                                                                                                                                                             >  8:10079/1‑134 (MQ=255)
                                                                                                                                                              aTTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCaaa                                                                                                                                                                                                                           <  7:340481/149‑1 (MQ=255)
                                                                                                                                                                ttAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAc                                                                                                                                                                                                                                                   >  5:326642/1‑123 (MQ=255)
                                                                                                                                                                ttAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAc                                                                                                                                                                                                                                                   <  6:326642/123‑1 (MQ=255)
                                                                                                                                                                   aaTCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACAt                                                                                                                                                                                                                      <  6:231587/149‑1 (MQ=255)
                                                                                                                                                                    aTCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACAtt                                                                                                                                                                                                                     <  6:383761/149‑1 (MQ=255)
                                                                                                                                                                       acaTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGat                                                                                                                                                                                                                  >  5:215503/1‑149 (MQ=255)
                                                                                                                                                                        caTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCaaa                                                                                                                                                                                                                           >  2:28599/1‑139 (MQ=255)
                                                                                                                                                                         aTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCGACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATTAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCACAAATAGTGACGCACTTACCAGTCCTACTACCGGCACCAGTAACGATAACGGTGCAACCCGCCAgg                                                                                  >  1:70469/1‑275 (MQ=255)
                                                                                                                                                                          tGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGAtata                                                                                                                                                                                                               <  6:450227/149‑1 (MQ=255)
                                                                                                                                                                          tGGGTACTGGCAGTGATAAGGCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGAtata                                                                                                                                                                                                               <  5:396187/149‑1 (MQ=255)
                                                                                                                                                                            ggTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACa                                                                                                                                                                                                                       <  4:28599/139‑1 (MQ=255)
                                                                                                                                                                              tACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGa                                                                                                                                                                                              >  2:467092/1‑162 (MQ=255)
                                                                                                                                                                               aCTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACAt                                                                                                                                                                                                                      >  6:304374/1‑137 (MQ=255)
                                                                                                                                                                               aCTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACAt                                                                                                                                                                                                                      <  5:304374/137‑1 (MQ=255)
                                                                                                                                                                                   gCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCaa                                                                                                                                                                                                                                 >  8:80992/1‑122 (MQ=255)
                                                                                                                                                                                   gCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCaa                                                                                                                                                                                                                                 >  6:422984/1‑122 (MQ=255)
                                                                                                                                                                                   gCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCaa                                                                                                                                                                                                                                 <  7:80992/122‑1 (MQ=255)
                                                                                                                                                                                   gCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCaa                                                                                                                                                                                                                                 <  5:422984/122‑1 (MQ=255)
                                                                                                                                                                                       tGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAc                                                                                                                                                                                                              >  1:56461/1‑137 (MQ=255)
                                                                                                                                                                                         aTAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAAt                                                                                                                                                                                                <  8:430464/149‑1 (MQ=255)
                                                                                                                                                                                           aaGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGcc                                                                                                                                                                                                          <  3:56461/137‑1 (MQ=255)
                                                                                                                                                                                             gTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgc                                                                                                                                                                                            >  8:267526/1‑149 (MQ=255)
                                                                                                                                                                                                 gATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCa                                                                                                                                                                                                                                       >  2:78583/1‑102 (MQ=255)
                                                                                                                                                                                                 gATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCCAGTCCTACTACGGGCACCAGTAACGATAACGGTGCAAcc                                                                                         >  2:59367/1‑244 (MQ=255)
                                                                                                                                                                                                             gTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAAtt                                                                                                                                                                            <  6:215503/149‑1 (MQ=255)
                                                                                                                                                                                                               gCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGc                                                                                                                                                                          <  7:230087/149‑1 (MQ=255)
                                                                                                                                                                                                                cccTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCa                                                                                                                                                                         >  8:453446/1‑149 (MQ=255)
                                                                                                                                                                                                                cccTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCCATTGCAAGTCCTACTACGGGCACCAGTAACGATAACGGTGCAACCCGCCAGGTTTCCTAGCGGCCCAGTAACGGCCCCCAGATCCCATaaa                                           >  4:326857/1‑274 (MQ=255)
                                                                                                                                                                                                                  cTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCGCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGa                                                                                                                                                                       <  7:293309/149‑1 (MQ=255)
                                                                                                                                                                                                                     ttatAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGAccc                                                                                                                                                                    >  6:174192/1‑149 (MQ=255)
                                                                                                                                                                                                                     ttatAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGAccc                                                                                                                                                                    <  7:230643/149‑1 (MQ=255)
                                                                                                                                                                                                                       ataATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGt                                                                                                                                                                  >  7:307972/1‑149 (MQ=255)
                                                                                                                                                                                                                           tCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACATCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAg                                                                                                                                                              >  6:408085/1‑149 (MQ=255)
                                                                                                                                                                                                                             aaCACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCg                                                                                                                                                            <  5:108767/149‑1 (MQ=255)
                                                                                                                                                                                                                               cacGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCCAGTCCtactac                                                                                                                     >  2:267985/1‑186 (MQ=255)
                                                                                                                                                                                                                                   cTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGc                                                                                                                                                                          >  6:69521/1‑129 (MQ=255)
                                                                                                                                                                                                                                   cTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATc                                                                                                                                                      <  5:440110/149‑1 (MQ=255)
                                                                                                                                                                                                                                   cTATGTAGTTTGTTCGGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCCAGTCCTACTACGGGc                                                                                                                 <  4:267981/186‑1 (MQ=255)
                                                                                                                                                                                                                                    tATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATcc                                                                                                                                                     >  5:433978/1‑149 (MQ=255)
                                                                                                                                                                                                                                    tATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATcc                                                                                                                                                     <  6:433978/149‑1 (MQ=255)
                                                                                                                                                                                                                                           tttGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATa                                                                                                                                              <  8:307972/149‑1 (MQ=255)
                                                                                                                                                                                                                                            ttGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGAcc                                                                                                                                                                     >  4:540359/1‑125 (MQ=255)
                                                                                                                                                                                                                                             tGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCCAGTCCTACTACGGGCACCAGTAACGATAACGGTGCAAcc                                                                                         >  3:17271/1‑200 (MQ=255)
                                                                                                                                                                                                                                              gTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACtt                                                                                                                                   >  2:191086/1‑157 (MQ=255)
                                                                                                                                                                                                                                               ttCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCCAGTCCTACTACGGGCACCAGTAACGATAACGGTGCAACCCGCCAGGTTTCATAGCGTCCCAGTAACGTCCCCCAGATCCCATAACCAACAATTGTCGCCACAAACGCCAGATACAt            <  1:143963/275‑1 (MQ=255)
                                                                                                                                                                                                                                                tCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGtt                                                                                                                                                                 <  2:540364/125‑1 (MQ=255)
                                                                                                                                                                                                                                                tCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAAAAATAGTGCCGCACTTGCCAGTCCTACTACGGGCACCAGTAACTATAACGGTGCAAACCGCCAGGTTTCCTAGCGTCCCAGTACCGTCCCCCAGATCCCATAACCAACACTTATCGCCCCAAACGcccagtacacc           >  2:329059/1‑266 (MQ=255)
                                                                                                                                                                                                                                                ataagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGAccc                                                                                                                                                                    >  7:72179/5‑122 (MQ=255)
                                                                                                                                                                                                                                                ataagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGAccc                                                                                                                                                                    <  8:72179/118‑1 (MQ=255)
                                                                                                                                                                                                                                                 cTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCCAGTCCTACTACGGGCACCAGTAACGATAACGGTGCAACCCGcc                                                                                     <  1:17271/200‑1 (MQ=255)
                                                                                                                                                                                                                                                  tGGCTCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCCAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGCTCAACCAACAATAGTCCCGCAc                                                                                                                                     >  4:296578/1‑151 (MQ=255)
                                                                                                                                                                                                                                                  tGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCCa                                                                                                                               <  4:191084/157‑1 (MQ=255)
                                                                                                                                                                                                                                                  aagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGc                                                                                                                                       >  8:460801/3‑149 (MQ=255)
                                                                                                                                                                                                                                                    gCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCAc                                                                                                                                     >  5:124756/1‑149 (MQ=255)
                                                                                                                                                                                                                                                    gCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCAc                                                                                                                                     >  6:138061/1‑149 (MQ=255)
                                                                                                                                                                                                                                                    gCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGccc                                                                                                                                                                                                         >  4:92156/1‑81 (MQ=255)
                                                                                                                                                                                                                                                      cccGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGc                                                                                                                                 <  2:296582/151‑1 (MQ=255)
                                                                                                                                                                                                                                                      cccGAAATCTCGGTGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCCAGTCCTACTACGGGCACCAGTAACGATAACGGTGCAACCCGCCAGGTTTCATAGCGTCCCAGTAACGTCCCCCAGATCCCATAACCAACAATTGTCGCCACAAACGCCAGATACATCAGAGAc     <  2:326861/275‑1 (MQ=255)
                                                                                                                                                                                                                                                        cGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTc                                                                                                                                                                                                     <  2:92157/81‑1 (MQ=255)
                                                                                                                                                                                                                                                        cGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGc                                                                                                                                 >  7:82641/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         gACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCCTGTCCTACTACGGGCACCAGTAACGATAACGGTGCAACCCGCCAGGTTTCATAGCGTCCCAGTAACGTCCCCCAGATTCCCTAACCAACAATTTTCGCCACAAACGCCAGATACATCCGAGAAAAg  >  2:342994/1‑275 (MQ=255)
                                                                                                                                                                                                                                                         gACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCCAGTCCTACTACGGGCACCAGTAACGATAACGGTGCAACCCGCCAGGTTTCATAGCGTCCCAGTAACGTCCCCCAGATCCCATAACCAACAATTGTCGCCCCAAACGCCAGATAAATCAGAGCCaa   >  2:433379/1‑274 (MQ=255)
                                                                                                                                                                                                                                                         gACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTAATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATATTGCCGCACTTGCCAGTCCTACTACGCGCCCCAGTAACGATAACGGTGCAACCACCCAGGTTTCCTCGCTTCCCAGTAACGCCCCCCAGATCCCAAAACCAACCATGTTCGCCACAAACCCCAGCTACATTAGAAACaaa  >  4:101304/1‑274 (MQ=255)
                                                                                                                                                                                                                                                             tctcGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCCAGTCCTACTACGGGCACCAGTAACGATAACGGTGCAACCCGCCAGGGTTCATAGCGTCCCAGTAACGTCCCCCAGATCCCATAACCAACAATTGTCGCCACAAACGCCAGATACAt            <  2:101305/261‑1 (MQ=255)
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                             
CGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCCAGTCCTACTACGGGCACCAGTAACGATAACGGTGCAACCCGCCAGGTTTCATAGCGTCCCAGTAACGTCCCCCAGATCCCATAACCAACAATTGTCGCCACAAACGCCAGATACATCAGAGACAAG  >  NC_000913/1619973‑1620496

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