Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,625,170 C→T H92H (CAC→CAT)  serB → phosphoserine phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,625,1700CT97.2% 189.4 / ‑2.3 71H92H (CAC→CAT) serBphosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (43/26);  total (45/26)
Fisher's exact test for biased strand distribution p-value = 5.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01

GTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑ACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACC‑GCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCA  >  NC_000913/4624931‑4625415
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       
gTCTCTTTATGGCCGGGTTTTCCTCCTGTATTAAGTGGTGAGGAAGGGATGACACGGGATTACCACGCAGGTCGTAGCGGCGGGCTGCGGTATGTGCGTGGGCGGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTTGTACGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCTCTCCCCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCaact                                                                                                                                                                                                                                              <  2:368476/251‑4 (MQ=255)
    ctTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGAt                                                                                                                                                                                                                                          >  1:201023/1‑251 (MQ=255)
        aTGGCCGGGTCTTCCTCTTTCATTAAGTGGTGATGGATTGATTCCACTGGATTACCACGCAGGTCGTTGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGGCCCAATACCAGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATggcg                                                                                                                                                                                                                                      <  2:180288/251‑4 (MQ=255)
        aTGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGTGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTc                                                                                                                                                                                                                                       <  1:665292/250‑1 (MQ=255)
               ggTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCt                                                                                                                                                                                                                               >  1:640237/1‑251 (MQ=255)
                   tGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGGCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCTCTCACCGCACGGGATCC‑GCGCCTGGCCCATGAACCGCAACTGGATGTCGCCCCGCTgggg                                                                                                                                                                                                                           >  2:116378/1‑251 (MQ=255)
                    gcatctTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGAAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGACTGACCCAATACCAGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCCCCTGGATCGTGGAAAATT‑CAAGGAGGATTCGTCTGGCAAGTTCACGCACCGCACGGGCTACGAC‑CCTGGACAATGAAGCGCACATGGATGTCGCCCCGCGgggg                                                                                                                                                                                                                           >  2:365558/4‑250 (MQ=255)
                    gCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTgggg                                                                                                                                                                                                                           >  1:541493/1‑250 (MQ=255)
                                  ttggagatgaAGTGATGCCACTGGATTAACACGCAGGTCGTAGCGGCTGGCTGCTGAATGGTCGTGGGCTGTATAAACAGCGTCTGAACCAATACCAGAGCAACCTGAGTGCGGCAGTG‑GGAATTGTTGCCGCCTGGTGCGTGGAAGATTATCATG‑AGAGTAGTCTGGCAGGTTCACTCACCGCACGGGATAC‑GCGCCCAACCCATGAAGCGCCACTGGATGGCGTCCCGCTggggtaacacccgaccctg                                                                                                                                                                                                            >  2:395847/6‑236 (MQ=255)
                                  gtggtgATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTg                                                                                                                                                                                                           >  1:91454/1‑251 (MQ=255)
                                   tggtgATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATACCGCACCTgc                                                                                                                                                                                                          >  1:300967/1‑251 (MQ=255)
                                         tgaAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGACGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑AGATTCGTCTGGCAGATTCACTCACCGCCCGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCGGCACCTGCGCaagc                                                                                                                                                                                                    >  2:572740/1‑248 (MQ=255)
                                           aaGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATG‑GAGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTTACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAAATGGATGTCGCCCCGCTGGGG‑AAACACCCGCCACTGCGCACGCCg                                                                                                                                                                                                  >  2:93958/1‑251 (MQ=255)
                                            agggtggccccgggagtACCACGCAGGCCGTATCGGCGTGCTGCGGTATGGTCGTTGGCTTGAGAAACAGCGTCTGACCCAATCCCAGAGCAAACTGGGTGCGGCGAGG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCCCCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCCGAGCCCGCCgg                                                                                                                                                                                                 <  2:375592/235‑1 (MQ=255)
                                               gATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGACGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGtt                                                                                                                                                                                              >  1:575090/1‑251 (MQ=255)
                                                aTGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTTACCCAATACCAGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTAGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCTCGCTGGGG‑AAAATCCAGCCCCTGCGCACGCCgggctt                                                                                                                                                                                             >  1:13450/1‑248 (MQ=255)
                                                             aCCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCCGAGCAAACTGGGTGCGGCGATG‑GTGCTTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGGTCACTCACCGCCCGGGCTCC‑GCGCCTTGCCCATGAAGCGCGACTGGAGGTCGCCCCGCTGGGG‑CAAATCCCGCACACGCGCACCCCGGGTGTGCTGGGGATggg                                                                                                                                                                                 >  2:544003/1‑249 (MQ=255)
                                                               ccctcAGGTCGTAGCGGCTGGCGTCGGTAGGTTCGGGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAATCTGGGTGCGTCGATT‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGACTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGacat                                                                                                                                                                              <  2:516705/247‑4 (MQ=255)
                                                               cACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGacat                                                                                                                                                                              >  1:516705/1‑248 (MQ=255)
                                                               cACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCATGGGCTGGATCTACAGCGTCTGAGCCAATACCCGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGGGGACGATTATCAGG‑TGATTCGTCTGGCAGGTTCATCCACCGCCCGGGGTAC‑GCGCCAGGCCCATGAAGCGCAACAGGCTGTCGCACCGCGGGGGCAAAATCC‑GCCCCTGCGCTCGCCGGGCTTgtcggcgttggacc                                                                                                                                                                               >  2:305300/1‑237 (MQ=255)
                                                                aCGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATg                                                                                                                                                                             <  1:583119/251‑1 (MQ=255)
                                                                        ggTAGCGGATGGCTGCTTGATGTGCGTGGGCTGGATAAACAGCGGCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCTCCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGTATTCCAc                                                                                                                                                                     <  2:601993/250‑1 (MQ=255)
                                                                               gctggctgCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCCAACTGGGTGCGGCGATG‑GCGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACC‑GCCATc                                                                                                                                                             >  1:330072/1‑251 (MQ=255)
                                                                                 tggctgCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACC‑GCCATcc                                                                                                                                                            >  1:341839/1‑250 (MQ=255)
                                                                                     tgctgTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑CAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGat                                                                                                                                                                          >  1:222322/1‑232 (MQ=255)
                                                                                      gctgTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACC‑GCCATCCAGagt                                                                                                                                                       <  1:394999/250‑3 (MQ=255)
                                                                                           aTGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACC‑GCCATCCAGATTGAATg                                                                                                                                                  <  1:15169/250‑1 (MQ=255)
                                                                                                         gATAAACCGCGTCTGACCCAATAAAAGAGCAAACTGGGTGCGGCGATG‑GGGATTGTTGCCGCCTGGGGCATGGAAGATTATCAGG‑TGATTCGCCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAAATGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTGTGCTGGAGATGGACATGGATTCCACC‑GCCATCCAGATTGAATGTCTTGATGAAATTGc                                                                                                                                   >  2:269002/1‑251 (MQ=255)
                                                                                                           tAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACC‑GCCATCCAGATTGAATGTATTGATGAAATTGcc                                                                                                                                  >  1:287203/1‑250 (MQ=255)
                                                                                                           tAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑GGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGGGATGGACATGGATTCCCCC‑GCCATCCAGATTGAATGTATTGATGAAATTGccc                                                                                                                                 >  2:141575/1‑250 (MQ=255)
                                                                                                            aaacagcatCTGACCCACTACCAGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATg                                                                                                                                                                             >  2:89365/9‑207 (MQ=255)
                                                                                                            aaaCAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATg                                                                                                                                                                             <  1:89365/207‑1 (MQ=255)
                                                                                                            aaaCAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCAcc                                                                                                                                                                    >  1:70579/1‑216 (MQ=255)
                                                                                                            aaaCAGCGTCTGAACCAAGACCAGAGCAAACTGGGTGCGGCGATT‑GTGATTTGTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCAcc                                                                                                                                                                    <  2:70579/216‑1 (MQ=255)
                                                                                                                 gcGTCTGACCCAATTCCAGAGCAAACTGGGTGCGGCGATT‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACC‑GCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTgg                                                                                                                           <  1:474025/250‑1 (MQ=255)
                                                                                                                  cGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACC‑GCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGc                                                                                                                          <  1:439117/251‑1 (MQ=255)
                                                                                                                  cGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACC‑GCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGc                                                                                                                          <  1:337986/251‑1 (MQ=255)
                                                                                                                      tGACCCAATACCAGAGCAAACTTGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCGCTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCTCCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACC‑GCCATCCAGATTGAATg                                                                                                                                                  <  2:376899/223‑1 (MQ=255)
                                                                                                                      tGACCCAATACCAGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACC‑GCCATCCAGATTGAATg                                                                                                                                                  >  1:376899/1‑223 (MQ=255)
                                                                                                                       gACCCAACACCAGACCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCATCTGGATGGCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACC‑GCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCg                                                                                                                        <  1:450520/248‑1 (MQ=255)
                                                                                                                         cccAATACCAGATCAAACTGGTTGCGGCATTG‑GTGAGTGTTGCCGCCTGGTGCGTTTAAGATTATCAGG‑TGTGTCGTCTGGCATGGTCACTCACCGCGCGGGCGAC‑GCGCCTGGCCCATTAAGCGCAACGGGATGTCGCCCCGCTGGGG‑AAAAGCCCGCTCCTGTGCACGCCGGGTTGGCTGGTGATGGACATTGATTCCTCC‑GCCATCCAGATTGATTGTATTGATGAAATTGCCAAACTGGCCGGAACg                                                                                                                   <  1:610084/251‑1 (MQ=255)
                                                                                                                          tcAATACCAGAGCAAACTGGGTGCGGCGATG‑TTGATTGGTGCCGACTGGCGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCTCTCACCGCACGGGCTAC‑GCGACTGGCCCATGAAGCGCAAATGGCTGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGACGGGCATGGGTTCCACC‑GCCCTCCAGGTTGAATGTATTGGACCAATTGCCAAACCGGCCGGAACgg                                                                                                                  >  2:167454/2‑251 (MQ=255)
                                                                                                                          ccAATACCAGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACC‑GCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACgg                                                                                                                  >  1:366310/1‑251 (MQ=255)
                                                                                                                                 cAGAGAAAATTGGGTGCGGTG‑TGTGTGATTGTTGCCGCCTGGGTCGTGGAAGATTACCAGG‑GGTTTCGTCGGGCAGTGTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAATCGCAACTGGATGTCGCCACGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACC‑GCCAGCCAGATTGAATGTATTGATGACATTGCCAAACTGGCCGGAACGGGCGAGAt                                                                                                           <  2:541493/251‑1 (MQ=255)
                                                                                                                                      cAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACC‑GCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGAtggtgg                                                                                                      >  1:114889/1‑251 (MQ=255)
                                                                                                                                       aaaCTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCCCTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGCCCTGGCTTCCACC‑GCCATCCAGATTGAATGTATTGTTGAAATTGCCAAACTGGCCGGAACGGGCGAGAgggtggg                                                                                                     >  2:223462/1‑250 (MQ=255)
                                                                                                                                        aaCTGGGTGCGGCGCAG‑GTGCTTGTTGCAGCCTGGTGCGAGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGACCATGAAGCGCAACTGGATGTCGCACCGCTGAGG‑AAATTCCCGCACCTGCGCACGCCGGGTTTGCTGGCGATGGCCATGGATTCCACC‑GCCATCCAGATTGACTGTATTGAGGAACTTGCCAAACTGGACGGAAAGGGCCGGATGGTGGCg                                                                                                    >  2:65147/1‑251 (MQ=255)
                                                                                                                                                ggggcgAGG‑GTGATTGTGGACGCCTGGTGCGTGTAAGATTATCAGG‑TGATGCGGATGGCAGTTTCCCTCACCGCACGGGCGAC‑GCGCCTGGCCAATGAAGCGCAACTGGATGTCGCACCGCTGGGG‑AAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGACATGGATTTCACC‑GCAAGCCAGAGTGAATGTATTGAGGAAATTGCCAAAATGGCCGGAACGGGCGAGATGGTGGCGGAAGTaa                                                                                             <  2:267977/248‑1 (MQ=255)
                                                                                                                                                    cgATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACC‑GCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGa                                                                                         >  1:114301/1‑250 (MQ=255)
                                                                                                                                                     gATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACC‑GCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAAc                                                                                       >  1:519704/1‑251 (MQ=255)
                                                                                                                                                       tG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACC‑GCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACgg                                                                                     >  1:200760/1‑251 (MQ=255)
                                                                                                                                                                 ttgCCGCCTGTTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCGCCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACC‑GCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATgcg                                                                             <  2:9239/250‑1 (MQ=255)
                                                                                                                                                                           gggCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGATAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACC‑GCCATCCAGATTGAATGTATTGATGAAAGTGCCAAACTGGCCGGAACGGGCGAGATGTTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCg                                                                  <  2:673057/249‑1 (MQ=255)
                                                                                                                                                                             gccgggAAGATTATCAGG‑TGATGCGTCGGGCAGGCTCTCTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCGGGGG‑AAAACCCCGCACCTGCGCACGCCGGGTTTGCTGGGGATGGACATGGATTCCACC‑GCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGAt                                                                <  2:366310/247‑1 (MQ=255)
                                                                                                                                                                              cGTGGAAGATTATCAGG‑TGATTAGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGCTTCCACC‑GCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGagcggaggagaactcgatt                                                               >  2:648411/1‑233 (MQ=255)
                                                                                                                                                                                tGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCAAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGACCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCAGGGTTTGCGGGTGAGGGACATAG‑CTTCACCCGCCAACCAGATTGAATGTATTGTTGTAATTGCCAAACGGGCCGGGACGGGAGAAAGGGGGGAGGAAGTAACCGAACGGGCTATGCGAGGCGAACACGAtttt                                                             >  2:369990/1‑251 (MQ=255)
                                                                                                                                                                                   aagcttATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑ACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGG‑ACAATCCCGCACCTGCGAACGCCGGGTTTGCTGGTGATGGATATGGACTCCACC‑GCCCTCCAGATTGCATGTATTGATGAAATTGCCAAACTGGCCGGACAGTGCGAGATGGTGGCGGAGTTAACCGAACGGGCGATGCGCGGCGACCTCGCTTTTAcc                                                          >  2:245472/5‑251 (MQ=255)
                                                                                                                                                                                      attattAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCGTGCCCATGAAGCGCAACGGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACC‑GCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTATCCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGcc                                                       <  2:114889/251‑1 (MQ=255)
                                                                                                                                                                                        tATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACC‑GCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGcc                                                       >  1:438685/1‑249 (MQ=255)
                                                                                                                                                                                            aGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACC‑GCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTg                                                 >  1:65851/1‑251 (MQ=255)
                                                                                                                                                                                                tGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGCTTCCACC‑GCCATCCAGATTGAATGTATTGATGAAATTGCAAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTAACGCCCGCCTggcc                                              >  2:4910/1‑248 (MQ=255)
                                                                                                                                                                                                         tGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACC‑GCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAACCTGCGCAGCCGcgtg                                     >  1:444659/1‑251 (MQ=255)
                                                                                                                                                                                                            cAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACC‑GCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCg                                  >  1:473570/1‑251 (MQ=255)
                                                                                                                                                                                                                ttCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCACGCTGGGG‑AAAATCCCGCACCTGCGAACGCCGGGTTTGCTGGTGATGGACATGGATTCCACC‑GCCATCCAGGTTGCATGTATTGATGTAATTGCCAAACAGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAAATCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGc                              >  1:37217/1‑251 (MQ=255)
                                                                                                                                                                                                                   actcacCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACC‑GCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGCTGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTg                            >  1:470659/1‑250 (MQ=255)
                                                                                                                                                                                                                      cacCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTTGTGATGGACATGGATTCCACC‑GCCATCCAGATTGAATGTAGTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGAGGCGCGGCGAACTCGATTTTACCGCCATCCTGCGCAGCCGCGTGGCGACGCTGAAAg                        <  1:273060/251‑1 (MQ=255)
                                                                                                                                                                                                                          gCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACC‑GCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGAtttt                                                             >  1:217075/1‑210 (MQ=255)
                                                                                                                                                                                                                          gCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACC‑GCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGAtttt                                                             <  2:217075/210‑1 (MQ=255)
                                                                                                                                                                                                                              gggCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACC‑GCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGAc                <  1:211941/251‑1 (MQ=255)
                                                                                                                                                                                                                                              cccATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACC‑GCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCa  >  1:441240/1‑250 (MQ=255)
                                                                                                                                                                                                                                              cccACGAAGCGCAGCTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACC‑GCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCa  >  2:131327/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       
GTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATG‑GTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑ACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACC‑GCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCA  >  NC_000913/4624931‑4625415

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