Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,041,937 T→G 8.8% 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 8.8% 521.4 / 36.2 158intergenic (+22/+93)appA/etkphosphoanhydride phosphorylase/tyrosine‑protein kinase, role in O‑antigen capsule formation
Reads supporting (aligned to +/- strand):  ref base T (69/75);  new base G (4/10);  total (73/85)
Fisher's exact test for biased strand distribution p-value = 2.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 160 total aligned reads displayed.

TACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGG  >  NC_000913/1041789‑1042085
                                                                                                                                                    |                                                                                                                                                    
  cgccgcCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTa                                                                                                                                                    <  3:883847/149‑1 (MQ=255)
   gccgccCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTAc                                                                                                                                                   <  1:403548/149‑1 (MQ=255)
    ccgccCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTAcc                                                                                                                                                  >  3:321727/1‑149 (MQ=255)
     cgccCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCt                                                                                                                                                 >  2:821324/1‑149 (MQ=255)
      gccCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTg                                                                                                                                                >  2:626825/1‑149 (MQ=255)
         cGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCATGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGGTAACTGAAt                                                                                                                                             <  1:1116926/149‑1 (MQ=255)
         cGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGGTAACTGAAt                                                                                                                                             <  1:1009241/149‑1 (MQ=255)
          ggAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGGTAACTGAATg                                                                                                                                            <  4:964791/149‑1 (MQ=255)
          ggAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGGTAACTGAATg                                                                                                                                            <  3:300034/149‑1 (MQ=255)
           gagaGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTAGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGGTAACTGAATGc                                                                                                                                           <  1:619047/149‑1 (MQ=255)
             gagGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGctc                                                                                                                                         <  2:311263/149‑1 (MQ=255)
             gagGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGctc                                                                                                                                         <  2:591242/149‑1 (MQ=255)
             gagGTGAAACTGAACCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGctc                                                                                                                                         <  3:232145/149‑1 (MQ=255)
              agGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGctct                                                                                                                                        <  3:1206173/149‑1 (MQ=255)
              agGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGctct                                                                                                                                        <  2:908920/149‑1 (MQ=255)
               ggTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGGTAACTGAATGctctt                                                                                                                                       <  1:1522619/149‑2 (MQ=255)
                 tGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAg                                                                                                                                     >  3:684672/1‑149 (MQ=255)
                 tGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAg                                                                                                                                     >  2:1041907/1‑149 (MQ=255)
                 tGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGGTAACTGAATGctctttt                                                                                                                                     <  4:544295/149‑4 (MQ=255)
                     aCTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCtg                                                                                                                                 >  1:1254974/1‑149 (MQ=255)
                      cTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGGGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCtga                                                                                                                                <  2:251480/149‑1 (MQ=255)
                       tGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCtgat                                                                                                                               <  2:262472/149‑1 (MQ=255)
                          cccTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGAc                                                                                                                            >  1:1383576/1‑149 (MQ=255)
                              ggCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGa                                                                                                                                      >  2:731898/1‑135 (MQ=255)
                              ggCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAAc                                                                                                                        >  4:114713/1‑149 (MQ=255)
                               gCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACg                                                                                                                       >  2:712252/1‑149 (MQ=255)
                                   gATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACgaaga                                                                                                                   >  1:348186/1‑149 (MQ=255)
                                   gATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACgaaga                                                                                                                   <  4:133069/149‑1 (MQ=255)
                                        gAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGt                                                                                                              >  3:359574/1‑149 (MQ=255)
                                         aaGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTc                                                                                                             <  3:1364371/149‑1 (MQ=255)
                                          agagCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCt                                                                                                            <  2:1688764/149‑1 (MQ=255)
                                            agCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTaa                                                                                                          <  3:1039383/149‑1 (MQ=255)
                                             gCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGctc                                                                                                                                         >  3:712121/1‑117 (MQ=255)
                                              cGAAATGCGCAGGGCATGTGTTCGTTGGCTGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATg                                                                                                        >  4:899534/1‑149 (MQ=255)
                                                        aGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCgg                                                                                              >  2:293083/1‑149 (MQ=255)
                                                            cATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGaaaa                                                                                          <  2:1068370/149‑1 (MQ=255)
                                                            cATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGaaaa                                                                                          >  4:71350/1‑149 (MQ=255)
                                                              tgtgTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAgg                                                                                        <  3:950374/149‑1 (MQ=255)
                                                               gtgtTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGc                                                                                       <  1:569677/149‑1 (MQ=255)
                                                                tgtTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCg                                                                                      <  3:1509796/149‑1 (MQ=255)
                                                                 gtTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGt                                                                                     <  1:855142/149‑1 (MQ=255)
                                                                 gtTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGt                                                                                     <  1:219816/149‑1 (MQ=255)
                                                                 gtTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGt                                                                                     <  1:1560027/149‑1 (MQ=255)
                                                                  ttCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGtt                                                                                    >  1:634712/1‑149 (MQ=255)
                                                                  ttCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGtt                                                                                    <  4:317246/149‑1 (MQ=255)
                                                                  ttCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGtt                                                                                    >  2:575917/1‑149 (MQ=255)
                                                                  ttCGTCGGCAGGTTTTACGCAAATCGCGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGtt                                                                                    >  4:1440973/1‑149 (MQ=255)
                                                                         gCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCAcgccg                                                                             >  1:1682694/1‑149 (MQ=255)
                                                                             gTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAGTGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATc                                                                         >  2:1011191/1‑149 (MQ=255)
                                                                              ttttACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATcc                                                                        >  3:54632/1‑149 (MQ=255)
                                                                                 tACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGc                                                                     <  2:1254974/149‑1 (MQ=255)
                                                                                    gCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCAc                                                                  <  2:634712/149‑1 (MQ=255)
                                                                                    gCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCAc                                                                  >  4:1256988/1‑149 (MQ=255)
                                                                                    gCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCAc                                                                  >  4:1396814/1‑149 (MQ=255)
                                                                                     cAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACt                                                                 >  4:1152701/1‑149 (MQ=255)
                                                                                      aaaTCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACtt                                                                <  1:293083/149‑1 (MQ=255)
                                                                                      aaaTCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACtt                                                                >  4:252322/1‑149 (MQ=255)
                                                                                      aaaTCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACtt                                                                >  3:908742/1‑149 (MQ=255)
                                                                                        aTCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTc                                                              >  3:808748/1‑149 (MQ=255)
                                                                                         tCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCa                                                             <  1:1316281/149‑1 (MQ=255)
                                                                                           gTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGt                                                           >  2:1707004/1‑149 (MQ=255)
                                                                                            tgaatgaaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGtt                                                          <  4:117628/149‑1 (MQ=255)
                                                                                            tgaatgaaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGtt                                                          >  1:109929/1‑149 (MQ=255)
                                                                                            tgaatgaaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGtt                                                          <  2:1535364/149‑1 (MQ=255)
                                                                                               atgaaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTa                                                       <  4:1346724/149‑1 (MQ=255)
                                                                                               atgaaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTa                                                       >  1:1092246/1‑149 (MQ=255)
                                                                                               atgaaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTa                                                       >  2:1362234/1‑149 (MQ=255)
                                                                                               atgaaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTa                                                       <  1:1362637/149‑1 (MQ=255)
                                                                                                    gCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACgaaga                                                                                                                   <  3:1173607/84‑1 (MQ=255)
                                                                                                    gCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCtt                                                  >  2:999208/1‑149 (MQ=255)
                                                                                                    gCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCtt                                                  >  3:804159/1‑149 (MQ=255)
                                                                                                      aCGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGTCCGGAAAAGGCGTTCACGCCGCGTCCGGCCACTTTCACTTTTACTCTTtc                                                <  3:514299/149‑1 (MQ=255)
                                                                                                         cATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCg                                             <  1:712252/149‑1 (MQ=255)
                                                                                                          aTACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCgg                                            >  3:317868/1‑149 (MQ=255)
                                                                                                          aTACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCg                                             <  3:1245126/148‑1 (MQ=255)
                                                                                                          aTACCGGCGTGCAGTGTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCg                                             >  4:1245126/1‑148 (MQ=255)
                                                                                                           ggtatgCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGa                                           <  1:592079/144‑1 (MQ=255)
                                                                                                           ggtatgCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGa                                           <  3:1402410/144‑1 (MQ=255)
                                                                                                                     cAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATaa                                 <  3:71350/149‑1 (MQ=255)
                                                                                                                       gTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAAcc                               <  1:1723120/149‑1 (MQ=255)
                                                                                                                        tttGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCa                                                             <  1:92659/118‑1 (MQ=255)
                                                                                                                        tttGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCg                              >  2:182661/1‑149 (MQ=255)
                                                                                                                           gTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTaa                           <  3:493958/149‑1 (MQ=255)
                                                                                                                            tAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAAt                          <  1:766646/149‑1 (MQ=255)
                                                                                                                             aaTGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCGCGGAGTAACTATAACCGTAATa                         <  1:1041907/149‑1 (MQ=255)
                                                                                                                                   tAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAAt                          >  3:1534973/1‑142 (MQ=255)
                                                                                                                                   tAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAAt                          <  4:1534973/142‑1 (MQ=255)
                                                                                                                                    aaaaaaGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTaa                                                                                                          >  4:536268/1‑61 (MQ=255)
                                                                                                                                    aaaaaaGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAAcc                               >  2:1274695/1‑136 (MQ=255)
                                                                                                                                    aaaaaaGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAg                  >  1:1433725/1‑149 (MQ=255)
                                                                                                                                    aaaaaaGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAg                  >  2:517383/1‑149 (MQ=255)
                                                                                                                                     aaaaaGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGc                 >  3:1336080/1‑149 (MQ=255)
                                                                                                                                      ctcagagCATTCAGGTAACTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGcc                >  2:666457/4‑149 (MQ=255)
                                                                                                                                        aaGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTtc                                                >  3:1569910/1‑115 (MQ=255)
                                                                                                                                        aaGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTtc                                                <  4:1569910/115‑1 (MQ=255)
                                                                                                                                         agagCATTCAGGTAACTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGa                                           >  1:1265171/1‑119 (MQ=255)
                                                                                                                                         agagCATTCAGGTAACTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTa             >  3:570425/1‑149 (MQ=255)
                                                                                                                                             cATTCAGGTAACTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGt                                                                                     >  2:1379687/1‑73 (MQ=39)
                                                                                                                                                   gTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCg   <  2:1510356/149‑1 (MQ=255)
                                                                                                                                                    ttACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCgg  <  2:1175265/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGG  >  NC_000913/1041789‑1042085

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: