Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,041,937 T→G 7.0% intergenic (+22/+93) appA → / ← etk phosphoanhydride phosphorylase/tyrosine‑protein kinase, role in O‑antigen capsule formation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,041,9370TG 7.0% 202.1 / 4.7 75intergenic (+22/+93)appA/etkphosphoanhydride phosphorylase/tyrosine‑protein kinase, role in O‑antigen capsule formation
Reads supporting (aligned to +/- strand):  ref base T (32/36);  new base G (4/2);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

CGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCG  >  NC_000913/1041791‑1042084
                                                                                                                                                  |                                                                                                                                                   
agcagccCGGAGAGGTGAAACAGAACAAAGAAAGATGAGAAGAGCGAAAGGCGCAGGGCAGGAGTTCGTTGGTAGGATATACGCAAACAGTGAATGAAGCACGCATACCGGCTTGCAGTAAAAAAAAAAAAAAAAAGAGCATTCAGTTa                                                                                                                                                   <  2:831735/145‑1 (MQ=255)
 gccgccCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTAc                                                                                                                                                  >  2:1009487/1‑149 (MQ=255)
   cgccCGGGGAGGTGAAACTGGCCCTGGCGGGATGTGAAGGGCGAAATGCGGAGGGCATGTGTTCGTTGGCAGTTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCt                                                                                                                                                <  2:63874/149‑1 (MQ=255)
   cgccCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCt                                                                                                                                                >  1:1020141/1‑149 (MQ=255)
   cgccCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCt                                                                                                                                                >  2:75352/1‑149 (MQ=255)
    gccCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAACAGAGCATTCAGGTACCTg                                                                                                                                               >  2:329243/1‑149 (MQ=255)
     cccGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGa                                                                                                                                              >  1:609812/1‑149 (MQ=255)
      ccGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGaa                                                                                                                                             <  1:1039502/149‑1 (MQ=255)
       cGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAAt                                                                                                                                            >  1:741854/1‑149 (MQ=255)
          gggggcgAAACGGGCCCTGGCAGGATGGGAAGCGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATTAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGct                                                                                                                                         <  2:41265/143‑1 (MQ=255)
          agagGTGAAACTGACCCTGGCAGGATGGGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGct                                                                                                                                         <  2:266258/149‑1 (MQ=255)
            agGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGAAAAGGTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGctct                                                                                                                                       <  1:839543/149‑1 (MQ=255)
              gTGAAACGGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGa                                                                                                                                     <  2:1124618/149‑1 (MQ=255)
               tGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGGTAACTGAATGctctttt                                                                                                                                    <  1:900884/149‑4 (MQ=255)
               ggAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAg                                                                                                                                    <  1:463212/148‑1 (MQ=255)
                gttaCTGACCCTGGCAGGATGTGGAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTGACGAAAAGCGCGAAGGGAGCACGCATACCGCCGTGCGGTTGGTAGGGCATAAAAAAGAGCACTCAGTTGCCTGACTGCTCTGAgg                                                                                                                                   >  2:591086/4‑149 (MQ=255)
                    ccgACCCTGGCAGGATGTGAAGAGCGAAATGCGCGGGGCAGGAGTTCGTTCGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTAAAAAAAAATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCtga                                                                                                                               <  2:946425/147‑1 (MQ=255)
                    cTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCtga                                                                                                                               <  1:92130/149‑1 (MQ=255)
                    cTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCtga                                                                                                                               <  1:822615/149‑1 (MQ=255)
                     tGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCtgat                                                                                                                              <  1:691537/149‑1 (MQ=255)
                       aCCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCtgatga                                                                                                                            >  1:20120/1‑149 (MQ=255)
                         ccTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGACTGATCTGAGGCTGATGACa                                                                                                                          >  2:149532/1‑149 (MQ=255)
                         ccTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACa                                                                                                                          >  2:60931/1‑149 (MQ=255)
                           tGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACaaa                                                                                                                        <  1:38946/149‑1 (MQ=255)
                           tGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACaaa                                                                                                                        <  2:916325/149‑1 (MQ=255)
                             gCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACg                                                                                                                      >  2:240135/1‑149 (MQ=255)
                               aGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACgaa                                                                                                                    <  1:1009487/149‑1 (MQ=255)
                                 gATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACgaaga                                                                                                                  >  2:538230/1‑149 (MQ=255)
                                   tgtgAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTAAGTTACCTGAATGCTCTGAGGCTGATGGCAAACgaagaaa                                                                                                                >  2:260166/1‑148 (MQ=255)
                                    gtgAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGACGCACGCATACCGGCGTGCAGTTGGTAATGCATAAAAACGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACGAACGAAGAACt                                                                                                               >  2:576038/1‑149 (MQ=255)
                                      gAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGGGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTAATGACAAACGAAGAACTGt                                                                                                             >  2:37591/1‑149 (MQ=255)
                                      gAAGAGCGAAATGAGCAGGGCATGTGTTCGTTGGCAGGTGTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGATCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAAAGAACAACTGt                                                                                                             >  2:704185/1‑149 (MQ=255)
                                       aaGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTc                                                                                                            <  1:504034/149‑1 (MQ=255)
                                           gCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAAt                                                                                                        <  2:370919/149‑1 (MQ=255)
                                              aaaTGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCg                                                                                                     >  1:365405/1‑149 (MQ=255)
                                               aaTGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGt                                                                                                    <  1:75352/149‑1 (MQ=255)
                                                 tGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAg                                                                                                  <  2:963090/149‑1 (MQ=255)
                                                          cATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGaaaa                                                                                         <  1:309859/149‑1 (MQ=255)
                                                            tgtgTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGCGGCTGATGACAAACGCGGACCTGTCTCATGTGTGGGCCGGACAAgg                                                                                       >  1:491757/1‑149 (MQ=255)
                                                               gtTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCgg                                                                                    >  1:300158/1‑148 (MQ=255)
                                                               gtTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGt                                                                                    <  1:897590/149‑1 (MQ=255)
                                                                ttCGTTGGCAGGTTATACGCAAATCGTGAAGGAAGCACGCCTACCGGCCTGCGGTTTGTAAGTAATAAAAAAGAGCATTCAGTTACCTGAATGATCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGtt                                                                                   <  1:579/149‑1 (MQ=255)
                                                                 tCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGAGTTc                                                                                  >  2:840943/1‑149 (MQ=255)
                                                                   gTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCAc                                                                                <  1:921695/149‑1 (MQ=255)
                                                                         aGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCa                                                                          >  2:1130765/1‑149 (MQ=255)
                                                                           gTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGCTTGTAATGCATAACAAGTAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATc                                                                        >  1:1023257/1‑149 (MQ=255)
                                                                              ttACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCgg                                                                     <  2:704768/149‑1 (MQ=255)
                                                                              ttACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCgg                                                                     >  1:804474/1‑149 (MQ=255)
                                                                                   cAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACt                                                                >  2:1093127/1‑149 (MQ=255)
                                                                                          tgaatgaaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGtt                                                         <  2:181718/149‑1 (MQ=255)
                                                                                           gaatgaaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGttt                                                        <  1:60931/149‑1 (MQ=255)
                                                                                           gaaggatGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGttt                                                        <  2:404013/149‑1 (MQ=255)
                                                                                             atgaaGCACGCATAGCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTa                                                      <  2:767573/149‑1 (MQ=255)
                                                                                              tgaaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACGTTCAGTTTTAc                                                     >  2:102018/1‑149 (MQ=255)
                                                                                              cgaaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTTATGACAAACGAAGAACTGTCCAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTAc                                                     <  1:614068/148‑1 (MQ=255)
                                                                                                   cACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCttt                                                <  1:1093127/149‑1 (MQ=255)
                                                                                                    aCGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTtc                                               >  1:449030/1‑149 (MQ=255)
                                                                                                         tACCGGCGTGCAGTTTGTAATGCATAAAAAATAGCATTCAGCTACCGGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGa                                          <  1:602014/149‑1 (MQ=255)
                                                                                                           ccGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAgg                                        >  2:914260/1‑148 (MQ=255)
                                                                                                            cGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTa                                       >  2:423443/1‑149 (MQ=255)
                                                                                                             cgCGTGAAGGTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTaa                                      <  1:329243/148‑1 (MQ=255)
                                                                                                              gCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATAACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTacc                                     >  2:1041109/1‑147 (MQ=255)
                                                                                                                gTGCAGTTTGTGCGGCATAAAAAAGGGCATTCCGGTACCTGAATGCTCTGAGGCTGATGACCAACGAAGAACTGTCTAATACGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTCTCAGTTTTACTCTTTCTCGGAGTAACta                                   <  2:648064/149‑1 (MQ=255)
                                                                                                                   taGTTTGTAATGCATAAAAAGGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCTCATCCGGCCACTTTCCGTTTTACTCTTTCTCGGAGTATCTCTaa                                <  2:394853/148‑1 (MQ=255)
                                                                                                                     gTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTTATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAAcc                              <  1:1125356/149‑1 (MQ=255)
                                                                                                                           aaTGCAAAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATa                        <  1:672613/149‑1 (MQ=255)
                                                                                                                                  aaaaaaGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAg                 <  2:20120/149‑1 (MQ=255)
                                                                                                                                      atgagCGTTGAGTTACCTGCATGCTGTGAGTTTGATGGCAAACGAAGGACTGTCTAGTGCGTTGACCGGAAAGGGCGTTCTGGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGt             <  2:258221/147‑1 (MQ=255)
                                                                                                                                      aaGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGt             >  2:284767/1‑149 (MQ=255)
                                                                                                                                        gagCATTCAGGTAACTGAATGCTCTGGCGATCGTCACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTaa           >  1:279657/1‑149 (MQ=255)
                                                                                                                                         agCATTCAGGTAACTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCAACTTCTGCTGGACGCGTGCACAGCGCACCTATCACCGTAATAGTTATAGCCGTGAc          >  1:1123897/1‑149 (MQ=255)
                                                                                                                                           cATTCAGGTAACTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATACCCGTCATAGTTATAACAGTAAcat        >  2:304017/1‑147 (MQ=255)
                                                                                                                                            aTTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCATTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTACCTGt       >  1:616000/1‑149 (MQ=255)
                                                                                                                                                aGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATCGTTATAGCCGTAACTGTAAGc   >  2:485910/1‑149 (MQ=255)
                                                                                                                                                 gTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCg  >  2:115632/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCG  >  NC_000913/1041791‑1042084

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