Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,017,217:1 +A 15.6% intergenic (+25/‑116) 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,2161.A15.6% 674.6 / 141.1 212intergenic (+24/‑117)udp/rmuCuridine phosphorylase/DNA recombination protein
Reads supporting (aligned to +/- strand):  ref base . (71/108);  new base A (16/17);  total (87/125)
Fisher's exact test for biased strand distribution p-value = 3.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01

Only 100 of 216 total aligned reads displayed.

AATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTG  >  NC_000913/4017118‑4017315
                                                                                                   |                                                                                                   
    cTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑Acgcg                                                                                                >  1:4267070/1‑100 (MQ=255)
     tGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGt                                                                                               <  2:3970575/100‑1 (MQ=255)
     tGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGt                                                                                               >  1:2842392/1‑100 (MQ=255)
       agaCGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTc                                                                                             <  2:1902269/100‑1 (MQ=255)
       agaCGATGAAACAAACCGAAAGCCATGCGGGGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTc                                                                                             <  1:1276840/100‑1 (MQ=255)
        gaCGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCg                                                                                            <  1:2145209/100‑1 (MQ=255)
        gaCGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCg                                                                                            <  2:504127/100‑1 (MQ=255)
         aCGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCgg                                                                                           <  2:1583461/100‑1 (MQ=255)
          cGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGc                                                                                          >  2:2762314/1‑100 (MQ=255)
          cGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGc                                                                                          >  2:3132832/1‑100 (MQ=255)
          cGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGCCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGc                                                                                          >  1:2679209/1‑100 (MQ=255)
           gATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGcc                                                                                         <  1:2059557/100‑1 (MQ=255)
            aTGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGGCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCt                                                                                        <  1:2181044/100‑1 (MQ=255)
              gAAACAAACCGAAAGCCCTGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTTAAGGTCT‑TCGCTTTTGGCCttt                                                                                      <  1:2562008/100‑1 (MQ=255)
               aaacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttca                                                                                     <  2:4044841/100‑3 (MQ=255)
               aaacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttca                                                                                     <  2:3277903/100‑3 (MQ=255)
               aaacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCtttt                                                                                     >  2:90020/1‑100 (MQ=255)
               aaacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGCCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttc                                                                                      <  2:3361501/99‑2 (MQ=255)
                aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttcag                                                                                    <  1:1120130/100‑4 (MQ=255)
                aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttcag                                                                                    <  2:54523/100‑4 (MQ=255)
                aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTg                                                                                    <  2:467687/100‑1 (MQ=255)
                aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTg                                                                                    >  2:2946362/1‑100 (MQ=255)
                 acaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttca                                                                                     <  1:1002661/98‑3 (MQ=255)
                 acaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGt                                                                                   >  2:1497494/1‑100 (MQ=255)
                  caaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTa                                                                                  >  2:4087133/1‑100 (MQ=255)
                   aaagCGAAAGGCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttca                                                                                     <  1:4053060/96‑3 (MQ=255)
                    aacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTAtt                                                                                <  2:1989390/100‑1 (MQ=255)
                     acCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttcagacagg                                                                               <  1:2186986/100‑9 (MQ=38)
                     acCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttcagacagg                                                                               <  2:1433323/100‑9 (MQ=38)
                     acCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTAttt                                                                               <  2:987399/100‑1 (MQ=255)
                      ccGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTAtttt                                                                              <  1:526254/100‑1 (MQ=255)
                       cGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTAttttt                                                                             <  2:801979/100‑1 (MQ=255)
                        gacaggcATGCGGTGAAAATCGTGGTGGAAGCTGCGCGTCGTCTGCGGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttcag                                                                                    <  1:2765090/86‑4 (MQ=37)
                        gAAAGCCATGCGGTGGAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTg                                                                            <  2:4337536/100‑1 (MQ=255)
                        gAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTg                                                                            >  1:3364677/1‑100 (MQ=255)
                         aaaGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTAGGCCTTTTGTATTTTTGc                                                                           <  2:1817559/100‑1 (MQ=255)
                          aaGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCg                                                                          <  2:2351399/100‑1 (MQ=255)
                          aaGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGAGTCGTCTGCTCTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTAtttttt                                                                            <  1:679246/98‑2 (MQ=255)
                            gCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTa                                                                        >  2:296077/1‑100 (MQ=255)
                               aTGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAgcg                                                                     <  1:1803983/100‑1 (MQ=255)
                                   ggTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTc                                                                 >  2:1213055/1‑100 (MQ=255)
                                    gTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCg                                                                <  2:3943569/100‑1 (MQ=255)
                                     tGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGc                                                               <  1:2496065/100‑1 (MQ=255)
                                     tGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGc                                                               >  2:1082944/1‑100 (MQ=255)
                                       aaaaTCGTGGTGGAAGCGGCGCGTCGTCTGCTGTTATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAg                                                             <  2:3545150/100‑1 (MQ=255)
                                        aaaTCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAgg                                                            >  2:2246596/1‑100 (MQ=255)
                                        aaaTCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAg                                                             <  2:1010254/99‑1 (MQ=255)
                                             gtggtgGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTCTTTTTGCGTAGCGCCTCGCAGGAAATg                                                       >  1:1605809/1‑100 (MQ=255)
                                               ggtggAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGcc                                                     <  1:3132815/100‑1 (MQ=255)
                                                  ggAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCttt                                                  >  2:1376648/1‑100 (MQ=255)
                                                   gAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTc                                                 <  2:1303760/100‑1 (MQ=255)
                                                     aGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCa                                               <  1:3313740/100‑1 (MQ=255)
                                                       cggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAAc                                             <  2:3327887/100‑1 (MQ=255)
                                                        ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAAc                                             <  1:46224/99‑1 (MQ=255)
                                                        ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACt                                            <  2:3832228/100‑1 (MQ=255)
                                                        ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACt                                            >  2:2646632/1‑100 (MQ=255)
                                                        ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACt                                            <  1:131093/100‑1 (MQ=255)
                                                        ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACt                                            >  1:106515/1‑100 (MQ=255)
                                                         gcgcGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGGAGGAAATGCCTTTCCAACTg                                           <  2:4140972/100‑1 (MQ=255)
                                                         gcgcGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTg                                           <  2:2454145/100‑1 (MQ=255)
                                                           gcgTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGa                                         <  1:3247104/100‑1 (MQ=255)
                                                            cgtcgtCTGCTGTAATTCTCTTCTCCTGTCGGAAGGCCG‑ACGCGTTCGGCCTTTTGGATTTTTGCGTAGCGCCTCGCAGGAAAGGCATTTcctgcgaggc                                        <  2:1690550/100‑9 (MQ=255)
                                                             gtcgtcTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACg                                       <  1:933292/100‑1 (MQ=255)
                                                                 tcTGCTGTAATTCTCTTCTCCTGTCTGAACGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTg                                   <  1:2641784/100‑1 (MQ=255)
                                                                  ctgctgTAATTCTCTTCTCCTGTCTGCAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGt                                  >  1:3606387/1‑100 (MQ=255)
                                                                  ctgctgTAATTCTCTTCTCCTGTCTGAAGGCTG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGt                                  <  1:2971813/100‑1 (MQ=255)
                                                                  cagctgTAATTCTATTCTCCTGTCTGAAGGCAG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGt                                  >  2:2971829/3‑100 (MQ=255)
                                                                    gctgTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTAc                                <  2:3199198/100‑1 (MQ=255)
                                                                    gctgTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTAc                                >  1:699425/1‑100 (MQ=255)
                                                                     ctgTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACa                               <  1:3453417/100‑1 (MQ=255)
                                                                      tgTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACa                               >  2:3429626/1‑99 (MQ=255)
                                                                      tgTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAg                              <  1:1458241/100‑1 (MQ=255)
                                                                       gTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGc                             <  1:1592031/100‑1 (MQ=255)
                                                                       gTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGc                             <  2:3472214/100‑1 (MQ=255)
                                                                        tAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGca                            <  2:917167/100‑1 (MQ=255)
                                                                          aTTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGcaca                          >  2:28403/1‑100 (MQ=255)
                                                                             ctcttctcCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAAtt                       >  2:2464929/1‑100 (MQ=255)
                                                                               cttctcCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCt                     >  1:2400837/1‑100 (MQ=255)
                                                                                   tcCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTAttt                 <  1:148555/100‑1 (MQ=255)
                                                                                   aaaatacaaaaGGCCGAACGCG‑TCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTAttt                 >  2:3209750/10‑100 (MQ=37)
                                                                                    ccTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTAtttt                >  1:1157482/1‑100 (MQ=255)
                                                                                     cTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTAtttt                <  2:3738318/99‑1 (MQ=255)
                                                                                       gTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTGGTACAGCACAAATCTATTTtgtg             >  1:422899/1‑100 (MQ=255)
                                                                                         caaaaGGCCGAACGCG‑TCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCg           >  1:3501540/4‑100 (MQ=255)
                                                                                         caaaaGGCCGAACGCG‑TCGGCCTATTGTATTTTTGAGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGAAGAGTGTACAACACAATTCTATTTtgtgag           >  1:3211728/4‑98 (MQ=255)
                                                                                          tGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTTCTACTGGTCGTTTGTACCGCACAATTCTATTTTGTGCgg          >  1:217717/1‑100 (MQ=255)
                                                                                          tGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCgg          >  2:2924966/1‑100 (MQ=255)
                                                                                          tGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCgg          >  2:1936913/1‑100 (MQ=255)
                                                                                          tGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCg           >  1:768899/1‑99 (MQ=255)
                                                                                           aaaGGCCGAACGCG‑TCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTTCggg         >  1:4018294/2‑100 (MQ=255)
                                                                                             aGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTa       <  2:2282108/100‑1 (MQ=255)
                                                                                              ggCCGAACGCG‑TCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTaa      >  2:659148/1‑100 (MQ=255)
                                                                                              ggCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTaa      <  1:3397003/100‑1 (MQ=255)
                                                                                              ggCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTa       >  1:1428646/1‑99 (MQ=255)
                                                                                               gCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAg     <  1:1213050/100‑1 (MQ=255)
                                                                                               gCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAg     <  2:1020549/100‑1 (MQ=255)
                                                                                                ccG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAgt    >  1:3859531/1‑100 (MQ=255)
                                                                                                  g‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAgttg  >  1:902429/1‑100 (MQ=255)
                                                                                                  g‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAACACAATTCTATTTTGTGCGGGTAAgttg  <  2:4238284/100‑1 (MQ=255)
                                                                                                  g‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGAAGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAgttg  >  2:154992/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
AATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTG  >  NC_000913/4017118‑4017315

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