Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 558,762 C→A 6.5% intergenic (+8/‑212) 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,7620CA 6.5% 266.0 / 7.9 77intergenic (+8/‑212)sfmA/sfmCFimA homolog, function unknown/putative periplasmic pilus chaperone
Reads supporting (aligned to +/- strand):  ref base C (34/38);  new base A (2/3);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.68e-01

GAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGT  >  NC_000913/558616‑558903
                                                                                                                                                   |                                                                                                                                             
gAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCAc                                                                                                                                             >  2:173194/1‑149 (MQ=255)
      cTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACgttgtt                                                                                                                                       >  2:284800/1‑149 (MQ=255)
           caaccaacCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGaa                                                                                                                                  <  2:681765/149‑1 (MQ=255)
           caaccaacCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGaa                                                                                                                                  <  3:18949/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                                                                                                                             <  4:132437/149‑1 (MQ=255)
                aaccaGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATtata                                                                                                                             <  3:357103/149‑1 (MQ=255)
                aaccaGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATtata                                                                                                                             <  1:173194/149‑1 (MQ=255)
                aaccaGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACAACGTGGTTTTGa                                                                                                                                   >  4:575307/1‑143 (MQ=255)
                aaccaGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACAACGTGGTTTTGa                                                                                                                                   <  3:575307/143‑1 (MQ=255)
                   caGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTatatat                                                                                                                          >  1:345949/1‑149 (MQ=255)
                   caGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTatatat                                                                                                                          >  2:37165/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                                                                                                                               >  4:171739/1‑138 (MQ=255)
                         ttGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATta                                                                                                                               <  3:171739/138‑1 (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                                                                                                                             >  1:685253/1‑114 (MQ=38)
                                              cttcATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATTCAAAACAACGTGGTTTTGattattc                                                                                                                             <  2:685253/120‑7 (MQ=38)
                                              cttcATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCg                                                                                               >  3:294857/1‑149 (MQ=255)
                                               ttcATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCgg                                                                                              >  1:75501/1‑149 (MQ=255)
                                                      ttCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTaa                                                                                                      >  4:304932/1‑134 (MQ=255)
                                                      ttCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTaa                                                                                                      <  3:304932/134‑1 (MQ=255)
                                                         tGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAg                                                                                    >  4:549306/1‑149 (MQ=255)
                                                            aCGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCgg                                                                                 >  3:419855/1‑149 (MQ=255)
                                                            aCGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCgg                                                                                 >  2:264920/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                                                                                                                 >  4:258174/1‑102 (MQ=255)
                                                                           cGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTta                                                                                                                 <  3:258174/102‑1 (MQ=255)
                                                                              tACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGCGGCAATGTAAAt                                                               <  1:223908/149‑1 (MQ=255)
                                                                                 aaGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGa                                                            >  2:616457/1‑149 (MQ=255)
                                                                                 aaGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGa                                                            >  2:94574/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)
                                                                                                                                                   |                                                                                                                                             
GAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGT  >  NC_000913/558616‑558903

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: