Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,998,060 C→G 7.7% Q140H (CAG→CAC)  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,0600CG 7.7% 258.3 / 12.1 78Q140H (CAG→CAC) yecSABC family putative inner membrane permease
Reads supporting (aligned to +/- strand):  ref base C (40/32);  new base G (2/4);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

CATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAGCGTTTCGGCAGCATAGGCGGCAGTGTTCAACGACAGACCAATCATT  >  NC_000913/1997916‑1998207
                                                                                                                                                |                                                                                                                                                   
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                                                                                                                  >  1:440987/1‑149 (MQ=255)
                               aTCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACgg                                                                                                                  >  4:467922/1‑149 (MQ=255)
                               aTCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACgg                                                                                                                  >  1:526767/1‑149 (MQ=255)
                               aTCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACgg                                                                                                                  <  3:148139/149‑1 (MQ=255)
                                tCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGTCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGc                                                                                                                 >  2:741310/1‑149 (MQ=255)
                                  aaCTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGt                                                                                                               <  4:34471/149‑1 (MQ=255)
                                  aaCTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGt                                                                                                               <  1:357914/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:266980/149‑1 (MQ=255)
                                                 aaCAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGc                                                                                                >  2:99569/1‑149 (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                                                                 <  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                                                       >  4:113602/1‑149 (MQ=255)
                                                                                          ttttACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATggcgg                                                       >  1:133109/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)
                                                                                                                                               ccTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAGCGTTTCGGCAGCATAGGCGGCAGTGTTCAACGACAGACCAATCAtt  <  4:523990/149‑1 (MQ=255)
                                                                                                                                               ccTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAGCGTTTCGGCAGCATAGGCGGCAGTGTTCAACGACAGACCAATCAtt  <  3:680563/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTGTCGCGGCCAGCGACGTATCTTTTACCAGGCTAATAAAGCTGTTCGACAGCGGCGGCAACGCCACACGCGCCGCCTGTGGCAAAATCGCGCGACGCATGGTCTGCCACGGCGTCATACCAATACTGGCTGCCGCTTCCCACTGACCTTTATCAATAGAAGAAATGGCGGCGCGCAGCGTTTCGGCAGCATAGGCGGCAGTGTTCAACGACAGACCAATCATT  >  NC_000913/1997916‑1998207

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: ‑