Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,637,751 G→A L352L (CTG→CTA creC → sensory histidine kinase CreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,637,7510GA97.0% 181.5 / NA 68L352L (CTG→CTAcreCsensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (31/33);  minor base C (2/0);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GCGGCGATTCGTGGAGCGGCGGAAATTTTACGCGAAGGTCCGCCGCCGGAAGTGGTGGCTCGTTTTACTGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGCTGTTGATGTGGCGGCATTATTCCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAA  >  NC_000913/4637509‑4637993
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  
gcggcgATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTat                                                                                                                                                                                                                                            <  1:563694/251‑2 (MQ=255)
          gaggcgctgcggAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGTCAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTTATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAg                                                                                                                                                                                                                                  <  2:355346/243‑1 (MQ=255)
          gTGGCGCGGCGGAAATTTTACGCGCAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGGCGCAACATGCGCGTATGCAGGCATTGGTAGACACGTTACTACGCAAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCAATCTTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGATGGCGGAAAAAGACATCACTCTCCATGTTATGccccctg                                                                                                                                                                                                                                    >  1:254122/1‑245 (MQ=255)
          gTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCGGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTACTTCGCCGCGTCCGCGCAGCGCGCACCGTGGAGTTGGCCGGAAAAAACATCACTCTACATGGTATGCCCActg                                                                                                                                                                                                                                    >  2:503000/1‑247 (MQ=255)
           tggcgcggcgGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCCTTATTTCGCCGCGTCTGCGAAGCGCGCCCCGTGCAGTTGGCAGAAAAACACATCCCTCTACATGttttgccccctgagg                                                                                                                                                                                                                                 >  1:222434/5‑238 (MQ=255)
            ggcgcggcgGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACCTGGTATGCCCACTGAGGt                                                                                                                                                                                                                                >  1:607296/4‑251 (MQ=255)
                cggcggAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGGTaat                                                                                                                                                                                                                            >  2:181235/1‑250 (MQ=255)
                    ggAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGAGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCCGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGACGTTGGGAGAAAAAAAAATCACCCTACATGTTATGACCACTGAAGTTaatgatg                                                                                                                                                                                                                        >  1:330188/1‑246 (MQ=255)
                              cgcgAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGCAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAAcc                                                                                                                                                                                                              >  2:456790/1‑251 (MQ=255)
                                   aGGTCCGCCACCGGAAGTGGTGGCTCGTTTTGCCGACAACATTCTGACGCAGAATGCGCGTATGCAGGCATTGGTAGAAGCGTTACTACGCCAGGCAAGACTGGAGAATCGTCATGAAGGCGGTCTGACTGTTGTTGATGTGGTGGCATTAGTTCGCCGCGTGTGCGAAGCGCGCGCCGTGCGTTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGt                                                                                                                                                                                                         <  1:309990/251‑1 (MQ=255)
                                    ggTCCGCCACCGGAAGTGGTGGCTCGGTTTAACGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGCTAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTTCATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGtt                                                                                                                                                                                                        <  2:375677/251‑1 (MQ=255)
                                         gcaaccgGAAGTGGTTGCTCGTTTTAACGAAAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTTTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTgg                                                                                                                                                                                                   <  2:445172/247‑1 (MQ=255)
                                              cgGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCTTTACttgggccg                                                                                                                                                                                              >  2:248915/1‑244 (MQ=255)
                                                gAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCCGCGAAGCGCGCACCGTGCAGTTGGCAGAAACAAACATCACTCTACATGTTATGCCCCCTGCGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAggg                                                                                                                                                                                            >  2:130333/1‑250 (MQ=255)
                                                 aaGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCGGGCATTGGTAGAAACGTTACTCCTCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTGCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGGCAGCGAAGCGCGCTCCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAAGGTTGCTGCTGAACCGTCGTTATTGGAGCAGgcg                                                                                                                                                                                           <  2:557048/251‑1 (MQ=255)
                                                      gtggCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCCCTGAGGTTAATGTTGCTGCTGCACCGGCGTTACTGGAGCAGGCGCTggg                                                                                                                                                                                      >  1:388765/1‑251 (MQ=255)
                                                       tggCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCAACGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACAGGCGTTACTGGAGCCGGCGCTgggg                                                                                                                                                                                     >  2:466440/1‑250 (MQ=255)
                                                              ttttatcgaaaactttCTGACGCAATGTGCGCGAATTCAGGCATTGGTAGAAACGTTACTTCGCAAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCCTCCGCGGAGCGCGCACCGTGCCGGTTGCAGAAAAAAACATCACCCTACCTGTTCAGCCCACTGCGGTTACTTTTACTGCCGATCCGGCGTTACTGGATAAGGCGCTgggggatttct                                                                                                                                                                              >  2:361622/15‑243 (MQ=255)
                                                              ttttACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGCAACGTTACTACGCCAGGCAAGACTGGAGACTCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCACCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTTCATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGTGAATTTAc                                                                                                                                                                              >  2:140084/1‑251 (MQ=255)
                                                                ttACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCCCTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGCGCAGGCGCTGGGGAATTTACTg                                                                                                                                                                            >  1:238888/1‑251 (MQ=255)
                                                                ttACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTAAGGAAGTCGTAATGACTGTTGTTGATGAGGCGGCATTATTTCGCCGCGTAAGCGGAGAGCGCACCGTGCAGTTGGAAGAAACAAACACAACTCTACATGTTTTGCCCCCTGCGGTTAAAGTTGCTGCTGAACAGGCGTTACTGGAGCAGGGGCTggggagattaacg                                                                                                                                                                            >  2:226168/1‑241 (MQ=255)
                                                                   ccgACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGAt                                                                                                                                                                         <  1:456790/249‑1 (MQ=255)
                                                                       caacaTTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACg                                                                                                                                                                     >  1:602001/1‑251 (MQ=255)
                                                                         acaTTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGcc                                                                                                                                                                   >  1:65720/1‑251 (MQ=255)
                                                                                gACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTGGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTTAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCg                                                                                                                                                               <  1:421242/248‑1 (MQ=255)
                                                                                 aCGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTTGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGAttt                                                                                                                                                           <  1:302929/251‑1 (MQ=255)
                                                                                 aCGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTTGTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGGATTTACTGGTTAACGCCTTCGGttt                                                                                                                                                           >  2:231555/1‑251 (MQ=255)
                                                                                  cGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCTTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCAGTATTTCGCCGCGTCAGCGAAGCGCGGACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACGGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGTAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGAtttt                                                                                                                                                          <  2:275530/251‑1 (MQ=255)
                                                                                      aaaTGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTAccc                                                                                                                                                      <  1:231555/251‑1 (MQ=255)
                                                                                                aTGCAGGCATTGGTAGAAACGTTTCTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCgg                                                                                                                                            <  1:521880/251‑1 (MQ=255)
                                                                                                        aTTGGTAGAAACGTTACTACGCCAGGCAAGACGGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGGTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCGGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCAtaa                                                                                                                                    <  2:85916/251‑3 (MQ=255)
                                                                                                              aGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTa                                                                                                                               <  2:447361/250‑1 (MQ=255)
                                                                                                                aaaCGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGAGCAGCGCGCACCGGGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGGTAATGTTGCTGCTGAACCGGCGGTACTGGCGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAgc                                                                                                                            >  1:280968/1‑251 (MQ=255)
                                                                                                                 aaCGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACGGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACCGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAgcg                                                                                                                           <  1:627929/251‑1 (MQ=255)
                                                                                                                   cGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGcc                                                                                                                         <  1:181235/251‑1 (MQ=255)
                                                                                                                            gCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGAt                                                                                                                <  1:46855/251‑1 (MQ=255)
                                                                                                                                          gcgaATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGGTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGGATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGGGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTcgc                                                                                                  >  2:17545/3‑249 (MQ=255)
                                                                                                                                              aTCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCTCCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTcgccctt                                                                                              <  1:396591/251‑7 (MQ=255)
                                                                                                                                              aTCGTCAGAAAGTCGATCTGACTGTTGTTGATGTGGCGGCCTTATTTCGCCGCGTCAGCGAAGCGGGAACCGTGCAGGTGACCGAAAAAAACATCACTCTACATTTTATGCCCACTGCGGTTAATGTTGCTGCTGAACCGGAGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGCCGAATAACGCTCAGAGCCGTAGTGGCTCAGGCACACGTccccctt                                                                                              >  2:379221/1‑245 (MQ=255)
                                                                                                                                                       aaGTCGTTCTGACTGTTGTTGATGTGGCGGCCTTATTTCGCCGCGTCAGCGCAAGGCGCACCGTGCAGTTGGCAGCACAAAACATACCTCTACAAGTTATGCCCACTGAGGTTAATGTTGCTGCTTCACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACCCTAAGCCCCGGAGTGGGTCAGGAACAGGTGGCCCTTAAGGTGCTg                                                                                     >  2:81394/1‑251 (MQ=255)
                                                                                                                                                           cGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTTCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGAt                                                                                  <  1:313093/250‑1 (MQ=255)
                                                                                                                                                                        gTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACGCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTa                                                                     <  2:194143/250‑1 (MQ=255)
                                                                                                                                                                          tGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTAtt                                                                   <  2:248124/250‑1 (MQ=255)
                                                                                                                                                                              gtgGCGGCATTCTTTAGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGGGGATCAGGAACACGTCGCCCTTAAGGGGCTGGATACCGGTAGTGGTATTCCTg                                                               >  1:483271/1‑250 (MQ=255)
                                                                                                                                                                              gtgGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGGGGTTAATGTTGCTGCTGAACCGGCGATACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGAACAGGAACACGTCGCCCTTAAGGAGCTGGATACCGGTAGGGGTATTCCTgg                                                              >  1:623454/1‑250 (MQ=255)
                                                                                                                                                                              gtgGCGGCAATATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGGAGTAAAAAACATCACTCTCCATGTTATGCCCCCTGAGGTTACTGTTGCTGCTGCACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGAATAACGCTAAGCGCCGAAGGGGGTCAGGAACACGTCGCCCATAAGGTGCTGGATACCGGTAGTGGGATTcatg                                                               >  1:363722/1‑247 (MQ=255)
                                                                                                                                                                                        tattTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCTCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCt                                                     <  1:259783/250‑1 (MQ=255)
                                                                                                                                                                                        tattTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCt                                                     <  1:536973/250‑1 (MQ=255)
                                                                                                                                                                                            tcgccgcGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCa                                                 <  1:248915/249‑1 (MQ=255)
                                                                                                                                                                                             cgccgcGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCAc                                                <  1:408735/250‑1 (MQ=255)
                                                                                                                                                                                               ccgcGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAACACCTCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGCTTTTACCCCCGAGAGCGGTCGCATAACGCTACGCGCCGAAGTGGCTCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACga                                              >  2:57617/1‑249 (MQ=255)
                                                                                                                                                                                               ccgcGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTa                                             <  1:348878/251‑1 (MQ=255)
                                                                                                                                                                                               ccgcGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGGAGTGGCATTCCTGACTCCGCGCTGTCACGt                                              >  2:244836/1‑250 (MQ=255)
                                                                                                                                                                                                cgcgTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTAt                                            <  1:180501/251‑1 (MQ=255)
                                                                                                                                                                                                     cAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTAttttt                                        <  1:221011/250‑1 (MQ=255)
                                                                                                                                                                                                        cGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGACGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCGGACTACGCGCTGTCACGTATTTTTGAAc                                    <  2:222434/251‑1 (MQ=255)
                                                                                                                                                                                                        cGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCAGCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAAc                                    <  1:590521/251‑1 (MQ=255)
                                                                                                                                                                                                              cgcgcACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCttt                               >  1:265195/1‑251 (MQ=255)
                                                                                                                                                                                                                 gcACCGTGCAGTTGGCAGAAAACCACCTCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGACGGGGATCAGGACCACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGACAGCTTTTAc                            >  2:397785/1‑251 (MQ=255)
                                                                                                                                                                                                                 gcACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTAc                            >  1:391613/1‑251 (MQ=255)
                                                                                                                                                                                                                   aCCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCCCTGAGGTTAATGTTGCTGCTGAACCGACGTTACTGGAGCAGGCGCTGGGGAATTTCCTGGATAACGCCATCGATTTTAACCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACAAGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTActc                          >  1:289193/1‑251 (MQ=255)
                                                                                                                                                                                                                    ccGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTActct                         <  2:388765/251‑1 (MQ=255)
                                                                                                                                                                                                                          aGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGCACGCTTTTACTCTTTGccg                   <  1:473832/251‑2 (MQ=255)
                                                                                                                                                                                                                                cAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCa             <  2:375002/251‑1 (MQ=255)
                                                                                                                                                                                                                                        aacaTCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGc      >  1:329851/4‑250 (MQ=255)
                                                                                                                                                                                                                                         acaTCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCa     >  1:620401/3‑250 (MQ=255)
                                                                                                                                                                                                                                          caTCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCaa    >  1:218688/2‑250 (MQ=255)
                                                                                                                                                                                                                                           aTCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCaaaa  <  1:130333/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  
GCGGCGATTCGTGGAGCGGCGGAAATTTTACGCGAAGGTCCGCCGCCGGAAGTGGTGGCTCGTTTTACTGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGCTGTTGATGTGGCGGCATTATTCCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAA  >  NC_000913/4637509‑4637993

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