Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,167,940 C→T G114G (GGC→GGT ycfJ → PF05433 family protein YcfJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,167,9400CT98.1% 155.7 / NA 53G114G (GGC→GGTycfJPF05433 family protein YcfJ
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (27/25);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

ATACGCTCAGGTTGTTTCTGCAACCCCAATCAAGGAAACGGTTAAAACACCGCGTCAGGAGTGTCGCAACGTCACAGTGACCCATCGTCGACCGGTGCAGGATGAAAATCGCATTACCGGGTCGGTGCTCGGCGCTGTTGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGCTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTG  >  NC_000913/1167697‑1168153
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                     
aTACGCTCAGGTTGTTTCTGCAACCCCAATCAAGGAAACGGTTAAAACACCGCGTCAGGAGTGTCGCAACGTCACAGTGACCCATCGTCGACCGGTGCAGGATGAAAATCGCATTACCGGGTCGGTGCTCGGCGCTGTTGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTcc                                                                                                                                                                                                                >  1:15218/1‑251 (MQ=255)
       cAGGTTGTTTCTGCAACCCCAATCAAGGAAACGGTTAAAACACCGCGTCAGGAGTGTCGCAACGTCACAGTGACCCATCGTCGACCGGTGCAGGATGAAAATCGCATTACCGGGTCGGTGCTCGGCGCTGTTGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGaa                                                                                                                                                                                                           >  1:154600/1‑249 (MQ=255)
         ggTTGTTTCTGCAACCCCAATCAAGGAAACGGTTAAAACACCGCGTCAGGAGTGTCGCAACGTCACAGTGACCCATCGTCGACCGGTGCAGGATGAAAATCGCATTACCGGGTCGGTGCTCGGCGCTGTTGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCg                                                                                                                                                                                                       >  2:241511/1‑251 (MQ=255)
                     aaCCCCAATCAAGGAAACGGTTAAAACACCGCGTCAGGAGTGTCGCAACGTCACAGTGACCCATCGTCGACCGGTGCAGGATGAAAATCGCATTACCGGGTCGGTGCTCGGCGCTGTTGCTGGCGGCGTGCCAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCAACTGTTGTGGGGGCGCTGGGTGGTGGATATGCAGGTAACCAAATCCAGGGTTCTCTCCAGGAAAGGGATACTTaca                                                                                                                                                                                              >  2:387998/1‑248 (MQ=255)
                     aaCCCCAAACAAGGAAACGGTTAAAACACCGCGTCAGGAGTGTCGCAACGTCACAGTGACCCATCGTCGACCGGTGCAGGATGAAAATCGCATTACCGGGTCGGTGCTCGGCGCGGTTGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGa                                                                                                                                                                                           <  2:154600/251‑1 (MQ=255)
                       ccccAATCAAGGAAACGGTTAAAACACCGCGTCAGGCGTGTCGCAACGTCACAGTGACCCATCGTCGACCGGTGCAGGATGAAAATCGCATTACCGGGTCGGTGCTCGGCGCTGTTGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGCGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGGTATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACt                                                                                                                                                                                         >  2:216995/1‑251 (MQ=255)
                          cAATCAAGGAAACGGTTAAAACACCGCGTCAGGAGTGTCGCAACGTCACAGTGACCCATCGTCGACCGGTGCAGGATGAAAATCGCATTACCGGGTCGGTGCTCGGCGCTGTTGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACg                                                                                                                                                                                      <  1:575770/251‑1 (MQ=255)
                            aTCAAGGAAACGGTTAAAACACCGCGTCAGGAGTGTCGCAACGTCACAGTGACCCATCGTCGACCGGTGCAGGATGAAAATCGCATTACCGGGTCGGTGCTCGGCGCTGTTGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTTGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCa                                                                                                                                                                                    <  2:278509/251‑1 (MQ=255)
                                                   gcgTCAGGAGTGTCGCAACGTCACAGTGACCCATCGTCGACCGGTGCAGGATGAAAATCGCATTACCGGGTCGGTGCTCGGCGCTGTTGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATg                                                                                                                                                             <  1:439926/251‑1 (MQ=255)
                                                    cgTCAGGAGTGTCGCAACGTCACAGTGACCCATCGTCGACCGGTGCAGGATGAAAATCGCATTACCGGGTCGGTGCTCGGCGCTGTTGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGa                                                                                                                                                            <  1:147395/251‑1 (MQ=255)
                                                                  cAACGTCACAGTCACTCATCGTCGACCGATGCAGGATGAAAATCGCATTACCGGGTCGGAGCTCGGCGCTGTTGCTGTCGGCGTGATAGTGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGGGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGGTCAAGGGTTCTCTCCAGGAAAGGGATACTTACACGACTACGCAACAGCGTTGTAAAACGGAGTATGACAAGGCAGaaaaa                                                                                                                                               >  2:379991/1‑250 (MQ=255)
                                                                     cGTCACAGTGACCCATCGTCGACCGGTGCAGGATGAAAATCGCATTACCGGGTCGGTGCTCGGCGCTGTTGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGc                                                                                                                                           <  1:323668/251‑1 (MQ=255)
                                                                         acaGTGACCCATCGTCGACCGGTGCAGGATGAAAATCGCATTACCGGGTCGGTGCTCGGCGCTGTTGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCgg                                                                                                                                       >  1:296774/1‑251 (MQ=255)
                                                                              taCCCATCGTCGACCGGTGCAGGATGAAAATCGCATTACCGGGTCGGTGCTCGGCGCTGTTGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTat                                                                                                                                   <  2:6742/249‑1 (MQ=255)
                                                                                  cATCGTCGACCGGTGCAGGATGAAAATCGCATTACCGGGTCGGTGCTCGGCGCTGTTGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGAtgt                                                                                                                              >  2:568588/1‑251 (MQ=255)
                                                                                        cgACCGGTGCAGGGTGAAAATCGCATTACCGGGTCGGTGCTCGGCGCTGTTGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCt                                                                                                                         >  2:656723/1‑250 (MQ=255)
                                                                                           ccGGTGCAGGATGAAAATCGCATTACCGGGTCGGTGCTCGGCGCTGTTGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCtata                                                                                                                      >  2:581347/1‑250 (MQ=255)
                                                                                              gTGCAGTAAGAAAATCGCATTACCGGGTCGGTGCTCAGCACTGTCTCTGGCCGAGTGATAGCGCATCAGTTTGGGGGAGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGGGTTGGATATGCCGGTAACCCAATCCAGGGTTCTCTCCAGGACAGCGATACTTACACGCCTACTCCACAGCGTTGTAACACGGTTGATACCAAGTAAGCAAAAATCCTCGGGTATGATGTGACCtatc                                                                                                                      >  2:318205/1‑246 (MQ=255)
                                                                                                gCAGGATGAAAATCGCATTACCGGGTCGGTGCTCGGCGCTGTTGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTg                                                                                                                >  1:644323/1‑251 (MQ=255)
                                                                                                  aGGATGAAAATCGCATTACCGGGTCGGTGCTCGGCGCTGTTGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGc                                                                                                              >  1:424861/1‑251 (MQ=255)
                                                                                                   ggATGAAAATCGCATTACCGGGTCGGTGCTCGGCGCTGTTGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCg                                                                                                             <  2:450048/251‑1 (MQ=255)
                                                                                                     aTGAAAATCGCATTACCGGGTCGGTGCTCGGCGCTGTTGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCAACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAAAGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGAt                                                                                                           >  1:2109/1‑251 (MQ=255)
                                                                                                                aTTACCGGGTCGGTGCTCGGCGCTGTTGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCa                                                                                                 >  1:6620/1‑250 (MQ=255)
                                                                                                                    ccGGGTCGGTGCTCTGCGCTGTTGCTGGCGGCGTGATAGGGCATCAGATTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGAGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTACAGGAAAGCGAGACTTACACGACTACGAAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATAAGCAGGGCAAAATc                                                                                            >  2:619721/1‑251 (MQ=255)
                                                                                                                     cGGGTCGGTGCTCGGCGCTGTTGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATcc                                                                                           <  1:241511/251‑1 (MQ=255)
                                                                                                                      gggacggggCTCGGCGCTGTTGCTGGCGGCGTGATAGGGCATCAGGTTGGTGGTGGTCGCGGTAAAGATGGCGCCACTGTTGTGGGGGAGCTGGGTGGTGGGTATGCCGGTAACCAGATCCAGGGTGCTCTACAGGAAAGCGATACTGACACGACTACGCAAAAGCGTTGTAAAACGGTGTATGACAAGTCAGAAACAATGCTCGGTTATGATGTAAACTATAAGCATGGCGATCAGCAGGACAAAATCCg                                                                                          >  2:430383/9‑251 (MQ=255)
                                                                                                                        gTCGGTGCTCGGCGCTGTTGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCa                                                                                        <  1:603447/251‑1 (MQ=255)
                                                                                                                         gcGGTGGTCGGCGCTGTTGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCAt                                                                                       <  2:149800/250‑1 (MQ=255)
                                                                                                                                tCGGCGCTGTTGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACcgc                                                                                >  1:477722/1‑251 (MQ=255)
                                                                                                                                         ttGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCAt                                                                                       <  2:114001/235‑1 (MQ=255)
                                                                                                                                         ttGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCAt                                                                                       >  1:114001/1‑235 (MQ=255)
                                                                                                                                           gCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGGAACCAGATCCAGGGTTCTCTCCAGGAAAGAGATACTTACACGACTACGCAACAGCGTTGTAAAACGGCGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTAc                                                                     >  1:222445/1‑251 (MQ=255)
                                                                                                                                             tGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGc                                                                   >  2:212063/1‑251 (MQ=255)
                                                                                                                                             gggcggcGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGc                                                                   <  1:568588/250‑1 (MQ=255)
                                                                                                                                                  gcgTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGAGCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATc                                                              <  1:281537/251‑1 (MQ=255)
                                                                                                                                                       aTAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGGTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTGCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAAGCCGCATGGACCGCGATCCGGGTACGCAGAGCCCGCt                                                         <  2:457969/251‑1 (MQ=255)
                                                                                                                                                            gCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTTCACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATa                                                    <  1:28041/251‑1 (MQ=255)
                                                                                                                                                            gCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATa                                                    <  2:170465/251‑1 (MQ=255)
                                                                                                                                                             cATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGAt                                                     >  2:307466/1‑249 (MQ=255)
                                                                                                                                                             cATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGAt                                                     >  1:411032/1‑249 (MQ=255)
                                                                                                                                                                cAGTTTGGTGGTGGTCGCGGTAAATATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCa                                                                                                                                                                                                               <  2:654461/92‑1 (MQ=255)
                                                                                                                                                                      ggtggtggtCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATAGGTCGGTAACCAGAGCGCGGGTTCGCTCCATGAAAGGGTTACTTACACGACTGAGCAACAGGGTTTCAAAACGGTGTATGAGAAGTCAGAAAAGATGTTCGGTTATGATTGGACCTAGAAGAGTGGCGATCAGTAGGGTAAAGTACGCATGGTCCGGGGTCCGGGGACggcgagcccggtggatagtaatgggc                                           <  2:318063/250‑26 (MQ=255)
                                                                                                                                                                         ggtggtCGCGATCAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTTGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGACAGCGATACTTACCCGACTACGCAACAGCGTTGTAAAACGGTGTCTGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACt                                       >  2:402535/1‑251 (MQ=255)
                                                                                                                                                                           tggtCGGGGTAAAGATGTCGCCACTGTTGTGGGGGAGCTGGGTGGAGGATATGCAGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGGGATACTTACACAACTACGCCACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGAGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGAATGGACCGCGATCCGGGTACGCAGATCCAGCTAGATAGCCACGGGCAANTg                                      >  2:588690/1‑250 (MQ=255)
                                                                                                                                                                               cgcgGTAAAGATGTCGCCACTGTTGTGGGGGCGCTTGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCATGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACTGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCACCGGGCAACTGAttt                                  <  2:399682/250‑1 (MQ=255)
                                                                                                                                                                                       aGATGTCGCCACTGTTGTGGGGGCGCTGTGTGGTGGATATGCCTGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACa                         <  1:379991/251‑1 (MQ=255)
                                                                                                                                                                                              gCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGtata                  <  1:30665/251‑1 (MQ=255)
                                                                                                                                                                                                cACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGCTCTCTCCAGGAAAGCGATACTTACACGACTAAGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCAGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGAATAAc                >  1:75308/1‑251 (MQ=255)
                                                                                                                                                                                                   tgttgtGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACa               >  2:448838/1‑249 (MQ=255)
                                                                                                                                                                                                    gttgtGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAg             <  1:162267/250‑1 (MQ=255)
                                                                                                                                                                                                        tgGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGt        >  1:308251/1‑251 (MQ=255)
                                                                                                                                                                                                         gggggAGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGAGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGATAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTa       <  2:244752/251‑1 (MQ=255)
                                                                                                                                                                                                              cgcTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGTTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACTCAGATCCCTCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTg  <  1:142492/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                     
ATACGCTCAGGTTGTTTCTGCAACCCCAATCAAGGAAACGGTTAAAACACCGCGTCAGGAGTGTCGCAACGTCACAGTGACCCATCGTCGACCGGTGCAGGATGAAAATCGCATTACCGGGTCGGTGCTCGGCGCTGTTGCTGGCGGCGTGATAGGGCATCAGTTTGGTGGTGGTCGCGGTAAAGATGTCGCCACTGTTGTGGGGGCGCTGGGTGGTGGATATGCCGGTAACCAGATCCAGGGCTCTCTCCAGGAAAGCGATACTTACACGACTACGCAACAGCGTTGTAAAACGGTGTATGACAAGTCAGAAAAAATGCTCGGTTATGATGTGACCTATAAGATTGGCGATCAGCAGGGCAAAATCCGCATGGACCGCGATCCGGGTACGCAGATCCCGCTAGATAGCAACGGGCAACTGATTTTGAATAACAAAGTATAACAAGGCTGTACTCTG  >  NC_000913/1167697‑1168153

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