Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,615,345 C→T G67S (GGC→AGC)  yjjW ← putative glycyl‑radical enzyme activating enzyme YjjW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,615,3450CT98.2% 187.9 / ‑6.1 55G67S (GGC→AGC) yjjWputative glycyl‑radical enzyme activating enzyme YjjW
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (29/25);  total (29/26)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTGATAGCAGTAAACAGCGCCACCACAAACGGCAGCTGGGTCGTGGCTTCACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCCGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTG  >  NC_000913/4615099‑4615588
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
gaaaGGTGAGAGGGCGCAGTTGGGGATCGTTGTTGATAGCAGTAAACAGCGCCACCACAAACGGCAGCGGGGTCGTGGCTTCACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTca                                                                                                                                                                                                                                                  <  1:41778/248‑1 (MQ=255)
   aGGTGAGATGGCGCAGTTGCGGATCGTTTTTGATAGCAGTAAACAGCGCCACCACAAACGGCAGCTGGGTCGTGGCTTCACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTcaac                                                                                                                                                                                                                                                >  2:112302/1‑250 (MQ=255)
    ggTGAGATGGCGCAGTTGCGGATCTTTTTGGATAGCAGTAAACAGCGCCACCACAAACGGCAGCTGGGTCGTGGCTTCACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTcaacaa                                                                                                                                                                                                                                              <  1:49396/251‑1 (MQ=255)
        agatgggtcAGTTGTGGTTCTTTTTTGATAGTAGGAACCACCGCCACCACAAACGGCAGCGGTGTCGTGCCTCCACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACGGCTCGCAAACCCCAGCGTTCCACACCACTTTGCTGTCAACAATCTg                                                                                                                                                                                                                                          <  1:129273/243‑1 (MQ=255)
        agaTGGCGCAGTTGCGGATCGTTTTTGATAGCAGTAAACAGCGCCACCACAAACGGCAGCTGGGTCGTGGCTTCACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTg                                                                                                                                                                                                                                          >  1:204298/1‑251 (MQ=255)
           tGGCGCAGTTGCGGATCGTTTTTGATAGCAGTAAACAGCCCCACCACAAACGGCAGCTGGGTCGGGGCTTCACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCa                                                                                                                                                                                                                                        <  1:196078/250‑1 (MQ=255)
                    tGCGGATCGTTTTTGATAGCAGTAAACAGCGCCACCACAAACGGCAGCTGGGTCGTGGCTTCACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATg                                                                                                                                                                                                                              >  1:148964/1‑251 (MQ=255)
                        gATCGTTTTTGATAGCAGTAAACAGCGCCACCACAAACGGCAGCTGGGTCGTGGCTTCACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCgg                                                                                                                                                                                                                           >  1:186994/1‑250 (MQ=255)
                              ttttgagatctgTAAACACCGCCCCCACAAACTGCCGCTGGTCCGTGGGTTCACCGCCACTCACCGTTATCCCTTCGATAAACAGCATTGCGTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGAGACCTCTTCAGGCAGGTATCGCACTGCTCGCAAACCACAGCGTTGCACGCCACGTTGCCGTCAACAATCTGCAACGCCTGATGCGGACACtg                                                                                                                                                                                                                     <  1:111491/239‑1 (MQ=255)
                                 ggATAGCAGAAAACAGCGCCACCAAAAAGGGCAGCGGGGTGGGGGCTTAACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGgc                                                                                                                                                                                                                 <  1:116444/250‑1 (MQ=255)
                                  gATAGCAGTAAACAGCGCCACCACAAACGGCAGCTGGGTCGTGGCTTCACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGgc                                                                                                                                                                                                                 >  1:7200/1‑250 (MQ=255)
                                   aTAGCAGTAAACAGCGCCACCACAAACGGCAGCTGGGTCGTGGCTTCACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCCCACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGgcac                                                                                                                                                                                                               >  2:219782/1‑251 (MQ=255)
                                    acgCAGGAAACGGCGCCCCCACAACCGGCTGGGGGGTCTTGTCTGCACCGCCACTCACCGTTATCCCGTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGgcacg                                                                                                                                                                                                              <  2:80329/249‑1 (MQ=255)
                                            aaCAGCCCCACCACAAACGGCAGCGGGGTCGGGGCTACACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACGGCTCGCAAACCACAGCGTTCCACACCACTTTTCTGTCATCAATCTGCAACTCCTGATGCGGACACTGTGGCACGCACTcccc                                                                                                                                                                                                     <  1:179522/251‑1 (MQ=255)
                                                          aagcgggtgctgggcggtgctttcgccGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTGGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTGCCACACCACTTTGCTGTCAACAAGCTGCAACGCCTGATGCGGACACTGTGGCACGCACTACCCACAGTCATTGCAAc                                                                                                                                                                                       <  2:249752/226‑1 (MQ=255)
                                                           aaCGGCAGCTGGGTCGTGGCTTCACCGACACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCCAACCACAGCGTTCCACACCCCTTTGCTGTCAACCATCTGAAAAGCCGGATGCGGAAACTGTGGCAACCACTCCCCACCGTCCTCGCAACg                                                                                                                                                                                      >  2:192774/1‑251 (MQ=255)
                                                               gCAGCTGGGGCGTGGCTTCACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTccc                                                                                                                                                                                  <  2:106029/251‑1 (MQ=255)
                                                                 aGCTGGGTCGTGGCTTCACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGCCATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCCCAGTCATTGCAACGTCCCAt                                                                                                                                                                                >  2:111862/1‑251 (MQ=255)
                                                                 aGCTGGGTCGTGGCTTCACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCAt                                                                                                                                                                                >  1:227113/1‑251 (MQ=255)
                                                                      ggTCGTGGCTTCACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTcc                                                                                                                                                                           <  2:10070/251‑1 (MQ=255)
                                                                         tgTGGCCTCACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACATCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCAGCGTCCAc                                                                                                                                                                         <  2:160652/249‑1 (MQ=255)
                                                                         gggggCTCCACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACg                                                                                                                                                                        <  2:13485/248‑1 (MQ=255)
                                                                           tGGCTTCACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACgg                                                                                                                                                                       >  1:39472/1‑250 (MQ=255)
                                                                            ggCTTCACCGCCACTCACCGTTATCCCTTCGATAACCAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTAATGGATTGAGCCCTGGGCGTGGCATGTTGCGGACCCCTCTTCCGGCAGGTATAACCCTGCTCGCAACCCACAGCGTTCCAAACCACTTTGCTGTCAACAATCTGCAACGCCTGATCCGGACACTGTGGAACGCACTCCCCACAGTCATTGCAACGTCCCATAGTCCACGGAt                                                                                                                                                                     >  1:253976/1‑251 (MQ=255)
                                                                               ttCACCGCCACTCACCGTTATCCCTTCGATAAACAAAACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGAGGACACCTCTTCAGGCAGGCATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACACCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTACCATCGTCCACGGATtgt                                                                                                                                                                  >  1:275538/1‑251 (MQ=255)
                                                                                            aCCGATATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGAATTGTGACAGATTTTGc                                                                                                                                                      >  2:202633/1‑250 (MQ=255)
                                                                                                         cGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGc                                                                                                                                          >  1:246459/1‑249 (MQ=255)
                                                                                                          gATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCCACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGc                                                                                                                                       >  2:187863/1‑251 (MQ=255)
                                                                                                                                    ggCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCCCGCACTCCCCACAGTCATTGCAACGTCCCATCGTCAACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTg                                                                                                             >  1:246394/1‑251 (MQ=255)
                                                                                                                                        aaGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACAACTCTTCAGGCAGGTATCACACTGCTCGCAAACCCCAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTg                                                                                                         >  2:259807/1‑251 (MQ=255)
                                                                                                                                          gCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGAGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCAACGGCTTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCCAAAACAGAGCCAGACGACGGccggg                                                                                                        >  2:209583/1‑247 (MQ=255)
                                                                                                                                            aCTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACAACACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGccc                                                                                                     >  1:40330/1‑251 (MQ=255)
                                                                                                                                                       gCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCa                                                                                                                                 >  1:269454/1‑212 (MQ=255)
                                                                                                                                                       gCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCAGCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCa                                                                                                                                 <  2:269454/212‑1 (MQ=255)
                                                                                                                                                          cATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACGCCCCACAGTCATTGCAACGTCCCAGCTTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCGGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGgaga                                                                                       <  1:172552/251‑1 (MQ=255)
                                                                                                                                                             ggATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAg                                                                                     <  2:204298/250‑1 (MQ=255)
                                                                                                                                                              gATTGCGCCATGGGCGTGGCATGTTGCTGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTGCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATGGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGTCAgac                                                                                                                  <  2:223624/220‑1 (MQ=255)
                                                                                                                                                              gATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAAGCAGGTATCACACTGCTAGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAgac                                                                                                                  >  1:223624/1‑220 (MQ=255)
                                                                                                                                                                  gcgcCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACGGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATa                                                                               <  2:5538/251‑1 (MQ=255)
                                                                                                                                                                      cATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCt                                                                           >  1:244156/1‑251 (MQ=255)
                                                                                                                                                                       aTGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCtt                                                                          <  1:38216/251‑1 (MQ=255)
                                                                                                                                                                              tGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACt                                                                   >  1:113032/1‑251 (MQ=255)
                                                                                                                                                                                       gCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCa                                                           <  2:227113/250‑1 (MQ=255)
                                                                                                                                                                                                  ttCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCAt                                               >  2:112193/1‑251 (MQ=255)
                                                                                                                                                                                                   tcatgcaggTATCACACTGCTCGCAAACCACAGCGTTCCACACCTCTTTGCTGTCAACCATCTGCAACGCCTGAGGCGTACACTGTGGCACCACCGCCCCACAGTCATTGCATCTTCCCAACGTCCACGGATTGTGACAGTTTGTGCAGCTCAGATTGTAGCCCTGCAAGAACACAGCCAGACGACTGCCTGGCCCGTCAACCCAGAAGAAGGGTATAAACTTACTGACGATATCGCATGTGCTGtaca                                                <  2:171508/245‑4 (MQ=255)
                                                                                                                                                                                                           gTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAg                                                                                     >  2:183655/1‑204 (MQ=255)
                                                                                                                                                                                                           gTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAg                                                                                     <  1:183655/204‑1 (MQ=255)
                                                                                                                                                                                                                           aaaCCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTc                      >  2:179487/1‑251 (MQ=255)
                                                                                                                                                                                                                             aCCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTGTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCATACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAAACTTACTGACTAAATCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCTTTCCa                    <  1:138340/251‑1 (MQ=255)
                                                                                                                                                                                                                              ccACAGCGTTCCCCACCCCTTTGCTGTCAAAAAGCTGCAACGCCTTATGCGGCCACTGTGGCACGCCCTCCCCACTGTCATTTCAACGACCCATCGTCCACGGATTGTGACAGTTTTTGCTGCGCAGATTGCAGCCCGGCAAAAACAGAGCCATACGACTGCCTGGTCCGGCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTAGCACGCGCGGCTGGCGTTCCAg                   <  2:52904/251‑1 (MQ=255)
                                                                                                                                                                                                                                      ttCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATTCGGACACTGTTGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGGTTTGGCAGCGCAGATGGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGTATAATCGTACTGACTAAAACGCAGCTGCTGTTCATGGCTTAGCGCGCTCTGCTGGCGGACCAGAATACg             <  2:253976/249‑1 (MQ=255)
                                                                                                                                                                                                                                       tCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGACATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGCCGACTGACTGGCCCGGCAAAGCAGGAGCAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTGCATGGCTTATCACGCGCGGCGGGCGGTCCAGAATAagagg          >  2:123803/1‑250 (MQ=255)
                                                                                                                                                                                                                                        ccACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCCGGCGTTCCAGAATACGAgtg         >  1:115151/1‑251 (MQ=255)
                                                                                                                                                                                                                                           caccacCTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTtg      <  2:66581/251‑1 (MQ=255)
                                                                                                                                                                                                                                               acTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTtgcgtg  >  2:202166/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
GACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTGATAGCAGTAAACAGCGCCACCACAAACGGCAGCTGGGTCGTGGCTTCACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCCGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCACTCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTG  >  NC_000913/4615099‑4615588

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