Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,620,224 C→G 5.4% intergenic (+17/+14) 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,2240CG 5.4% 387.6 / 10.5 111intergenic (+17/+14)marB/eamAperiplasmic mar operon regulator/cysteine and O‑acetyl‑L‑serine efflux system
Reads supporting (aligned to +/- strand):  ref base C (53/50);  new base G (4/2);  total (58/53)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.54e-01

Only 100 of 111 total aligned reads displayed.

GCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCCAGTCCTACTACGGGCACCAGTAACGATAACGGTGCAACCCGCCAGGTTTCATAGCGTCCCAGTAACGTCCCCCAGATCCCATAACCAACAATTGTCGCCACAAACGCCAGATACATCAGAGACAAG  >  NC_000913/1619953‑1620496
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                
gCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTAATTCTTGTGCCAATGTCGTGGTTGTTCCCACATCGCAGGAACACGAAACGTTTGATTTAAATCACATGGATACTGGCAGTGATAAGTCGGATGAGCTCGGCGTTGACTATTATACTAAACCAGCTATGTAGTTTGTTCTGGCCCAGAAAta                                                                                                                                                                                                                                                                               >  1:137956/1‑274 (MQ=255)
                    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                                                                                                                                                                                                                                                                             >  6:355777/1‑144 (MQ=255)
                                                                                                                                     ccAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACAtctc                                                                                                                                                                                                                                                                             <  5:355777/144‑1 (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:223317/149‑1 (MQ=255)
                                                                                                                                                                            cGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCaa                                                                                                                                                                                                                                 <  7:216055/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                                                                                                                                                                                                                                                   <  6:326642/123‑1 (MQ=255)
                                                                                                                                                                                    ttAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAc                                                                                                                                                                                                                                                   >  5:326642/1‑123 (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                                                                                                                                                                                                                      <  5:304374/137‑1 (MQ=255)
                                                                                                                                                                                                   aCTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACAt                                                                                                                                                                                                                      >  6:304374/1‑137 (MQ=255)
                                                                                                                                                                                                       gCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCaa                                                                                                                                                                                                                                 >  6:422984/1‑122 (MQ=255)
                                                                                                                                                                                                       gCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCaa                                                                                                                                                                                                                                 <  7:80992/122‑1 (MQ=255)
                                                                                                                                                                                                       gCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCaa                                                                                                                                                                                                                                 >  8:80992/1‑122 (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                                                                                                                                                                    <  7:230643/149‑1 (MQ=255)
                                                                                                                                                                                                                                         ttatAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGAccc                                                                                                                                                                    >  6:174192/1‑149 (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                                                                                                                                     >  6:138061/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                        gCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCAc                                                                                                                                     >  5:124756/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)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                
GCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCCAGTCCTACTACGGGCACCAGTAACGATAACGGTGCAACCCGCCAGGTTTCATAGCGTCCCAGTAACGTCCCCCAGATCCCATAACCAACAATTGTCGCCACAAACGCCAGATACATCAGAGACAAG  >  NC_000913/1619953‑1620496

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