Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 558,767 T→G 6.5% intergenic (+13/‑207) sfmA → / → sfmC FimA homolog, function unknown/putative periplasmic pilus chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913558,7670TG 6.5% 267.6 / 10.2 77intergenic (+13/‑207)sfmA/sfmCFimA homolog, function unknown/putative periplasmic pilus chaperone
Reads supporting (aligned to +/- strand):  ref base T (33/39);  new base G (2/3);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.02e-01

CTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAG  >  NC_000913/558622‑558915
                                                                                                                                                  |                                                                                                                                                    
cTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACgttgtt                                                                                                                                                   >  2:284800/1‑149 (MQ=255)
     caaccaacCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGaa                                                                                                                                              <  3:18949/149‑1 (MQ=255)
     caaccaacCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGaa                                                                                                                                              <  2:681765/149‑1 (MQ=255)
      aaccaaccAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAAt                                                                                                                                             >  4:357103/1‑149 (MQ=255)
        ccaaccaGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATta                                                                                                                                           >  4:288621/1‑149 (MQ=255)
        ccaaccaGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACAACGTGGTTTTGattat                                                                                                                                           <  1:232411/149‑5 (MQ=255)
          aaccaGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATtata                                                                                                                                         <  1:173194/149‑1 (MQ=255)
          aaccaGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATtata                                                                                                                                         <  3:357103/149‑1 (MQ=255)
          aaccaGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATtata                                                                                                                                         <  4:132437/149‑1 (MQ=255)
          aaccaGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACAACGTGGTTTTGa                                                                                                                                               >  4:575307/1‑143 (MQ=255)
          aaccaGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACAACGTGGTTTTGa                                                                                                                                               <  3:575307/143‑1 (MQ=255)
             caGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTatatat                                                                                                                                      >  2:37165/1‑149 (MQ=255)
             caGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTatatat                                                                                                                                      >  1:345949/1‑149 (MQ=255)
              aGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATc                                                                                                                                     <  1:98819/149‑1 (MQ=255)
               gAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTCTAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCa                                                                                                                                    <  3:439727/149‑1 (MQ=255)
                  cTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGt                                                                                                                                 <  2:691816/149‑1 (MQ=255)
                   ttGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATta                                                                                                                                           <  3:171739/138‑1 (MQ=255)
                   ttGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATta                                                                                                                                           >  4:171739/1‑138 (MQ=255)
                   ttGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTc                                                                                                                                <  1:160381/149‑1 (MQ=255)
                       tCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTtat                                                                                                                            >  2:141964/1‑149 (MQ=255)
                        ccccGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTtata                                                                                                                           <  1:37165/149‑1 (MQ=255)
                             ggACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACaaa                                                                                                                      <  4:290667/149‑1 (MQ=255)
                                 cAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTaa                                                                                                                  <  4:478566/149‑1 (MQ=255)
                                   aCGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATg                                                                                                                >  1:59667/1‑149 (MQ=255)
                                   aCGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATg                                                                                                                >  3:603852/1‑149 (MQ=255)
                                        cttcATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATTCAAAACAACGTGGTTTTGattattc                                                                                                                                         <  2:685253/120‑7 (MQ=38)
                                        cttcATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATTCAAAACAACGTGGTTTTGattattc                                                                                                                                         >  1:685253/1‑114 (MQ=38)
                                        cttcATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCg                                                                                                           >  3:294857/1‑149 (MQ=255)
                                         ttcATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCgg                                                                                                          >  1:75501/1‑149 (MQ=255)
                                                ttCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTaa                                                                                                                  <  3:304932/134‑1 (MQ=255)
                                                ttCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTaa                                                                                                                  >  4:304932/1‑134 (MQ=255)
                                                   tGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAg                                                                                                >  4:549306/1‑149 (MQ=255)
                                                      aCGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCgg                                                                                             >  2:264920/1‑149 (MQ=255)
                                                      aCGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCgg                                                                                             >  3:419855/1‑149 (MQ=255)
                                                             aaaGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTAtggtg                                                                                      <  4:603852/149‑1 (MQ=255)
                                                              aaGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAAt                                                                                                                 <  2:680452/121‑1 (MQ=255)
                                                              aaGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAAt                                                                                                                 >  1:680452/1‑121 (MQ=255)
                                                               aGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGATGGc                                                                                    >  2:7157/1‑149 (MQ=255)
                                                                ggCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTaa                                                                                                                  >  2:312572/1‑118 (MQ=255)
                                                                ggCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTaa                                                                                                                  <  1:312572/118‑1 (MQ=255)
                                                                ggCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCa                                                                                   >  3:513495/1‑149 (MQ=255)
                                                                    ccGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGt                                                                               <  3:549306/149‑1 (MQ=255)
                                                                     cGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTta                                                                                                                             <  3:258174/102‑1 (MQ=255)
                                                                     cGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTta                                                                                                                             >  4:258174/1‑102 (MQ=255)
                                                                        tACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGCGGCAATGTAAAt                                                                           <  1:223908/149‑1 (MQ=255)
                                                                           aaGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGa                                                                        >  2:94574/1‑149 (MQ=255)
                                                                           aaGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGa                                                                        >  2:616457/1‑149 (MQ=255)
                                                                            aGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGaa                                                                       <  2:345949/149‑1 (MQ=255)
                                                                                  tCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATg                                                                 >  3:300518/1‑149 (MQ=255)
                                                                                    agcagGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGtt                                                               <  2:75501/149‑1 (MQ=255)
                                                                                     gcagGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACaaa                                                                                                                      >  1:607622/1‑93 (MQ=255)
                                                                                     gcagGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACaaa                                                                                                                      <  2:607622/93‑1 (MQ=255)
                                                                                        gggCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGTATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACt                                                           >  4:71407/1‑149 (MQ=255)
                                                                                         ggCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACtt                                                          <  1:94574/149‑1 (MQ=255)
                                                                                          gcaagcaaATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTGCACttt                                                         <  3:673572/149‑1 (MQ=255)
                                                                                            aagcaaATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGt                                                       <  4:294857/149‑1 (MQ=255)
                                                                                                aaaTGCTGACGCGACTTTTATTATGAGATATGTATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCa                                                   <  3:71407/149‑1 (MQ=255)
                                                                                                    gCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGcc                                               <  1:302578/149‑1 (MQ=255)
                                                                                                         cgcgACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCttt                                          <  1:616457/149‑1 (MQ=255)
                                                                                                             aCTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTttatta                                      >  1:457104/1‑149 (MQ=255)
                                                                                                               ttttATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGt                                                       >  2:674842/1‑130 (MQ=255)
                                                                                                               ttttATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGt                                                       <  1:674842/130‑1 (MQ=255)
                                                                                                               ttttATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATc                                                    <  3:342145/133‑1 (MQ=255)
                                                                                                               ttttATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATc                                                    >  4:342145/1‑133 (MQ=255)
                                                                                                               ttttATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTaaa                                    <  1:141964/149‑1 (MQ=255)
                                                                                                               ttttATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTaaa                                    >  2:479748/1‑149 (MQ=255)
                                                                                                               ttttATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTaaa                                    <  1:479748/149‑1 (MQ=255)
                                                                                                                    ttatGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGaaa                               >  3:445693/1‑149 (MQ=255)
                                                                                                                    ttatGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGaaa                               >  1:477738/1‑149 (MQ=255)
                                                                                                                      atGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAg                             <  3:326698/149‑1 (MQ=255)
                                                                                                                       tGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAgg                            <  4:419855/149‑1 (MQ=255)
                                                                                                                         agaTATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGa                          <  4:445693/149‑1 (MQ=255)
                                                                                                                               gAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAatgatg                    <  3:288621/149‑1 (MQ=255)
                                                                                                                                   aatCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGtgtt                <  4:513495/149‑1 (MQ=255)
                                                                                                                                    atCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGtgttg               >  4:381921/1‑149 (MQ=255)
                                                                                                                                      tCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGtgttgt              >  2:371730/1‑149 (MQ=255)
                                                                                                                                              aCGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTaaa                                    >  4:290809/1‑119 (MQ=255)
                                                                                                                                              aCGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTaaa                                    <  3:290809/119‑1 (MQ=255)
                                                                                                                                                  tgttTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAg  <  1:371730/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAG  >  NC_000913/558622‑558915

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