Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,017,222 Δ1 bp 7.6% intergenic (+30/‑111) udp → / → rmuC uridine phosphorylase/DNA recombination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,017,2220T. 7.6% 633.0 / 41.9 183intergenic (+30/‑111)udp/rmuCuridine phosphorylase/DNA recombination protein
Reads supporting (aligned to +/- strand):  ref base T (87/82);  new base . (7/7);  total (94/89)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-01

Only 100 of 183 total aligned reads displayed.

GTAGCGGGTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGC  >  NC_000913/4017076‑4017367
                                                                                                                                                   |                                                                                                                                                 
gTAGCGGGTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTc                                                                                                                                                 >  3:215418/1‑149 (MQ=255)
     gggTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCt                                                                                                                                            >  1:671028/1‑149 (MQ=255)
      ggTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCtt                                                                                                                                           <  2:896890/149‑1 (MQ=255)
      ggTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCtt                                                                                                                                           >  4:872388/1‑149 (MQ=255)
      ggTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCtt                                                                                                                                           >  4:736924/1‑149 (MQ=255)
       gtgtTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttc                                                                                                                                          <  4:236313/149‑2 (MQ=255)
        tgtTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCtttt                                                                                                                                         <  4:1025894/149‑1 (MQ=255)
        tgtTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCtttt                                                                                                                                         >  4:131106/1‑149 (MQ=255)
          ttATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGt                                                                                                                                       >  3:1433184/1‑149 (MQ=255)
          ttATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGt                                                                                                                                       <  3:85936/149‑1 (MQ=255)
           tATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTa                                                                                                                                      >  3:609176/1‑149 (MQ=255)
             tCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTAtt                                                                                                                                    <  2:950027/149‑1 (MQ=255)
             tCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTAtt                                                                                                                                    >  2:368124/1‑149 (MQ=255)
                ttAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTAttttt                                                                                                                                 >  4:655585/1‑149 (MQ=255)
                ttAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTAttttt                                                                                                                                 >  2:493917/1‑149 (MQ=255)
                ttAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTAttttt                                                                                                                                 <  3:776253/149‑1 (MQ=255)
                 tAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTg                                                                                                                                <  3:355255/149‑1 (MQ=255)
                    ccGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGt                                                                                                                             >  4:507800/1‑149 (MQ=255)
                       cACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAgc                                                                                                                          <  4:684310/149‑1 (MQ=255)
                            aGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTc                                                                                                                     <  4:609176/149‑1 (MQ=255)
                            aGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTc                                                                                                                     >  3:1029217/1‑149 (MQ=255)
                             gCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttc                                                                                                                                          >  4:578995/1‑126 (MQ=255)
                             gCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttc                                                                                                                                          <  3:578995/127‑2 (MQ=255)
                               aaGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCa                                                                                                                  <  3:1136122/149‑1 (MQ=255)
                                   gaTCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGaa                                                                                                              >  1:960625/1‑149 (MQ=255)
                                    aTCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCtttt                                                                                                                                         <  1:1074934/121‑1 (MQ=255)
                                    aTCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCtttt                                                                                                                                         >  2:1074934/1‑121 (MQ=255)
                                       ccGAATGCTGAGCCGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGc                                                                                                          >  3:523678/1‑149 (MQ=255)
                                       ccGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGc                                                                                                          <  1:336741/149‑1 (MQ=255)
                                       ccGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGc                                                                                                          <  4:1433184/149‑1 (MQ=255)
                                       ccGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCAGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGc                                                                                                          <  4:1130052/149‑1 (MQ=255)
                                           aTGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCttt                                                                                                      <  3:1221890/149‑1 (MQ=255)
                                            tGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTc                                                                                                     <  2:1355980/149‑1 (MQ=255)
                                               tGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCTTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCaa                                                                                                  >  4:624622/1‑149 (MQ=255)
                                                gagaCGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAAc                                                                                                 <  4:1235119/149‑1 (MQ=255)
                                                     gATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGAc                                                                                            >  2:931515/1‑149 (MQ=255)
                                                       tGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGt                                                                                          <  4:432378/149‑1 (MQ=255)
                                                        gAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCCACTGGACGtt                                                                                         >  2:624389/1‑149 (MQ=255)
                                                        gAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGtt                                                                                         >  4:35271/1‑149 (MQ=255)
                                                         aaacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGttt                                                                                        <  3:768136/149‑1 (MQ=255)
                                                         aaacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACCGGACGttt                                                                                        <  4:330449/149‑1 (MQ=255)
                                                            caaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTa                                                                                     <  2:1537794/149‑1 (MQ=255)
                                                             aaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGcc                                                                                                         >  3:989997/1‑128 (MQ=255)
                                                             aaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTAc                                                                                    >  1:761563/1‑149 (MQ=255)
                                                             aaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTAc                                                                                    >  4:503772/1‑149 (MQ=255)
                                                              aacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACa                                                                                   >  2:575118/1‑149 (MQ=255)
                                                               acCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAg                                                                                  >  3:1210229/1‑149 (MQ=255)
                                                                ccGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGttt                                                                                        >  3:945004/1‑142 (MQ=255)
                                                                 cGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGca                                                                                >  4:220393/1‑149 (MQ=255)
                                                                        cATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCt                                                                         >  3:49622/1‑149 (MQ=255)
                                                                           gCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTAtt                                                                      <  2:1154053/149‑1 (MQ=255)
                                                                           gCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTAtt                                                                      <  4:853669/149‑1 (MQ=255)
                                                                                 aaaaTCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGc                                                                <  2:1426770/149‑1 (MQ=255)
                                                                                   aaTCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCgg                                                              >  1:1111232/1‑149 (MQ=255)
                                                                                   aaTCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCgg                                                              <  2:960625/149‑1 (MQ=255)
                                                                                        tggtggAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAg                                                         >  1:1175608/1‑149 (MQ=255)
                                                                                            ggAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAgttgt                                                     >  1:149470/1‑149 (MQ=255)
                                                                                             gAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAgttgtt                                                    <  3:1252903/149‑1 (MQ=255)
                                                                                                  ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTc                                               <  3:872388/149‑1 (MQ=255)
                                                                                                   gcgcGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCa                                              <  2:172570/149‑1 (MQ=255)
                                                                                                       gtcgtcTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCaggag                                          >  2:234418/1‑149 (MQ=255)
                                                                                                            ctgctgTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGtt                                     <  1:183999/149‑1 (MQ=255)
                                                                                                               ctgTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTtgtg                                  <  2:671028/149‑1 (MQ=255)
                                                                                                                tgTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTgg                                 >  1:679444/1‑149 (MQ=255)
                                                                                                                 gTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGa                                <  4:179435/149‑1 (MQ=255)
                                                                                                                    aTTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCGACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGAttt                             <  1:1403070/149‑1 (MQ=255)
                                                                                                                      tctcttctctTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTCCGTCAGGAGGCGTTGTGGATTtct                           <  3:1324180/149‑1 (MQ=255)
                                                                                                                      tctcttctcCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGTGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTtct                           >  1:1288746/1‑149 (MQ=255)
                                                                                                                      tctcttctcCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTtct                           <  1:1150502/149‑1 (MQ=255)
                                                                                                                        tcttctcCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAgttgtt                                                    >  3:888470/1‑122 (MQ=255)
                                                                                                                        tcttctcCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCa                         >  2:207135/1‑149 (MQ=255)
                                                                                                                        tcttctcCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCa                         >  1:1223115/1‑149 (MQ=255)
                                                                                                                        tcttctcCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCa                         >  2:655138/1‑149 (MQ=255)
                                                                                                                         cttctcCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCaa                        <  3:460358/149‑1 (MQ=255)
                                                                                                                          ttctcCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAAt                       >  3:656449/1‑149 (MQ=255)
                                                                                                                             tcCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGTGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCAt                    <  2:1277561/149‑1 (MQ=255)
                                                                                                                              ccTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATg                   <  1:794819/149‑1 (MQ=255)
                                                                                                                               cTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCaggag                                          <  2:1207087/125‑1 (MQ=255)
                                                                                                                               cTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCaggag                                          >  1:1207087/1‑125 (MQ=255)
                                                                                                                               cTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATgg                  >  4:1365573/1‑149 (MQ=255)
                                                                                                                                tGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGt                 >  1:915274/1‑149 (MQ=255)
                                                                                                                                 gTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGtt                >  2:734792/1‑149 (MQ=255)
                                                                                                                                    tGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCt                                                                                                        >  4:1109422/1‑58 (MQ=255)
                                                                                                                                    aaaaGGCCGAACGCG‑TCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTAc             >  1:447455/3‑149 (MQ=255)
                                                                                                                                     gAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACg            >  3:68384/1‑149 (MQ=255)
                                                                                                                                     aaaGGCCGAACGCG‑TCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTa                                                                        >  1:12253/2‑89 (MQ=255)
                                                                                                                                     aaaGGCCGAACGCG‑TCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTa                                                                        <  2:12253/88‑1 (MQ=255)
                                                                                                                                      aaGGCCGAACGCG‑TCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCag                                             <  1:308622/115‑1 (MQ=255)
                                                                                                                                      aaGGCCGAACGCG‑TCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCag                                             >  2:308622/1‑115 (MQ=255)
                                                                                                                                      aaGGCCGAACGCG‑TCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCa                                              >  1:623527/1‑114 (MQ=255)
                                                                                                                                      aaGGCCGAACGCG‑TCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCa                                              <  2:623527/114‑1 (MQ=255)
                                                                                                                                      aaGGCCGAACGCG‑TCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGc           >  3:1310104/1‑149 (MQ=255)
                                                                                                                                      aaGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGc           >  1:1519232/1‑149 (MQ=255)
                                                                                                                                       aGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCa          >  4:925349/1‑149 (MQ=255)
                                                                                                                                        ggCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACATGTCTATTTGGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAACCATGGTTTACGCAg         <  2:255600/149‑1 (MQ=255)
                                                                                                                                           cG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTa      >  1:1579657/1‑149 (MQ=255)
                                                                                                                                           cG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTa      >  3:609376/1‑149 (MQ=255)
                                                                                                                                              aCGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTAtt    <  1:734792/149‑1 (MQ=255)
                                                                                                                                                gcgTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGc  <  1:953088/149‑1 (MQ=255)
                                                                                                                                                  gTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACaa                                                                             >  4:982565/1‑72 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GTAGCGGGTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGC  >  NC_000913/4017076‑4017367

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