Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,512,293 T→C H41R (CAC→CGC)  fecD ← ferric citrate ABC transporter membrane subunit FecD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,512,2930TC98.6% 208.7 / ‑4.5 73H41R (CAC→CGC) fecDferric citrate ABC transporter membrane subunit FecD
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (38/34);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-01

CAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGG  >  NC_000913/4512046‑4512538
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     
cacccccGCCAGGCCCCCCGCAAAGGCCAGCGGCGGCAGCACCACCACGGGCAGTGACGGCAAAAGAAGTAGAGCCCCCACGGAGGCCAGGCGGGCGGCATGGTTACCCCCGAGAATACCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTc                                                                                                                                                                                                                                                    <  2:117884/248‑1 (MQ=255)
       cccaggccgccCGCAAAGGCCAGCAGCGGCAGCCCCATCCCGGGCAGTGACGGCAAAAGAAGTAGAGCCCCCCCAGAGGCCAGGCTGGCGGCATGGTTAACCCCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTg                                                                                                                                                                                                                                             <  1:407466/246‑1 (MQ=255)
              cgcccgcAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCCCGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCCAGCGCGGCAGTCGGGACTCCATCAATACATACTAATGCTCGCGTCCGCCCTGcccagcc                                                                                                                                                                                                                                      >  2:345535/1‑246 (MQ=255)
               gcccgcAAAGGCCAGCAGCGGCAGCAACATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCgg                                                                                                                                                                                                                                     >  1:159348/1‑251 (MQ=255)
                          ccAGCAGCGGTAGCACCATCCGGGGCAGTGACGGCAAAAGAAGAGGAGCCCCCACAGAGGCCAGTCGTGCGGCATGGTAACCACCAAAAAGACCCGGTTATGCCGGAGGGGTGCGCACAACCCCCGTTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAgcgcg                                                                                                                                                                                                                          <  2:40986/251‑1 (MQ=255)
                           cagcagCGGCAGCACCATCACGGGCAGTGCGGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACCCCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGcgctgcaattggtcgcgcaaattc                                                                                                                                                                                                                         <  1:209533/251‑24 (MQ=255)
                           cagcagCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACCCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAgcgcgc                                                                                                                                                                                                                         >  1:268333/1‑251 (MQ=255)
                              cagcGGCAGCACCATCACGGGCCGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGCCAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATCAATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGcc                                                                                                                                                                                                                      >  2:329019/1‑251 (MQ=255)
                                 cGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCcacg                                                                                                                                                                                                                   >  1:344933/1‑251 (MQ=255)
                                          cATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGGACTCCATCAAAACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACggggg                                                                                                                                                                                                          >  2:486352/1‑250 (MQ=255)
                                             cACGGGCAGTGACGGCATAAGAAGGAGAGCCCCAACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTGGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCa                                                                                                                                                                                                       <  1:468804/251‑1 (MQ=255)
                                                 ggCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCGGGCGGCATGGTTAACCCCGAGAATATCCGGTGATGCCAGAGGTTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTcc                                                                                                                                                                                                   <  1:448446/251‑1 (MQ=255)
                                                     gTGAGGGCATAAGAAGGAGGCCCCCAAAAGAGGCCGGGCGGGCGGCGTGGTTAACCCCGAGAATACCCGGTGAGGCCAGAGGGTTGCGCCCAACCCCCTGTACCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCata                                                                                                                                                                                               <  2:270762/251‑1 (MQ=255)
                                                       gACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCCCGGGGGTCACCCCCCTATg                                                                                                                                                                                             >  2:274267/1‑251 (MQ=255)
                                                             atagaaagTAGACCCCCCACATAGGCCGGGATGGGGGCTGCGTTAGCCCCGATAATCCCCGGTGATGCCAGAGGGTCTCTCACAACCCCCTTCATCACCCCGCCCGACACGGCGAGGGCTCCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGGCGGTACTCCAGCAATACATAATAATGCTCGCATCCGGCCTTCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGa                                                                                                                                                                                       <  2:417268/246‑1 (MQ=255)
                                                              aaaaaagTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGAt                                                                                                                                                                                      <  1:424741/247‑1 (MQ=255)
                                                                 gaagTAGAGCCCCCACGGAGGCGTGGCTGGCGGCAGGTTTAACCCGGAGAATCCCGGGTGTTCCCAGGGGGTTGCGCACAATCCCTTGTATCAGCCCGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAAc                                                                                                                                                                                   <  2:123861/251‑1 (MQ=255)
                                                                        agCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACACTCCCCTGTATCAGCACGCCCGCCACGGCGGGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCCGTCGGGACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGGGCCCGCCACGGCACGGGGGACACTCCCCTCTGGAGTGATAAAAGCGcac                                                                                                                                                                            >  2:126669/1‑251 (MQ=255)
                                                                            cccACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGTTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGcc                                                                                                                                                                        >  1:295479/1‑251 (MQ=255)
                                                                                 agagGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGcc                                                                                                                                                                    >  1:406016/1‑250 (MQ=255)
                                                                                 agagGCCAGGCTGGCGGCAGGGTTAACACCGAGAATACCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCCCGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCa                                                                                                                                                                   <  1:298132/251‑1 (MQ=255)
                                                                                    gggCCAGGCTGGCGGCAGGGTTAACCCCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAACCCCCTGTATAAGCCCGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAgg                                                                                                                                                                 <  1:486352/250‑1 (MQ=255)
                                                                                    agGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGACCGCCCGCCACGGCGAGGGCAGCACCGACAACCAGTGCCAGCAGCAAGCGCGGCAGTCAGTAATCCATCAATACCTAATAATGCTCGCGTCCAGCCTGCCAGCCGGTCAGCAGCGCGAGCACCGGCAAGGGGATCCCACCCATATGGGGTGAttgaagcacacacccccgcagg                                                                                                                                                                 >  2:469009/1‑230 (MQ=255)
                                                                                         ggggtGGCGGCATGGTTAACACCGAGAAAACCCGGGGATGCCAGAGGGTGGCGCACAACCCCCTGTATCAGCCCGCCCGCCCCGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCAGATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAg                                                                                                                                                            <  2:390378/247‑1 (MQ=255)
                                                                                         aggcaggcGGCATGGTTAACACCGAGACTATCCGGTGATGCCAGAGGGTTGAGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGAGCGGCAGTCGGTACTCCATCAAAACATAATAATGCACGCGTCCGGCCTGCACGTCGGGCAGCAGCGCGCACCACGGAGCGGGGAACACTCCCATCTGGAGTGATTACAgcacaacgcacgcaagggca                                                                                                                                                              >  2:370400/1‑231 (MQ=255)
                                                                                             gtgcggcATTGTTAACACGAAGAAAATCGGGTGACTCAAGGGGGTTGCCCACAATCCCCTGTAACAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCGCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACATCCCGCCAGGGCAAGGGtg                                                                                                                                                        <  2:327593/249‑1 (MQ=255)
                                                                                              ggcggcATGGTTAACACCGAGAATATCCGATGATGCCAGAGGGTTGCGCACAACCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGt                                                                                                                                                         <  1:389477/249‑1 (MQ=255)
                                                                                                        gTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGAAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGCTCACTCCCATATGGAGTGATAACAGCGAACAGCCCGCCAGGGCAAGGGTGATGAAAATa                                                                                                                                               >  1:483339/1‑250 (MQ=255)
                                                                                                           aaCACCGAGAATATACGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGAACAGACCGCCAGGGCAAGGATGATGAAAATTAcc                                                                                                                                            >  1:400026/1‑250 (MQ=255)
                                                                                                               ccGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGAGATGAAAATAACCAgcgc                                                                                                                                       >  1:437736/1‑251 (MQ=255)
                                                                                                                 gagaATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGGGATGAAAATAACCAGCGCa                                                                                                                                      >  2:227934/1‑250 (MQ=255)
                                                                                                                   gaATATCCGGTGATGCCAGAGGGTGGCGCACAACCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAAtt                                                                                                                                   <  1:274267/251‑1 (MQ=255)
                                                                                                                           ggTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATcgt                                                                                                                           >  2:123359/1‑249 (MQ=255)
                                                                                                                            gTGATGCCAGAGGGTTGCGCACAACCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTc                                                                                                                          <  1:236849/251‑1 (MQ=255)
                                                                                                                             ggctgccagggggtggccaaaaacccccTGTTGCAGCACCCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACTGCACGGGGATCACTCCCATATGGAGTGATAACAGCTCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCg                                                                                                                         <  2:196180/227‑1 (MQ=255)
                                                                                                                                    agggggTTGCGCACAACCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCctcctc                                                                                                                  <  1:329019/248‑1 (MQ=255)
                                                                                                                                     gagGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCctcctc                                                                                                                  >  1:397692/1‑250 (MQ=255)
                                                                                                                                      agGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCGAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTcac                                                                                                                >  1:224703/1‑251 (MQ=255)
                                                                                                                                             cgcACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATGACCAGCGCAATTTTCATCCTCGCCGCCTCACAAGCCa                                                                                                          >  1:335163/1‑250 (MQ=255)
                                                                                                                                             cgcACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGcc                                                                                                           >  1:28929/1‑249 (MQ=255)
                                                                                                                                               cacaATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCCAGGGTGATGAAAATAACCAGCGCAATTTTCCTCCTCGCCTCCTCCCAAGCCCGAc                                                                                                       >  2:260891/1‑251 (MQ=255)
                                                                                                                                                acaACCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACa                                                                                                      <  1:345535/251‑1 (MQ=255)
                                                                                                                                                        cTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAgg                                                                                              >  1:310013/1‑251 (MQ=255)
                                                                                                                                                           tATCAGCACGCCAGCCACGGCGAGGGCTGCACCGCCAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCt                                                                                           >  1:377677/1‑251 (MQ=255)
                                                                                                                                                               aGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCACGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCCGAGCAATTTTCCTCCTCGCCTCCTCCCAAGCCACACAAAGCCAGGGCTGcc                                                                                        >  2:189931/1‑250 (MQ=255)
                                                                                                                                                                    gcccgccACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAAACAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATc                                                                                    >  1:113089/1‑249 (MQ=255)
                                                                                                                                                                              gCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACt                                                                        <  1:259968/251‑1 (MQ=255)
                                                                                                                                                                               cGAGGGCTGCACCGATAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTg                                                                       <  1:205257/251‑1 (MQ=255)
                                                                                                                                                                                     cggacccgcaaaccAGTGCCAGCCGCAAGCGCGGCCGTCGGACCTCCATCAATACATAATAAGACCCGCGTCCTCCCTGCCAGTCGGTCAGCACCGCGCGCCCCGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAAGAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGCCAAAGCAAGGGCTGCCAAGCAGCGCCAGCACTGCGCCTg                                                                 <  2:413178/238‑1 (MQ=255)
                                                                                                                                                                                       gCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTg                                                                 <  1:381096/249‑1 (MQ=255)
                                                                                                                                                                                            gACAAACAGTGCCAGCAGCAAGCGCGGCAGGCGGTACTCCATCAATACATAATAAGGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCACCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAg                                                          <  2:159348/251‑1 (MQ=255)
                                                                                                                                                                                             aCAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCAATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAAGGCTGCCAATCAGCGACAGCACTGCGCCTGCGGGc                                                            >  1:270921/1‑248 (MQ=255)
                                                                                                                                                                                              cAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGAt                                                        <  1:493456/251‑1 (MQ=255)
                                                                                                                                                                                                    gTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCGCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCgg                                                  <  1:62078/251‑1 (MQ=255)
                                                                                                                                                                                                         agcagcAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGcc                                                                                                                       >  2:249911/1‑177 (MQ=255)
                                                                                                                                                                                                         agcagcAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGcc                                                                                                                       <  1:249911/177‑1 (MQ=255)
                                                                                                                                                                                                            agcaAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGcc                                          <  1:374735/251‑1 (MQ=255)
                                                                                                                                                                                                             gcaAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGcc                                          >  1:318623/1‑250 (MQ=255)
                                                                                                                                                                                                               aaGCGCGGCAGTCGGTACTCCATCAATACATAATACTGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGCGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGGGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGa                                                                                                        >  2:47917/1‑186 (MQ=255)
                                                                                                                                                                                                               aaGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGa                                                                                                        <  1:47917/186‑1 (MQ=255)
                                                                                                                                                                                                               aaGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAgt                                       >  1:239578/1‑250 (MQ=255)
                                                                                                                                                                                                               aaGCGCGGCAGTCGGTACTCCATCAATACATAATAATGATCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAgt                                       >  1:354773/1‑250 (MQ=255)
                                                                                                                                                                                                                      gCAGTCGGTAATCCATCAATACATACTAATGCTCGCGTCAGGCCAGCCAGTCGGACAGCAGAGCGCGCCACGGCACGGGGATCACTCCAATATGGAGTGATACCAGCGCACAGCCCGCCAGGGCAAGGGTGATGAACATAACCAGCGCAATTTTCATCCTAGCCTCCTCACCATCCCGCAAAAGCAAGGGCTGCCAATCAGCGCCATCACTGCGCCTGCGGGCAGATCTCCGGGGCAGTACAGGGCGCGTg                                >  2:199165/1‑251 (MQ=255)
                                                                                                                                                                                                                           cGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGt                           <  1:496885/251‑1 (MQ=255)
                                                                                                                                                                                                                              tACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGCTGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCTCTGCGCCTGCGGGCAGATCTCCAGGGAAGGCCAGTGCGCGTGCGAGTAc                         >  2:349777/1‑250 (MQ=255)
                                                                                                                                                                                                                              tACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGGGCGCGTGCGAGTAc                         >  1:275679/1‑250 (MQ=255)
                                                                                                                                                                                                                                  ccATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGCGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCt                     <  1:355242/250‑1 (MQ=255)
                                                                                                                                                                                                                                        aTACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTg                                                                                          <  2:394810/175‑1 (MQ=255)
                                                                                                                                                                                                                                        aTACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTg                                                                                          >  1:394810/1‑175 (MQ=255)
                                                                                                                                                                                                                                            ataataATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGCGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAgcagcac          >  2:144380/1‑250 (MQ=255)
                                                                                                                                                                                                                                            ataataAGGCCCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCAt          <  2:215487/251‑1 (MQ=255)
                                                                                                                                                                                                                                              aataatGCTCGCGTCCGGCCTGCCAGGCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATc         <  1:70726/250‑1 (MQ=255)
                                                                                                                                                                                                                                                    gCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTgg  >  1:164347/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     
CAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGG  >  NC_000913/4512046‑4512538

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: