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

GAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGT  >  NC_000913/558616‑558903
                                                                                                                                                   |                                                                                                                                             
gAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCAc                                                                                                                                             >  2:173194/1‑149 (MQ=255)
      cTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACgttgtt                                                                                                                                       >  2:284800/1‑149 (MQ=255)
           caaccaacCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGaa                                                                                                                                  <  2:681765/149‑1 (MQ=255)
           caaccaacCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGaa                                                                                                                                  <  3:18949/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                                                                                                                             <  4:132437/149‑1 (MQ=255)
                aaccaGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATtata                                                                                                                             <  3:357103/149‑1 (MQ=255)
                aaccaGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATtata                                                                                                                             <  1:173194/149‑1 (MQ=255)
                aaccaGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACAACGTGGTTTTGa                                                                                                                                   >  4:575307/1‑143 (MQ=255)
                aaccaGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACAACGTGGTTTTGa                                                                                                                                   <  3:575307/143‑1 (MQ=255)
                   caGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTatatat                                                                                                                          >  1:345949/1‑149 (MQ=255)
                   caGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTatatat                                                                                                                          >  2:37165/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                                                                                                                               >  4:171739/1‑138 (MQ=255)
                         ttGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATta                                                                                                                               <  3:171739/138‑1 (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                                                                                                                             >  1:685253/1‑114 (MQ=38)
                                              cttcATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATTCAAAACAACGTGGTTTTGattattc                                                                                                                             <  2:685253/120‑7 (MQ=38)
                                              cttcATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCg                                                                                               >  3:294857/1‑149 (MQ=255)
                                               ttcATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCgg                                                                                              >  1:75501/1‑149 (MQ=255)
                                                      ttCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTaa                                                                                                      >  4:304932/1‑134 (MQ=255)
                                                      ttCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTaa                                                                                                      <  3:304932/134‑1 (MQ=255)
                                                         tGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAg                                                                                    >  4:549306/1‑149 (MQ=255)
                                                            aCGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCgg                                                                                 >  3:419855/1‑149 (MQ=255)
                                                            aCGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCgg                                                                                 >  2:264920/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                                                                                                                 >  4:258174/1‑102 (MQ=255)
                                                                           cGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTta                                                                                                                 <  3:258174/102‑1 (MQ=255)
                                                                              tACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGCGGCAATGTAAAt                                                               <  1:223908/149‑1 (MQ=255)
                                                                                 aaGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGa                                                            >  2:616457/1‑149 (MQ=255)
                                                                                 aaGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGa                                                            >  2:94574/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)
                                                                                                                                                   |                                                                                                                                             
GAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAA‑TCAAAACCACGTTGTTTTGAATTATATATCACGTCTTATAACAAAGTAATGTACCGGTTGTCTGAAGCGGTATGGTGGCAATGTAAATCGAAATCATGTTCACTTTGTATCATGCCGCTTTATTAAATGAAAAGGGAATGATGTGTTGT  >  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: ‑