Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,620,227 C→G 6.3% 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 6.3% 402.9 / 14.6 114intergenic (+20/+11)marB/eamAperiplasmic mar operon regulator/cysteine and O‑acetyl‑L‑serine efflux system
Reads supporting (aligned to +/- strand):  ref base C (49/57);  new base G (4/3);  total (54/60)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-02

Only 100 of 126 total aligned reads displayed.

CTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCCAGTCCTACTACGGGCACCAGTAACGATAACGGTGCAACCCGCCAGGTTTCATAGCGTCCCAGTAACGTCCCCCAGATCCCATAACCAACAATTGTCGCCACAAACGCCAGATACATCAGAGACAAGATGGT  >  NC_000913/1619954‑1620501
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                  
cTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCGTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGACAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACAtct                                                                                                                                                                                                                                                                                   <  2:92802/275‑1 (MQ=255)
                  aCGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACGAACAACGCAGCCAGTTGGTACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAGATCACATGGGTACTGGAAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGCGCTTATTAACTTcc                                                                                                                                                                                                                                                                 <  2:377616/275‑1 (MQ=255)
                                             tttCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCac                                                                                                                                                                                                                                      >  2:243272/1‑274 (MQ=255)
                                             tttCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCaa                                                                                                                                                                                                                                      >  2:159134/1‑275 (MQ=255)
                                                cATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATGGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGcc                                                                                                                                                                                                                                   <  2:270924/275‑1 (MQ=255)
                                                cATCCGCAATAGCAGCTGCGCTGATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATGTAAATCACATGGGTACTGGCAGTGATAAGTCGGAGGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGcc                                                                                                                                                                                                                                   <  1:151713/275‑1 (MQ=255)
                                                     gctaTAGTAGCTGCGCTTATTCTCTTTTCCGGGAAGGGCGTTAATGAACAAACCCCGCAGCCAGTTGTTACTTCGTGTGCCAATGTCTTGGTTGTTCCCCCAGCGCAGGAACACCCAACGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATa                                                                                                                                                                                                                              <  2:529075/272‑1 (MQ=255)
                                                     gCAATAGCAGCTGCGCTTATTTTCTTTTCCGCGACGGGCGTTGACGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATa                                                                                                                                                                                                                              <  4:451568/275‑1 (MQ=255)
                                                     gCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATa                                                                                                                                                                                                                              <  2:539185/275‑1 (MQ=255)
                                                       aaatgcccggtcgcTTTTTCTCTTTTCCGCGCAGGGCGTTTCGGAACAAACCACCCCGCCAGTTGTTACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACa                                                                                                                                                                                                                            <  4:531899/264‑1 (MQ=255)
                                                       aaTAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCa                                                                                                                                                                                                                                       >  2:492143/1‑264 (MQ=255)
                                                       aaTAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACa                                                                                                                                                                                                                            <  2:57944/275‑1 (MQ=255)
                                                                    cTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGc                                                                                                                                                                                                                >  2:365157/1‑274 (MQ=255)
                                                                       tttttcgtgtCCGCGCAGGGCGTTGCGGAAACAACCACTAAGGCTGTTGTTACTTATTGGTCTAATGCCGTGGGTGTTGCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACGGGCAGTAATAAGGCGGACGCGTGCGGTGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACAGCTCGGGACTTATGAACTTCCCACCTTCACCGCTGTACGCCACCGCAAGCCAAATACATTGATATACAGCCCgg                                                                                                                                                                                                            <  4:243978/266‑1 (MQ=255)
                                                                           tctTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTc                                                                                                                                                                                                          <  2:118886/273‑1 (MQ=255)
                                                                                         ggCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgcac                                                                                                                                                                                          <  2:498923/275‑1 (MQ=255)
                                                                                                                                             tGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGAGGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGc                                                                                                                                      <  4:109880/275‑1 (MQ=255)
                                                                                                                                                     cccccATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCAc                                                                                                                                          >  2:67747/1‑263 (MQ=255)
                                                                                                                                                       cccATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACAtc                                                                                                                                                                                                                                                                                    >  7:925842/1‑123 (MQ=255)
                                                                                                                                                         cATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGc                                                                                                                                      <  4:67746/263‑1 (MQ=255)
                                                                                                                                                                         accGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCCAGTCCTACTACGGGCACCAGTAACGAt                                                                                                          <  4:365156/275‑1 (MQ=255)
                                                                                                                                                                          ccGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCGGGGcc                                                                                                                                                                                                                                                                            >  8:1076807/1‑111 (MQ=255)
                                                                                                                                                                           cGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGc                                                                                                                                                                                                                                                                             <  6:849734/110‑1 (MQ=255)
                                                                                                                                                                           cGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGc                                                                                                                                                                                                                                                                             >  5:849734/1‑110 (MQ=255)
                                                                                                                                                                              ttGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCgg                                                                                                                                                                                                                                                                                <  8:1142620/104‑1 (MQ=255)
                                                                                                                                                                              ttGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCgg                                                                                                                                                                                                                                                                                >  7:1142620/1‑104 (MQ=255)
                                                                                                                                                                              tcgATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCg                                                                                                                                                                                                                                                                                 >  6:993607/3‑103 (MQ=255)
                                                                                                                                                                              tcgATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCg                                                                                                                                                                                                                                                                                 <  5:993607/101‑1 (MQ=255)
                                                                                                                                                                                  tttAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCt                                                                                                                                                                                                                                                                            >  5:53966/1‑104 (MQ=255)
                                                                                                                                                                                  tttAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCt                                                                                                                                                                                                                                                                            <  6:53966/104‑1 (MQ=255)
                                                                                                                                                                                   ttAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACAtct                                                                                                                                                                                                                                                                                   <  7:486566/96‑1 (MQ=255)
                                                                                                                                                                                   ttAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACAtct                                                                                                                                                                                                                                                                                   <  7:84057/96‑1 (MQ=255)
                                                                                                                                                                                   ttAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACtt                                                                                                                                                                                                                                                                   >  7:704321/1‑112 (MQ=255)
                                                                                                                                                                                   ttAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACtt                                                                                                                                                                                                                                                                   <  8:704321/112‑1 (MQ=255)
                                                                                                                                                                                   ttAAATCACATGGGTACAGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACAtct                                                                                                                                                                                                                                                                                   >  8:84057/1‑96 (MQ=255)
                                                                                                                                                                                    tAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACAtctc                                                                                                                                                                                                                                                                                  >  8:888207/1‑96 (MQ=255)
                                                                                                                                                                                    tAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACAtctc                                                                                                                                                                                                                                                                                  <  7:888207/96‑1 (MQ=255)
                                                                                                                                                                                    tAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATa                                                                                                                                                                                                                              >  6:80701/1‑148 (MQ=255)
                                                                                                                                                                                       aTCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACAt                                                                                                                                                                                                                           <  5:80701/148‑1 (MQ=255)
                                                                                                                                                                                                gTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCg                                                                                                                                                                                                                                                                                 <  8:1047148/85‑1 (MQ=255)
                                                                                                                                                                                                gTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCg                                                                                                                                                                                                                                                                                 >  7:1047148/1‑85 (MQ=255)
                                                                                                                                                                                                    tGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGccc                                                                                                                                                                                                              >  6:167686/1‑148 (MQ=255)
                                                                                                                                                                                                    tGGCAGTGATAAGCCGGATTCGCTCCGCGTGCCCTATTATAACCAACACGCTTTGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTGTACCGCTTTACGACACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCCAGTCCTACTACGGGCACCAGTAACGATAACGGTGCAACCCGCCAGGTTTCATAg                                                                               <  4:243271/275‑1 (MQ=255)
                                                                                                                                                                                                      gCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCgg                                                                                                                                                                                                                                                                                <  8:241201/80‑1 (MQ=255)
                                                                                                                                                                                                      gCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCACCAATAGTGCCGCACTTGCCCTTCCTCCTACGGACCCCAGTAACCAAACCGGGGAACCCCGCCAGGTTTCAGAGc                                                                              >  4:147727/1‑274 (MQ=255)
                                                                                                                                                                                                         gTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACCATAGTGCCGCACTTGCCAGTCCTACTACGGGCACCAGTAACGATAACGGTGCAACCCGCCAGGTTTCATAGCGTcc                                                                          >  4:219187/1‑275 (MQ=255)
                                                                                                                                                                                                          tGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCACATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCCCACTTGCCAGTCCTACTCCGGGCACCAGTAACGATAACGGTGCACCCCGCCAGGTTTCATAGCGAccc                                                                         >  2:5839/1‑275 (MQ=255)
                                                                                                                                                                                                          tGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCCAGTCCTACTACGGGCACCAGTAACGATAACGGt                                                                                                    >  4:263660/1‑248 (MQ=255)
                                                                                                                                                                                                              aaGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCCAGTCCTACTACGGGCACCAGTAACGATAACGGTGCaa                                                                                                <  2:263661/248‑1 (MQ=255)
                                                                                                                                                                                                                gTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAg                                                                                                                                                                                                  >  7:126282/1‑148 (MQ=255)
                                                                                                                                                                                                                gTCGGATGCGCTCAGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCCAGTCCTACTACGGGCACCAGTAACGATAACGGTTCAACCCGCCAGGTTTCATTGCGTCCCAGTaa                                                                    >  2:197214/1‑274 (MQ=255)
                                                                                                                                                                                                                 tCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCtt                                                                                                                                                                                                                                                                           >  6:459475/1‑74 (MQ=255)
                                                                                                                                                                                                                 tCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgc                                                                                                                                                                                                 <  8:126282/148‑1 (MQ=255)
                                                                                                                                                                                                                  cGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCt                                                                                                                                                                                                                                                                            <  8:670865/72‑1 (MQ=255)
                                                                                                                                                                                                                        cgcTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgcac                                                                                                                                                                                          <  6:619049/148‑1 (MQ=255)
                                                                                                                                                                                                                               cgggCCATAATATAATCACCCCGCTATGTGGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTCCCCACCTTTACCAGTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGTCCCGTTAGTCGTTCATCCAACAATAGTGCCGCACTTGCCAGTCCTACTACGGGCACCAGTAACGATAACGGTGCAACCCGCCAGGTTTCATAGCGTCCCAGTAACGTCCCCCAGATCCCa                                                    <  4:159133/272‑1 (MQ=255)
                                                                                                                                                                                                                                gagaaagattatATTCAAAAAGCTATGTAGTTTGTTCTGGCCCCGACACCTCGGGGCTTATTAACTTCCAACCTGTACCGCTTGACGCCACCGCAAGCCAAATACATGGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCCAGTCCTACTACGGGCACCAGTAACGATAACGGTGCAACCCGCCAGGTTTCATAGCGTCCCAGTAACGTCCCCCAGATCCCa                                                    <  4:510552/267‑1 (MQ=255)
                                                                                                                                                                                                                                       attatAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCttt                                                                                                                                                                                                                                                  >  5:566650/1‑77 (MQ=255)
                                                                                                                                                                                                                                       attatAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTa                                                                                                                                                                                                                                                 >  5:54785/1‑78 (MQ=255)
                                                                                                                                                                                                                                       attatAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTa                                                                                                                                                                                                                                                 <  6:54785/78‑1 (MQ=255)
                                                                                                                                                                                                                                           taATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGa                                                                                                                                                                            >  8:27504/1‑143 (MQ=255)
                                                                                                                                                                                                                                           taATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGa                                                                                                                                                                            <  7:27504/143‑1 (MQ=255)
                                                                                                                                                                                                                                            aaTCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgc                                                                                                                                                                                                 <  5:672783/121‑1 (MQ=255)
                                                                                                                                                                                                                                            aaTCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgc                                                                                                                                                                                                 >  6:672783/1‑121 (MQ=255)
                                                                                                                                                                                                                                              tCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCCAGTCCTACTACGGGCACCAGTAACGATAACGGTGCAACCCGCCAGGTTTCATAGCGTCCCAGTAACGTCCCCCAGATCCCATAACCAACACTTGTc                                     >  2:114455/1‑275 (MQ=255)
                                                                                                                                                                                                                                                aaCACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGc                                                                                                                                               >  4:71208/1‑167 (MQ=255)
                                                                                                                                                                                                                                                 acacGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACa                                                                                                                                                                                                                            >  5:402387/1‑89 (MQ=255)
                                                                                                                                                                                                                                                 acacGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCg                                                                                                                                                                                                             >  6:956104/1‑104 (MQ=255)
                                                                                                                                                                                                                                                  cacGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAg                                                                                                                                                                                                                                     >  4:184088/1‑79 (MQ=255)
                                                                                                                                                                                                                                                  cacGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAAtt                                                                                                                                                                                 >  4:221603/1‑131 (MQ=255)
                                                                                                                                                                                                                                                  cacGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCCAGTCCTACTACGGGCACCAGTTACGATAACGGTGCAACCCGCCAGGTTTCATAGCGTCCCATTAACGTCCCCCAGAt                                                        >  2:510553/1‑252 (MQ=255)
                                                                                                                                                                                                                                                    ctcTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCCAGTCCTACTACGGGCACCAGTAACGATaa                                                                                                        <  1:131930/200‑1 (MQ=255)
                                                                                                                                                                                                                                                    cGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCa                                                                                                                                           <  2:71209/167‑1 (MQ=255)
                                                                                                                                                                                                                                                      cTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCaa                                                                                                                                                                                                                                 <  2:184089/79‑1 (MQ=255)
                                                                                                                                                                                                                                                      cTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAg                                                                                                                                                                             <  2:221604/131‑1 (MQ=255)
                                                                                                                                                                                                                                                      cTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCAt                                                                                                                                                            >  5:292886/1‑148 (MQ=255)
                                                                                                                                                                                                                                                       tATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACa                                                                                                                                                                                                                            >  7:133591/1‑83 (MQ=255)
                                                                                                                                                                                                                                                       tATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACa                                                                                                                                                                                                                            <  8:133591/83‑1 (MQ=255)
                                                                                                                                                                                                                                                            aGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAc                                                                                                                                                                                                                   >  7:771330/1‑87 (MQ=255)
                                                                                                                                                                                                                                                              tttGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGa                                                                                                                                                                                                   <  6:397990/101‑1 (MQ=255)
                                                                                                                                                                                                                                                              tttGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGa                                                                                                                                                                                                   >  5:397990/1‑101 (MQ=255)
                                                                                                                                                                                                                                                              tttGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTa                                                                                                                                                                     <  7:1076587/131‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ttCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAATGCCGCACTTGCCAGTCCTACTACGGGCACCAGTAACGATAACGGTGCAACCCGCCAGGTTTCATAGCGTCCCAGTAACGTCCCCCAGATCCCATAACCAACAATTGTCGCCACAAACGCCAGATACAt                 <  2:219188/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                   tCTGGCCCCGACATCTCGGGGCTTATTCACATCCCACCTTTACCGCTTTACGCCACCGCa                                                                                                                                                                                                                                       >  8:700588/1‑60 (MQ=255)
                                                                                                                                                                                                                                                                   tCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCa                                                                                                                                                                                                                                       <  7:700588/60‑1 (MQ=255)
                                                                                                                                                                                                                                                                   tCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATAc                                                                                                                                                                                                                             >  5:950417/1‑70 (MQ=255)
                                                                                                                                                                                                                                                                   tCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATTAGCACCGCCCCTAAAAATTGCAGACCCGTTAAGCGTTCATCca                                                                                                                                                         >  1:20312/1‑138 (MQ=255)
                                                                                                                                                                                                                                                                     tGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgca                                                                                                                                                                                                >  7:149792/1‑97 (MQ=255)
                                                                                                                                                                                                                                                                      agCCCTGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTg                                                                                                                                                                                >  8:430799/2‑112 (MQ=255)
                                                                                                                                                                                                                                                                       gCCCCGACATCTCGCGGCTTATTAACTTCCCACCTTTACCGCTTTAGTCCACCGCAAGCCAAAAACATTGATATACAGCCCGGTCATAATGAACACCGCACCTAAAAAGTGCAGACCCGTTAAGCGTTCATCcaacaa                                                                                                                                                     <  3:20312/138‑1 (MQ=255)
                                                                                                                                                                                                                                                                        ccccGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAAtt                                                                                                                                                                                 >  8:833278/1‑109 (MQ=255)
                                                                                                                                                                                                                                                                        ccccGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAg                                                                                                                                                                   <  5:420882/123‑1 (MQ=255)
                                                                                                                                                                                                                                                                        ccccGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAg                                                                                                                                                                   >  6:420882/1‑123 (MQ=255)
                                                                                                                                                                                                                                                                          ccGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGa                                                                                                                                                                                                   >  7:733627/1‑89 (MQ=255)
                                                                                                                                                                                                                                                                          ccGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGa                                                                                                                                                                                                   <  8:733627/89‑1 (MQ=255)
                                                                                                                                                                                                                                                                            gACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCCAGTCCTACTACGGGCACCAGTAACGATAACGGTGCAACCCGCCAGGTTTCATAGCGTCCCAGTAACGTCCCCCAGATCCCATAACCAACACTTTTCGCCACAAACGCCAGAACAATCAGAAACac        >  4:401196/1‑273 (MQ=255)
                                                                                                                                                                                                                                                                               aTCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCCAg                                                                                                                                   <  6:292886/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 ctcGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGa                                                                                                                                                                            <  8:352384/105‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 ctcGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGa                                                                                                                                                                            >  7:352384/1‑105 (MQ=255)
                                                                                                                                                                                                                                                                                 ctcGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCCAGTCCTACTACGGGCACCAGTAACGATAACGGTGCAACCCGCCAGGTTTCATAGCGTCCCAGTAACGTCCCCCAGATCCCATAACCAACAATTGTCGCCACAAACGCCAGAAACATCAGAGACAAGAtggt  >  1:131391/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                  
CTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCCAGTCCTACTACGGGCACCAGTAACGATAACGGTGCAACCCGCCAGGTTTCATAGCGTCCCAGTAACGTCCCCCAGATCCCATAACCAACAATTGTCGCCACAAACGCCAGATACATCAGAGACAAGATGGT  >  NC_000913/1619954‑1620501

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