Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,624,355 T→C T145T (ACA→ACG ytjB ← putative inner membrane protein Smp

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,624,3550TC96.4% 148.4 / ‑2.2 57T145T (ACA→ACGytjBputative inner membrane protein Smp
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (28/26);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.11e-02

CTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCG  >  NC_000913/4624127‑4624592
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              
cTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGCTTTCTGCCAGCGGGTACGTTTTCCCTGTTGCAGCGGGCGGGGCAGCACTACGCCGATTGCCAGTTAAAGCAGCCACATCAGGCGGCAAAAGTTAGTGGTGTTATCCCCCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGGCGGGCGGCAGGCGCAGATAGcccggcg                                                                                                                                                                                                                          >  2:99725/1‑246 (MQ=255)
  ttCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCgg                                                                                                                                                                                                                         <  1:336341/250‑1 (MQ=255)
   tCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAAACGGAAGGGCGGTTGCTGCCAGCGGGGACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGGGAAAGCAGCAACATCAGGCGCAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGGGGGCGTATGGGGGTCGGGGGGCAGGCGCAggtaagcgggcggtg                                                                                                                                                                                                                       >  1:3854/1‑237 (MQ=255)
          gtaATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTTACAGGGTGGGCGGTTTCTTCCCCCGGGGACCTTTTCCCTGTAGCAGCGTGCGGGTCCGCCATACGCCGGTTGCCCGGGGAACAAGCAACATCCGGCGGAAAAAGTTCATGGGGTTATCCACCTGGTGGGGTTCGGTGGGGAGGGTATAGGTGGCGTGCGTCAGgcgccggttacccgacgggcccttttt                                                                                                                                                                                                                >  2:159467/1‑228 (MQ=255)
               gggTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGAGTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGtttttcccgg                                                                                                                                                                                                           <  1:428319/251‑6 (MQ=255)
                          ttttCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTggttc                                                                                                                                                                                                 <  2:401757/250‑4 (MQ=255)
                          ttttCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGGGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGGAAAATGTTAGTGGTGTTATCCACCTGTTTGGCTTCGGTGGCGGGCGTATGGGGGTCGGGCGTCAGGCGCAGATAAACGAGCGGTCCGTTTTTCCCGGCCATTGGTTCg                                                                                                                                                                                                >  1:391370/1‑251 (MQ=255)
                              ctcGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGGACGTTTACCCTGTTGCAGCGTGCGGGTCAGCACTACGACGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCCCCTGTTGGGCTTCGGTGGCGGGCGGATGGGGGTCGGGCGTCCGGCGCAGATAGCCGAGCGGTCCGTTTTTTCCGGCCATTGGTTCGAcc                                                                                                                                                                                             >  1:433673/1‑249 (MQ=255)
                               tcGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTTCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGGGCGGGTCAGCACTACGCCGCTTGCCCGTGAAAACAGCAACATCCGGCGTAAAATGTTAGTGGTGTTATTCCACTGTTGGGCTTCGGTGGCGGGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGGCCGGTTTTCCCGGCAATTGGTTCGTCAAt                                                                                                                                                                                           >  1:439434/1‑251 (MQ=255)
                                   tttCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGGAGGGCCATTGCTGCCAGCGGGTAAGTTTACCCTGTAGCCACGGGCGGGGCCGCCCTACGCCCGTTGCCAGTGAAAGCCGCAACATCCGGCGGTAAATTTTAGTGGGGTTATCCCCCTGTTGGGGTTCCGTGGCGCGCGGATTGGGGTCGAGCGGCAGGCGCAGGTAACCGCGCGGGCCGttttttccggcccatggtttgcacctctgt                                                                                                                                                                                       >  2:480212/1‑226 (MQ=255)
                                       ttcctcttccGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTTCCCTGTAGCCGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGGAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGGGCGTATGGGGGTCGAGCGTCAGGCGCAGATAGCCGGGCGGTCCGTTTTTTCCGGCCATTGGTTCGACAATCTTCTGGt                                                                                                                                                                                   >  1:355266/1‑251 (MQ=255)
                                        tcctcttccGGCCCCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGt                                                                                                                                                                                   <  2:13241/250‑1 (MQ=255)
                                        tcctcttccGGCCCCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGt                                                                                                                                                                                   <  2:434238/250‑1 (MQ=255)
                                        tcctcttccGGCACCGGTTTGCTGGCGTTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCTAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATctgctg                                                                                                                                                                                     <  1:440692/248‑1 (MQ=255)
                                        tcctcctccGGCACCGGTTTGCTGGCGGGTAACAGGAAGGGCCATTGCTGCCAGCGGGTACGTTTACGCTGGAGCAGCGTGCGGGTCAGACCTACGCCGATTGCCAGTGAAAGCAGCAACAGCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTGTTCCCGGCAATTGGTTCGACAATCTGCTgg                                                                                                                                                                                    <  2:166404/249‑1 (MQ=255)
                                              tccGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAAt                                                                                                                                                                             <  1:78316/250‑1 (MQ=255)
                                                 ggCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGGACGTTTTCCCTGTAGCAGCGGGCGGGGCAGCCCTACGCCGGTTGCCCGTGGAAGCAGCAACCTCAGGCGGAAAATGTTAGTGGGGTTATCCACCTGTTGGGCTTCGGTGGCGGGCGTGTGGGGGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCAGCTTTTTCCGGCAATTGGGTCGCCAATCTGCTGGTTAAAATAgcc                                                                                                                                                                         >  2:241672/1‑251 (MQ=255)
                                                      cGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTTGCGAGCGTATGTGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAgccgcccg                                                                                                                                                                    <  1:285257/251‑3 (MQ=255)
                                                        gTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGGACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTAAGCCGATTGCCCGTGTAAGAAGCCACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTTGGGGGCGGGCGTATGGGTGTCGGGCGTCAGGCGCAGGTAGCCGAGAGGTCCGTTTTTCCCGGCAATTGGTTCGCCAAATTGCTGGTTAAACTAcccgcccgcc                                                                                                                                                                  >  1:437045/1‑247 (MQ=255)
                                                                  ggTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATCGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATATCCGCCCGCTTTTTTACCGt                                                                                                                                                        >  2:482843/1‑251 (MQ=255)
                                                                                  gATTGCTGCCAGCGGGTACGTTTACCCTTTAGCAGCGTTCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGGGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTc                                                                                                                                        <  2:355266/251‑1 (MQ=255)
                                                                                  gATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTc                                                                                                                                        >  1:165367/1‑251 (MQ=255)
                                                                                  gATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTc                                                                                                                                        <  1:109060/251‑1 (MQ=255)
                                                                                   aTTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGtttttcccg                                                                                                                                                                                                            >  1:202075/1‑178 (MQ=255)
                                                                                   aTTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTcg                                                                                                                                       >  1:261567/1‑251 (MQ=255)
                                                                                   aTTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGGGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGtttttcccg                                                                                                                                                                                                            <  2:202075/182‑5 (MQ=255)
                                                                                       ctgcCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCGGGCGCAGATAGCCGATCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCgcgc                                                                                                                                    >  1:424883/1‑250 (MQ=255)
                                                                                               gggTCCGTTTCCCCTGTAGCAGCGTGCGGGGCAGCACTCAGCCGATCGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGGCGAGCGTCAGGCGCAGATATCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGGCGCGCACTTCGACg                                                                                                                           <  2:179621/251‑1 (MQ=255)
                                                                                                        tACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGgc                                                                                                                                                                                                                                         <  2:177587/132‑1 (MQ=255)
                                                                                                        tACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGgc                                                                                                                                                                                                                                         >  1:177587/1‑132 (MQ=255)
                                                                                                         aCCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGcc                                                                                                                  >  1:68296/1‑249 (MQ=255)
                                                                                                                agcagcGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACgggc                                                                                                         <  1:66029/251‑4 (MQ=255)
                                                                                                                  cagcGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCg                                                                                                        >  2:140155/1‑250 (MQ=255)
                                                                                                                     cGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATaa                                                                                                    <  2:391370/251‑1 (MQ=255)
                                                                                                                              cAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCtt                                                                                           >  1:237649/1‑251 (MQ=255)
                                                                                                                                     tacgccGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCttttttcc                                                                                                                                                           >  2:360671/5‑180 (MQ=255)
                                                                                                                                     tacgccGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTAc                                                                                                                                                           <  1:360671/178‑1 (MQ=255)
                                                                                                                                      acgccGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGGCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTc                                                                                    >  2:439504/4‑251 (MQ=255)
                                                                                                                                                     gaaagcagcAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTAGGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTGTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTc                                                                     <  1:47798/248‑1 (MQ=255)
                                                                                                                                                      aaggcagcaACATCAGGCGAAAAATGTTAGTGGTGTTATCCACCTGTTGGACTTCGGAGGCGAGCGTATTGGAGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGATTTTCCCGGCAATTAGTTCGACAATCTGCTGGTTAAAATATCCACCCGCGTTTTTACCGTCGATCGCCACCCGGTCGCGCACTTCGCGACTTTCGCCCGAGCGGGCGATAAGATCACCTTGtacgaaataccagcccgactcg                                                                    >  2:72363/4‑230 (MQ=255)
                                                                                                                                                             gcaACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGCGGAt                                                               >  1:358684/1‑250 (MQ=255)
                                                                                                                                                                 cATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTAGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGc                                                          >  1:35358/1‑251 (MQ=255)
                                                                                                                                                                    cAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCttt                                                       <  1:391751/251‑1 (MQ=255)
                                                                                                                                                                         gTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATcc                                                  >  1:341866/1‑251 (MQ=255)
                                                                                                                                                                                gTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGTCAGAACGGGCGATAAGATCGCTTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAAc                                            <  1:406633/250‑1 (MQ=255)
                                                                                                                                                                                    ggggtgTTATCCGCCTGTTGGGCTTCGGTGGCGGGCGTTTGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGTTCCGTTTTCCCCGGCAATTGGTTCGACAACCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATc                                       <  2:261567/249‑1 (MQ=255)
                                                                                                                                                                                           tATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGAt                                 <  1:480212/250‑1 (MQ=255)
                                                                                                                                                                                           tATCCACCGGTTGGGCTCCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCCCTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGAt                                 <  2:165367/250‑1 (MQ=255)
                                                                                                                                                                                                     gTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAAt                        <  1:99725/249‑1 (MQ=255)
                                                                                                                                                                                                         ggCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGtt                   >  1:47189/1‑251 (MQ=255)
                                                                                                                                                                                                           cTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGtttt                 >  1:481573/1‑251 (MQ=255)
                                                                                                                                                                                                                 tGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATcc           >  2:129227/1‑251 (MQ=255)
                                                                                                                                                                                                                    cGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCCCTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGCTCCGGCTTTCATCCGTTCACTGATCGAGGATCGGCTGGATGCgatgttaatcccgc        >  2:365576/1‑238 (MQ=255)
                                                                                                                                                                                                                    cGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGcc                                                                                             >  2:310494/1‑166 (MQ=255)
                                                                                                                                                                                                                    cGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGcc                                                                                             <  1:310494/166‑1 (MQ=255)
                                                                                                                                                                                                                         gTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTc   <  1:241790/251‑1 (MQ=255)
                                                                                                                                                                                                                          tATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCCGCAGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGCTCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTCATCAAGGATCGCCGGCATGCGTTTTTCATCCggcggagcg  >  2:68209/1‑244 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              
CTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCG  >  NC_000913/4624127‑4624592

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: