Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,998,052 C→G 7.7% R143P (CGT→CCT)  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,0520CG 7.7% 266.0 / 12.0 78R143P (CGT→CCT) yecSABC family putative inner membrane permease
Reads supporting (aligned to +/- strand):  ref base C (42/30);  new base G (2/4);  total (44/34)
Fisher's exact test for biased strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

AGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAGCGTTTCGGCAGCATAGGCGGCAGTGTTCAACGACAGAC  >  NC_000913/1997904‑1998199
                                                                                                                                                    |                                                                                                                                                   
aGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCacac                                                                                                                                                     >  2:266515/1‑149 (MQ=255)
  ccgccAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACAcgc                                                                                                                                                   <  2:367161/149‑1 (MQ=255)
   cgccAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACAcgcg                                                                                                                                                  >  2:357914/1‑149 (MQ=255)
    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                                                                                                                            >  4:43011/1‑149 (MQ=255)
                         tCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATcgcgcg                                                                                                                            >  3:435664/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:693028/1‑149 (MQ=255)
                                          aaTCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACg                                                                                                           >  1:42891/1‑149 (MQ=255)
                                           aTCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACgg                                                                                                          <  3:148139/149‑1 (MQ=255)
                                           aTCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACgg                                                                                                          >  1:526767/1‑149 (MQ=255)
                                           aTCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACgg                                                                                                          >  4:467922/1‑149 (MQ=255)
                                           aTCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACgg                                                                                                          >  1:440987/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                                                                                        <  2:266980/149‑1 (MQ=255)
                                                             aaCAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGc                                                                                        >  1:824172/1‑149 (MQ=255)
                                                              aCAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCt                                                                                       >  3:21143/1‑149 (MQ=255)
                                                              aCAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCt                                                                                       >  2:637400/1‑149 (MQ=255)
                                                                aGCTCCGGCACCTGGATTGTCGCGGCCAGCGGCGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGc                                                                                     >  4:541635/1‑149 (MQ=255)
                                                                aGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGTCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGc                                                                                     >  2:727496/1‑149 (MQ=255)
                                                                aGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGc                                                                                     >  3:398060/1‑149 (MQ=255)
                                                                aGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGc                                                                                     >  4:475273/1‑149 (MQ=255)
                                                                 gCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGcc                                                                                    >  4:618527/1‑149 (MQ=255)
                                                                        cACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTccc                                                                             <  3:83774/149‑1 (MQ=255)
                                                                        cACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTccc                                                                             <  1:384817/149‑1 (MQ=255)
                                                                         aCCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCa                                                                            <  1:728951/149‑1 (MQ=255)
                                                                         aCCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCa                                                                            <  1:426845/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                                                         <  3:222297/149‑1 (MQ=255)
                                                                                            gCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATagaa                                                         <  4:129946/149‑1 (MQ=255)
                                                                                                 gTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATagaa                                                         >  1:155393/1‑144 (MQ=255)
                                                                                                 gTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATagaa                                                         <  2:155393/144‑1 (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                                         >  1:297197/1‑114 (MQ=255)
                                                                                                                                               ccACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCa                                         <  2:297197/114‑1 (MQ=255)
                                                                                                                                                cacaCGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAGCGTTTCGGCAGCATAGGCGGCAGTGTTCAACGACa     <  1:99569/149‑1 (MQ=255)
                                                                                                                                                cacaCGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAGCGTTTCGGCAGCATAGGCGGCAGTGTTCAACGACa     <  2:824172/149‑1 (MQ=255)
                                                                                                                                                   acGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAGCGTTTCGGCAGCATAGGCGGCAGTGTTCAACGACAGAc  <  3:58662/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAGCGTTTCGGCAGCATAGGCGGCAGTGTTCAACGACAGAC  >  NC_000913/1997904‑1998199

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