Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,998,056 C→G 7.9% A142P (GCG→CCG)  yecS ← ABC family putative inner membrane permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,998,0560CG 7.9% 247.7 / 12.2 77A142P (GCG→CCG) yecSABC family putative inner membrane permease
Reads supporting (aligned to +/- strand):  ref base C (39/31);  new base G (2/4);  total (42/35)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.73e-01

GCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAGCGTTTCGGCAGCATAGGCGGCAGTGTTCAACGACAGACCAAT  >  NC_000913/1997908‑1998203
                                                                                                                                                    |                                                                                                                                                   
gccAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACgcgc                                                                                                                                                     >  4:837171/1‑149 (MQ=255)
  cAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGcgccg                                                                                                                                                   >  1:543812/1‑149 (MQ=255)
        cATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCtgtg                                                                                                                                             <  2:7589/149‑1 (MQ=255)
          tGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGc                                                                                                                                           >  1:105680/1‑149 (MQ=255)
             tGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACAGGCGGCGCGTGTGGcgtt                                                                                                                                        <  2:820435/149‑4 (MQ=255)
             tGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACAGGCGGCGCGTGTGGcgtt                                                                                                                                        <  3:312943/149‑4 (MQ=255)
             tGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCaaa                                                                                                                                        >  3:54995/1‑149 (MQ=255)
               aaCACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACAGGCGGCGCGTGTGGcgttgc                                                                                                                                      <  2:171954/149‑6 (MQ=255)
                acacCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATc                                                                                                                                     <  1:72996/149‑1 (MQ=255)
                    cTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATcgcgc                                                                                                                                 >  2:821383/1‑149 (MQ=255)
                     tCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATcgcgcg                                                                                                                                >  3:435664/1‑149 (MQ=255)
                     tCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATcgcgcg                                                                                                                                >  4:43011/1‑149 (MQ=255)
                       cAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGAc                                                                                                                              <  3:439111/149‑1 (MQ=255)
                           gTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCAt                                                                                                                          >  1:281098/1‑149 (MQ=255)
                               tcgAGGTAATCAGCCGCGCCTGACGGTACAGCTCCGGCACCAGGACTGTCGCGGAGAGCGACGTATCCTTTACCAGCCTAATAAAGCAGTTCGACAGCGGGGGCAACGCCACACGCGCCAGCTGTGGCAAAATCGCGCGACGCaagctc                                                                                                                      <  4:105938/148‑6 (MQ=255)
                                cgAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCt                                                                                                                     <  1:328253/149‑1 (MQ=255)
                                     tAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCAc                                                                                                                >  4:222297/1‑149 (MQ=255)
                                      aaTCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACg                                                                                                               >  1:42891/1‑149 (MQ=255)
                                      aaTCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACg                                                                                                               >  1:693028/1‑149 (MQ=255)
                                       aTCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACgg                                                                                                              >  4:467922/1‑149 (MQ=255)
                                       aTCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACgg                                                                                                              <  3:148139/149‑1 (MQ=255)
                                       aTCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACgg                                                                                                              >  1:440987/1‑149 (MQ=255)
                                       aTCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACgg                                                                                                              >  1:526767/1‑149 (MQ=255)
                                        tCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGTCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGc                                                                                                             >  2:741310/1‑149 (MQ=255)
                                          aaCTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGt                                                                                                           <  1:357914/149‑1 (MQ=255)
                                          aaCTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGt                                                                                                           <  4:34471/149‑1 (MQ=255)
                                                    gACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATa                                                                                                 <  4:54995/149‑1 (MQ=255)
                                                       ggAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACATCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTg                                                                                              >  4:722011/1‑149 (MQ=255)
                                                       ggAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAAAGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTg                                                                                              >  4:613391/1‑149 (MQ=255)
                                                         aaCAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGc                                                                                            >  2:99569/1‑149 (MQ=255)
                                                         aaCAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGc                                                                                            >  1:824172/1‑149 (MQ=255)
                                                         aaCAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGc                                                                                            <  2:266980/149‑1 (MQ=255)
                                                          aCAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCt                                                                                           >  2:637400/1‑149 (MQ=255)
                                                          aCAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCt                                                                                           >  3:21143/1‑149 (MQ=255)
                                                            aGCTCCGGCACCTGGATTGTCGCGGCCAGCGGCGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGc                                                                                         >  4:541635/1‑149 (MQ=255)
                                                            aGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGTCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGc                                                                                         >  2:727496/1‑149 (MQ=255)
                                                            aGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGc                                                                                         >  4:475273/1‑149 (MQ=255)
                                                            aGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGc                                                                                         >  3:398060/1‑149 (MQ=255)
                                                             gCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGcc                                                                                        >  4:618527/1‑149 (MQ=255)
                                                                    cACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTccc                                                                                 <  1:384817/149‑1 (MQ=255)
                                                                    cACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTccc                                                                                 <  3:83774/149‑1 (MQ=255)
                                                                     aCCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCa                                                                                <  1:426845/149‑1 (MQ=255)
                                                                     aCCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCa                                                                                <  1:728951/149‑1 (MQ=255)
                                                                        tGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTg                                                                             >  3:258189/1‑149 (MQ=255)
                                                                         ggATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGa                                                                            <  1:464484/149‑1 (MQ=255)
                                                                           aTTGTCGCGGCCACCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGAcc                                                                          <  4:559349/149‑1 (MQ=255)
                                                                                  cgGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCa                                                                   >  3:301233/1‑149 (MQ=255)
                                                                                   ggCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTGATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCaa                                                                  <  2:281098/149‑1 (MQ=255)
                                                                                        gCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATagaa                                                             <  4:129946/149‑1 (MQ=255)
                                                                                        gCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATagaa                                                             <  3:222297/149‑1 (MQ=255)
                                                                                             gTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATagaa                                                             <  2:155393/144‑1 (MQ=255)
                                                                                             gTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATagaa                                                             >  1:155393/1‑144 (MQ=255)
                                                                                              tATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATg                                                       <  2:373651/149‑1 (MQ=255)
                                                                                                  ttttACCAGGCTGATAAAGTTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTTGCTGCCGCTACCCACTGACCTTTATCAATAGAAGAAATggcgg                                                   >  2:498216/1‑149 (MQ=255)
                                                                                                  ttttACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATggcgg                                                   >  1:133109/1‑149 (MQ=255)
                                                                                                  ttttACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATggcgg                                                   >  4:113602/1‑149 (MQ=255)
                                                                                                        cAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCa                                             >  4:189905/1‑149 (MQ=255)
                                                                                                         aGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAg                                            >  4:738346/1‑149 (MQ=255)
                                                                                                          ggCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATa                                                                <  3:555203/128‑1 (MQ=255)
                                                                                                          ggCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATa                                                                >  4:555203/1‑128 (MQ=255)
                                                                                                             taataaAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAGCGtt                                        >  4:470613/1‑149 (MQ=255)
                                                                                                               ataaAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAGCGTTTc                                      <  4:435664/149‑1 (MQ=255)
                                                                                                                  aaGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAGCGTTTCGgc                                   <  1:59549/149‑1 (MQ=255)
                                                                                                                   aGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAGCGTTTCGgca                                  >  2:119496/1‑149 (MQ=255)
                                                                                                                     cTGTTCGAAAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAGCGTTTCGgcagc                                <  4:783516/149‑1 (MQ=255)
                                                                                                                         tCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAGCGTTTCGGCAGCATAg                            >  2:725726/1‑149 (MQ=255)
                                                                                                                           gACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAGCGTTTCGGCAGCATAggc                          >  3:263732/1‑149 (MQ=255)
                                                                                                                                 ggcggcAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAGCGTTTCGGCAGCATAGGCGGCAgt                    >  3:523990/1‑149 (MQ=255)
                                                                                                                                       aaCGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAGCGTTTCGGCAGCATAGGCGGCAGTGTTCaa              <  1:498954/149‑1 (MQ=255)
                                                                                                                                        aCGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAGCGTTTCGGCAGCATAGGCGGCAGTGTTCAAc             <  3:113602/149‑1 (MQ=255)
                                                                                                                                          gCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAGCGTTTCGGCAGCATAGGCGGCAGTGTTCAACGa           >  3:154642/1‑149 (MQ=255)
                                                                                                                                           ccACAGGCGGCGCGTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAGCGTTTCGGCAGCATAggcgg                        >  4:435301/1‑135 (MQ=255)
                                                                                                                                           ccACAGGCGGCGCGTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAGCGTTTCGGCAGCATAggcgg                        <  3:435301/135‑1 (MQ=255)
                                                                                                                                           ccACAGGCGGCGCGTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAGCGTTTCGGCAGCATAGGCGGCAGTGTTCAACGAc          >  4:692079/1‑149 (MQ=255)
                                                                                                                                           ccACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCa                                             <  2:297197/114‑1 (MQ=255)
                                                                                                                                           ccACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCa                                             >  1:297197/1‑114 (MQ=255)
                                                                                                                                            cacaCGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAGCGTTTCGGCAGCATAGGCGGCAGTGTTCAACGACa         <  2:824172/149‑1 (MQ=255)
                                                                                                                                            cacaCGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAGCGTTTCGGCAGCATAGGCGGCAGTGTTCAACGACa         <  1:99569/149‑1 (MQ=255)
                                                                                                                                               acGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAGCGTTTCGGCAGCATAGGCGGCAGTGTTCAACGACAGAc      <  3:58662/149‑1 (MQ=255)
                                                                                                                                                  cgccgcCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAGCGTTTCGGCAGCATAGGCGGCAGTGTTCAACGACAGACCaa   <  4:516601/149‑1 (MQ=255)
                                                                                                                                                   gccgccTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAGCGTTTCGGCAGCATAGGCGGCAGTGTTCAACGACAGACCAAt  <  1:41562/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAGCGTTTCGGCAGCATAGGCGGCAGTGTTCAACGACAGACCAAT  >  NC_000913/1997908‑1998203

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑