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

CTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAG  >  NC_000913/558622‑558915
                                                                                                                                                  |                                                                                                                                                    
cTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACgttgtt                                                                                                                                                   >  2:284800/1‑149 (MQ=255)
     caaccaacCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGaa                                                                                                                                              <  3:18949/149‑1 (MQ=255)
     caaccaacCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGaa                                                                                                                                              <  2:681765/149‑1 (MQ=255)
      aaccaaccAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAAt                                                                                                                                             >  4:357103/1‑149 (MQ=255)
        ccaaccaGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATta                                                                                                                                           >  4:288621/1‑149 (MQ=255)
        ccaaccaGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACAACGTGGTTTTGattat                                                                                                                                           <  1:232411/149‑5 (MQ=255)
          aaccaGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATtata                                                                                                                                         <  1:173194/149‑1 (MQ=255)
          aaccaGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATtata                                                                                                                                         <  3:357103/149‑1 (MQ=255)
          aaccaGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATtata                                                                                                                                         <  4:132437/149‑1 (MQ=255)
          aaccaGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACAACGTGGTTTTGa                                                                                                                                               >  4:575307/1‑143 (MQ=255)
          aaccaGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACAACGTGGTTTTGa                                                                                                                                               <  3:575307/143‑1 (MQ=255)
             caGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTatatat                                                                                                                                      >  2:37165/1‑149 (MQ=255)
             caGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTatatat                                                                                                                                      >  1:345949/1‑149 (MQ=255)
              aGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATc                                                                                                                                     <  1:98819/149‑1 (MQ=255)
               gAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTCTAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCa                                                                                                                                    <  3:439727/149‑1 (MQ=255)
                  cTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGt                                                                                                                                 <  2:691816/149‑1 (MQ=255)
                   ttGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATta                                                                                                                                           <  3:171739/138‑1 (MQ=255)
                   ttGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATta                                                                                                                                           >  4:171739/1‑138 (MQ=255)
                   ttGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTc                                                                                                                                <  1:160381/149‑1 (MQ=255)
                       tCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTtat                                                                                                                            >  2:141964/1‑149 (MQ=255)
                        ccccGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTtata                                                                                                                           <  1:37165/149‑1 (MQ=255)
                             ggACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACaaa                                                                                                                      <  4:290667/149‑1 (MQ=255)
                                 cAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTaa                                                                                                                  <  4:478566/149‑1 (MQ=255)
                                   aCGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATg                                                                                                                >  1:59667/1‑149 (MQ=255)
                                   aCGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATg                                                                                                                >  3:603852/1‑149 (MQ=255)
                                        cttcATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATTCAAAACAACGTGGTTTTGattattc                                                                                                                                         <  2:685253/120‑7 (MQ=38)
                                        cttcATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATTCAAAACAACGTGGTTTTGattattc                                                                                                                                         >  1:685253/1‑114 (MQ=38)
                                        cttcATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCg                                                                                                           >  3:294857/1‑149 (MQ=255)
                                         ttcATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCgg                                                                                                          >  1:75501/1‑149 (MQ=255)
                                                ttCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTaa                                                                                                                  <  3:304932/134‑1 (MQ=255)
                                                ttCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTaa                                                                                                                  >  4:304932/1‑134 (MQ=255)
                                                   tGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAg                                                                                                >  4:549306/1‑149 (MQ=255)
                                                      aCGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCgg                                                                                             >  2:264920/1‑149 (MQ=255)
                                                      aCGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCgg                                                                                             >  3:419855/1‑149 (MQ=255)
                                                             aaaGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTAtggtg                                                                                      <  4:603852/149‑1 (MQ=255)
                                                              aaGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAAt                                                                                                                 <  2:680452/121‑1 (MQ=255)
                                                              aaGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAAt                                                                                                                 >  1:680452/1‑121 (MQ=255)
                                                               aGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGATGGc                                                                                    >  2:7157/1‑149 (MQ=255)
                                                                ggCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTaa                                                                                                                  >  2:312572/1‑118 (MQ=255)
                                                                ggCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTaa                                                                                                                  <  1:312572/118‑1 (MQ=255)
                                                                ggCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCa                                                                                   >  3:513495/1‑149 (MQ=255)
                                                                    ccGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGt                                                                               <  3:549306/149‑1 (MQ=255)
                                                                     cGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTta                                                                                                                             <  3:258174/102‑1 (MQ=255)
                                                                     cGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTta                                                                                                                             >  4:258174/1‑102 (MQ=255)
                                                                        tACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGCGGCAATGTAAAt                                                                           <  1:223908/149‑1 (MQ=255)
                                                                           aaGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGa                                                                        >  2:94574/1‑149 (MQ=255)
                                                                           aaGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGa                                                                        >  2:616457/1‑149 (MQ=255)
                                                                            aGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGaa                                                                       <  2:345949/149‑1 (MQ=255)
                                                                                  tCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATg                                                                 >  3:300518/1‑149 (MQ=255)
                                                                                    agcagGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGtt                                                               <  2:75501/149‑1 (MQ=255)
                                                                                     gcagGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACaaa                                                                                                                      >  1:607622/1‑93 (MQ=255)
                                                                                     gcagGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACaaa                                                                                                                      <  2:607622/93‑1 (MQ=255)
                                                                                        gggCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGTATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACt                                                           >  4:71407/1‑149 (MQ=255)
                                                                                         ggCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACtt                                                          <  1:94574/149‑1 (MQ=255)
                                                                                          gcaagcaaATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTGCACttt                                                         <  3:673572/149‑1 (MQ=255)
                                                                                            aagcaaATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGt                                                       <  4:294857/149‑1 (MQ=255)
                                                                                                aaaTGCTGACGCGACTTTTATTATGAGATATGTATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCa                                                   <  3:71407/149‑1 (MQ=255)
                                                                                                    gCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGcc                                               <  1:302578/149‑1 (MQ=255)
                                                                                                         cgcgACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCttt                                          <  1:616457/149‑1 (MQ=255)
                                                                                                             aCTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTttatta                                      >  1:457104/1‑149 (MQ=255)
                                                                                                               ttttATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGt                                                       >  2:674842/1‑130 (MQ=255)
                                                                                                               ttttATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGt                                                       <  1:674842/130‑1 (MQ=255)
                                                                                                               ttttATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATc                                                    <  3:342145/133‑1 (MQ=255)
                                                                                                               ttttATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATc                                                    >  4:342145/1‑133 (MQ=255)
                                                                                                               ttttATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTaaa                                    <  1:141964/149‑1 (MQ=255)
                                                                                                               ttttATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTaaa                                    >  2:479748/1‑149 (MQ=255)
                                                                                                               ttttATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTaaa                                    <  1:479748/149‑1 (MQ=255)
                                                                                                                    ttatGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGaaa                               >  3:445693/1‑149 (MQ=255)
                                                                                                                    ttatGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGaaa                               >  1:477738/1‑149 (MQ=255)
                                                                                                                      atGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAg                             <  3:326698/149‑1 (MQ=255)
                                                                                                                       tGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAgg                            <  4:419855/149‑1 (MQ=255)
                                                                                                                         agaTATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGa                          <  4:445693/149‑1 (MQ=255)
                                                                                                                               gAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAatgatg                    <  3:288621/149‑1 (MQ=255)
                                                                                                                                   aa‑tCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGtgtt                <  4:513495/149‑1 (MQ=255)
                                                                                                                                    a‑tCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGtgttg               >  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)
                                                                                                                                                  |                                                                                                                                                    
CTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGTAAGAAACCAAAG  >  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: ‑