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

GACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTGATAGCAGTAAACAGCGCCACCACAAACGGCAGCTGGGTCGTGGCTTCACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCCGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTG  >  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)
                                            aaCAGCCCCACCACAAACGGCAGCGGGGTCGGGGCTACACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACGGCTCGCAAACCACAGCGTTCCACACCACTTTTCTGTCATCAATCTGCAACTCCTGATGCGGACACTGTGGCACGCAC‑Tcccc                                                                                                                                                                                                     <  1:179522/251‑1 (MQ=255)
                                                          aagcgggtgctgggcggtgctttcgccGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTGGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTGCCACACCACTTTGCTGTCAACAAGCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TACCCACAGTCATTGCAAc                                                                                                                                                                                       <  2:249752/226‑1 (MQ=255)
                                                           aaCGGCAGCTGGGTCGTGGCTTCACCGACACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCCAACCACAGCGTTCCACACCCCTTTGCTGTCAACCATCTGAAAAGCCGGATGCGGAAACTGTGGCAACCAC‑TCCCCACCGTCCTCGCAACg                                                                                                                                                                                      >  2:192774/1‑251 (MQ=255)
                                                               gCAGCTGGGGCGTGGCTTCACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTccc                                                                                                                                                                                  <  2:106029/251‑1 (MQ=255)
                                                                 aGCTGGGTCGTGGCTTCACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGCCATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCCCAGTCATTGCAACGTCCCAt                                                                                                                                                                                >  2:111862/1‑251 (MQ=255)
                                                                 aGCTGGGTCGTGGCTTCACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCAt                                                                                                                                                                                >  1:227113/1‑251 (MQ=255)
                                                                      ggTCGTGGCTTCACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCATCGTcc                                                                                                                                                                           <  2:10070/251‑1 (MQ=255)
                                                                         tgTGGCCTCACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACATCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCAGCGTCCAc                                                                                                                                                                         <  2:160652/249‑1 (MQ=255)
                                                                         gggggCTCCACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCATCGTCCACg                                                                                                                                                                        <  2:13485/248‑1 (MQ=255)
                                                                           tGGCTTCACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCATCGTCCACgg                                                                                                                                                                       >  1:39472/1‑250 (MQ=255)
                                                                            ggCTTCACCGCCACTCACCGTTATCCCTTCGATAACCAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTAATGGATTGAGCCCTGGGCGTGGCATGTTGCGGACCCCTCTTCCGGCAGGTATAACCCTGCTCGCAACCCACAGCGTTCCAAACCACTTTGCTGTCAACAATCTGCAACGCCTGATCCGGACACTGTGGAACGCAC‑TCCCCACAGTCATTGCAACGTCCCATAGTCCACGGAt                                                                                                                                                                     >  1:253976/1‑251 (MQ=255)
                                                                               ttCACCGCCACTCACCGTTATCCCTTCGATAAACAAAACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGAGGACACCTCTTCAGGCAGGCATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACACCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTACCATCGTCCACGGATtgt                                                                                                                                                                  >  1:275538/1‑251 (MQ=255)
                                                                                            aCCGATATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCATCGTCCACGAATTGTGACAGATTTTGc                                                                                                                                                      >  2:202633/1‑250 (MQ=255)
                                                                                                         cGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGc                                                                                                                                          >  1:246459/1‑249 (MQ=255)
                                                                                                          gATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCCACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGc                                                                                                                                       >  2:187863/1‑251 (MQ=255)
                                                                                                                                    ggCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCCCGCAC‑TCCCCACAGTCATTGCAACGTCCCATCGTCAACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTg                                                                                                             >  1:246394/1‑251 (MQ=255)
                                                                                                                                        aaGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACAACTCTTCAGGCAGGTATCACACTGCTCGCAAACCCCAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTg                                                                                                         >  2:259807/1‑251 (MQ=255)
                                                                                                                                          gCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGAGGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCATCGTCAACGGCTTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCCAAAACAGAGCCAGACGACGGccggg                                                                                                        >  2:209583/1‑247 (MQ=255)
                                                                                                                                            aCTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACAACACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGccc                                                                                                     >  1:40330/1‑251 (MQ=255)
                                                                                                                                                       gCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCa                                                                                                                                 >  1:269454/1‑212 (MQ=255)
                                                                                                                                                       gCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCAGCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCa                                                                                                                                 <  2:269454/212‑1 (MQ=255)
                                                                                                                                                          cATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑GCCCCACAGTCATTGCAACGTCCCAGCTTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCGGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGgaga                                                                                       <  1:172552/251‑1 (MQ=255)
                                                                                                                                                             ggATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAg                                                                                     <  2:204298/250‑1 (MQ=255)
                                                                                                                                                              gATTGCGCCATGGGCGTGGCATGTTGCTGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTGCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATGGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGTCAgac                                                                                                                  <  2:223624/220‑1 (MQ=255)
                                                                                                                                                              gATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAAGCAGGTATCACACTGCTAGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAgac                                                                                                                  >  1:223624/1‑220 (MQ=255)
                                                                                                                                                                  gcgcCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACGGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATa                                                                               <  2:5538/251‑1 (MQ=255)
                                                                                                                                                                      cATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCt                                                                           >  1:244156/1‑251 (MQ=255)
                                                                                                                                                                       aTGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCtt                                                                          <  1:38216/251‑1 (MQ=255)
                                                                                                                                                                              tGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACt                                                                   >  1:113032/1‑251 (MQ=255)
                                                                                                                                                                                       gCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCa                                                           <  2:227113/250‑1 (MQ=255)
                                                                                                                                                                                                  ttCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCAt                                               >  2:112193/1‑251 (MQ=255)
                                                                                                                                                                                                   tcatgcaggTATCACACTGCTCGCAAACCACAGCGTTCCACACCTCTTTGCTGTCAACCATCTGCAACGCCTGAGGCGTACACTGTGGCAC‑CACCGCCCCACAGTCATTGCATCTTCCCAACGTCCACGGATTGTGACAGTTTGTGCAGCTCAGATTGTAGCCCTGCAAGAACACAGCCAGACGACTGCCTGGCCCGTCAACCCAGAAGAAGGGTATAAACTTACTGACGATATCGCATGTGCTGtaca                                                <  2:171508/245‑4 (MQ=255)
                                                                                                                                                                                                           gTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAg                                                                                     >  2:183655/1‑204 (MQ=255)
                                                                                                                                                                                                           gTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAg                                                                                     <  1:183655/204‑1 (MQ=255)
                                                                                                                                                                                                                           aaaCCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTc                      >  2:179487/1‑251 (MQ=255)
                                                                                                                                                                                                                             aCCACAGCGTTCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTGTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCATACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAAACTTACTGACTAAATCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCTTTCCa                    <  1:138340/251‑1 (MQ=255)
                                                                                                                                                                                                                              ccACAGCGTTCCCCACCCCTTTGCTGTCAAAAAGCTGCAACGCCTTATGCGGCCACTGTGGCACGCCC‑TCCCCACTGTCATTTCAACGACCCATCGTCCACGGATTGTGACAGTTTTTGCTGCGCAGATTGCAGCCCGGCAAAAACAGAGCCATACGACTGCCTGGTCCGGCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTAGCACGCGCGGCTGGCGTTCCAg                   <  2:52904/251‑1 (MQ=255)
                                                                                                                                                                                                                                      ttCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATTCGGACACTGTTGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGGTTTGGCAGCGCAGATGGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGTATAATCGTACTGACTAAAACGCAGCTGCTGTTCATGGCTTAGCGCGCTCTGCTGGCGGACCAGAATACg             <  2:253976/249‑1 (MQ=255)
                                                                                                                                                                                                                                       tCCACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGACATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGCCGACTGACTGGCCCGGCAAAGCAGGAGCAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTGCATGGCTTATCACGCGCGGCGGGCGGTCCAGAATAagagg          >  2:123803/1‑250 (MQ=255)
                                                                                                                                                                                                                                        ccACACCACTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCCGGCGTTCCAGAATACGAgtg         >  1:115151/1‑251 (MQ=255)
                                                                                                                                                                                                                                           caccacCTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTtg      <  2:66581/251‑1 (MQ=255)
                                                                                                                                                                                                                                               acTTTGCTGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTtgcgtg  >  2:202166/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
GACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTGATAGCAGTAAACAGCGCCACCACAAACGGCAGCTGGGTCGTGGCTTCACCGCCACTCACCGTTATCCCTTCGATAAACAGCACTGCTTTGCGGACATGGCTAAGCACTTCGTCCACGCTCATGGATTGCGCCATGGGCGTGGCATGTTGCGGACACCTCTTCAGGCAGGTATCACACTGCTCGCAAACCACAGCGTTCCACACCACTTTGCCGTCAACAATCTGCAACGCCTGATGCGGACACTGTGGCACGCAC‑TCCCCACAGTCATTGCAACGTCCCATCGTCCACGGATTGTGACAGTTTTTGCAGCGCAGATTGCAGCCCTGCAAAAACAGAGCCAGACGACTGCCTGGCCCGTCAACGCAGGAGAAGGGGATAATCTTACTGACTAAAGCGCATCTGCTGTTCATGGCTTATCACGCGCGGCTGGCGTTCCAGAATACGAGTGTTGCGTG  >  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: ‑