Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,895,527 G→C 15.0% A361G (GCG→GGG)  adeP ← adenine permease, high affinity, adenine:H+ symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,895,5270GC15.0% 142.9 / 2.5 69A361G (GCG→GGG) adePadenine permease, high affinity, adenine:H+ symporter
Reads supporting (aligned to +/- strand):  ref base G (28/28);  new base C (7/5);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 7.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGCGGGAATTTCCCCTTCTCATCCGC  >  NC_000913/3895293‑3895765
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       
aTCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCATGCACGTTTACCCGTGCCAGACTGGAAATCATCAGACCGCCACCGGAAATAACCGCAc                                                                                                                                                                                                                                                     >  2:833942/1‑238 (MQ=255)
     ccgagtAAACAGCACCCAAACGATTATTACGCCCGGGCTTAGATCCCGCAGCCGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTg                                                                                                                                                                                                                                              <  1:833942/234‑1 (MQ=255)
     cAGGATAAACAGCAGCGCAACGCTGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCACCAGCACGCCAACGTAAATCAGCGCACCAGCTg                                                                                                                                                                                                                                              <  2:750701/239‑1 (MQ=255)
          tAAACAGCATCGCAACGATGATTACGCACGGGCCTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATTCCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGc                                                                                                                                                                                                                                        <  2:598050/238‑1 (MQ=255)
          aaaaCAGCAGCCCAACGATGATTACGCCCGGGCTTAGAGCACGCAGACGTCCGGTACCGATTTTCATCACGGAGTAGGAAATAAAGCCCAGCTCGATACCCTCGGTAATTGAGAATCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACGGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCg                                                                                                                                                                                                                                       <  1:280126/238‑1 (MQ=255)
                gcaacccAACGATTATTACGCCCGGGCTGAGATCCCCCAGACGTCCGGTTCCGATTTTCATCCCGCAGGAGGAGATAAAGCCCAGCGGGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCt                                                                                                                                                                                                                                 <  1:501233/233‑1 (MQ=255)
                     gcAACGATGATTACGCACGGGCTTAGATCACGGAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCa                                                                                                                                                                                                                            <  2:277690/238‑1 (MQ=255)
                              atTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATTACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCGGAACGGCATCATCACGGCGGCAATAAACCCACGAACAGGTTCCATAAAAACCTGACAAGTCACGCGTGCCATACTGGAGGTCATCCACAAACCAAAGTAAACAAGCGAACCAGCCGCGCCGAAACCGGGCACCaacaccc                                                                                                                                                                                                                  >  1:758576/1‑232 (MQ=255)
                                    cacgggattagagtccggagcgttCCGGTGCCTATTGGCATCACGCGGGAGGAAAGAAAGCCCAGGGCGATCCCTGCGGAAATCGAGAAGCTGACGGGCAACATGACGGCGGTAATAAACGCAGAAACAGATTCAGGAAGAACCTGCGAGTTCACGCGTGCCAGGCTGGAAGTCATCAGCCCGCAAACGTAATTGAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAgcgg                                                                                                                                                                                                           <  1:785973/216‑1 (MQ=255)
                                    cACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACCCGTGCCAGGCTGGAAGTCATCAGCACCCCAACGTAAACAAACGCACCAGCTGCAGCGTAGCCTGGGCACCACCCCGCCAgcg                                                                                                                                                                                                            >  2:379072/1‑238 (MQ=255)
                                                   cGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTCAAATCCCGCCAGATCACGCGTGCCACACTGGAAGACACAACGACGCCCACCGAAATCAAGCCACCAACTCCAGCGAAGACCGGCACCCACCCCCCCCCCGGCGCCacacaaaaaaa                                                                                                                                                                                           >  1:567746/1‑229 (MQ=255)
                                                   cGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCCGGCACCATCCCCGCCACCGGCGACCGAAAGATAAc                                                                                                                                                                                           >  2:177540/1‑239 (MQ=255)
                                                   cGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTAAGTAAGATCCTGCCAGTTCACCCGTGCCAAACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTTCAGCGTACCCTGGCCCCACCCCCCCCAccggggccgagaagagaac                                                                                                                                                                                            >  2:590289/1‑224 (MQ=255)
                                                   cGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGAATCAGTAAGATCCTTCCCGTTTAACCGCGCCAGAATTGAAGGCCTACGACCGCCAACCGAAATCAGCGCACCAGCTGAAGCGCAACCTGGCACCAACccccccccgggggaaacaaaaaaaaa                                                                                                                                                                                            >  1:267674/1‑219 (MQ=255)
                                                   cGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAAATTCAGTAAGATCCTGCCCATTCACGCGTGCCAAACTGGAAAGTATCAGCACGCCCACCCAAATCAACCCCCCAACTCCAGCGTAGCCTGGCACCCACCCCCCCAgcggcgccgaaaaaaaaac                                                                                                                                                                                            >  1:756073/1‑226 (MQ=255)
                                                   cGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGGAATAAACGCCGGAAACGATTCAATAAGATCCTGCCAGTTCACCCGTGCCAGGCTGGAAGTCATCAGAACCCCAACGCAAATCAGCCGCCCAGCGTCAAGGGTACCTGGCACCAACccccccccccgcggcacgaagcaaac                                                                                                                                                                                            >  2:590348/1‑219 (MQ=255)
                                                   cGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGGAATAAACACCGGAACAGAATCAAGAAAAACCTGCCAGCTCACGCGTGCCCGACTGGAAGTCATACGCACGCCAACCGAAATAAGGCCACCCACTGCAACGAACCCTGGCACCAACCCCCCCCgcggcggacaaaagagaaa                                                                                                                                                                                            >  1:691727/1‑226 (MQ=255)
                                                         cGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGGACAGATTCAATAAGATCCTGCCAGTTCACCCGTGCCAAACTGGAAGGCATCAACACCCCAACGTAAATCAAGCCACCCCACGCAGCGCACCCTGGCCCCCCCCCCCCCCCCGGCGCCAGAAAAATAAcacacca                                                                                                                                                                                     >  1:771979/1‑233 (MQ=255)
                                                                    gATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAAAGCCCGAACAGATTCAGTAAAATCCTGCCAGGTCACGCGTGCCAGACTGGAAGTCCTCAGCACGCCAACGCAAATCAGCGCACCCCCTGCAGCCGAGCCTGGCCCCATCCCCGCCACCCGCGACCAAAAGAAAAAcacaagaaaaaacaaac                                                                                                                                                                          >  1:141818/1‑224 (MQ=255)
                                                                     aTTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAAAGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACCCGTGCCAGACTGGAAGTCATCCGCACGCCAACCTAAATCAGCGCACCAGCTGCCAGCGAGCCTTGCACCATCCCCGCCACCGGCGGCAGAAAGGTaaaccaccgggcaagcaaacc                                                                                                                                                                         >  2:27336/1‑220 (MQ=255)
                                                                               cGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACCCGTGCCAGACTGGAAGTCAACAGCACGCCCACGTAAAACAGCGCACCAACTGCCGCGTAGCCCGGCCCCACCCCCGCCACCGGCGGCCGAAAGATAAACACCAGGAACAGCAGACcaacaaac                                                                                                                                                                  >  2:830343/1‑234 (MQ=255)
                                                                                  aGTAGGAGAAAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCCCGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGaa                                                                                                                                                                                                  <  2:548609/201‑1 (MQ=255)
                                                                                    tAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGGAATAAACGCAGGGACAAAATCAGTAAGATCCTGCCAGCTCACACGGGCCAGACTGGGAGTCATCCACCAGCCAACGTAACTCAAACCACCCACTGCGCGGTAAACCGGGCCACTCCCCCCCAACGGCGAACGAAAAATAACCACCAGGAAAAGCAAAcacaaaaacaacaccgc                                                                                                                                                          >  1:393705/1‑223 (MQ=255)
                                                                                     aGGAGTAAAAGCCTAGCGCGATACCTTGGGTAATGGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAAAAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGt                                                                                                                                                         <  2:267674/238‑1 (MQ=255)
                                                                                            aaaGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCCCCATCCCCCCCAGCGGCGCCCGAAAGATAACCAGCAGGAACAGCAGACCACAACACACTGCCGCCAAACCgg                                                                                                                                                 >  2:221780/1‑239 (MQ=255)
                                                                                              aGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAAATCCTGCCAGTTCACGCGGGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCCCCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAAAAAACACTGCCGTCAGACCggg                                                                                                                                                >  2:580011/1‑237 (MQ=255)
                                                                                              aGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAAAGCCGGAACAGATTCAGTAAGATACTGCCCGTTCACCCGTGCCAGACTGGAAGGCATCAGCACGCCCACGCAAAACAGCCCACCCCGCTCCACCTAGCCCGGCACCATCCCCCCCCGCGGCGCACAAAAGAAAAACAGCCGGaacaacaaaacaacaaaacaccgccccaaaccggga                                                                                                                                                >  1:586843/1‑219 (MQ=255)
                                                                                                 ccgccgcattcccTTTGGTAATGGAGAACTGAACGGCATCCCCCGCGGGGTAAAAAAACCCAGAAAAAGATTCAGAAGATCCCGGCATGTTCCCGCGGCCAAGATGGGAAGCTCAAACCCCCCCAACTAAAATCAGCCCCCCGGCTGAAGGTTAGCCTGCCCCCACCCCCGCCAGCGGCGCCAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGa                                                                                                                                            <  1:22414/228‑1 (MQ=255)
                                                                                                                  gaaaTGGAAAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTAAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGc                                                                                                                            <  2:771979/236‑1 (MQ=255)
                                                                                                                      cccaaaaGCTGACCGGCAACATCCGGGGGTTAATAAACGCCGGAAAAGATTCAGTAAGACCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGaa                                                                                                                       <  2:25483/234‑1 (MQ=255)
                                                                                                                        ggaaatcgtaaCGGCCTCACCACGCGGTTAATAAACCCGGGACAAAATTCGGTACGACCCTGCCAGTCCAGGCGTGCCAGGCGGGAAGTTACCAGCACCCCAACGAAAAACGGGCCCCAAGCTGAAGCGTGGCCTGCCACCACCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAaga                                                                                                                     <  1:379072/230‑1 (MQ=255)
                                                                                                                        gagaAGCTGAACGGCACCATCACGGCGGTAATAACCGCCGGAACAGATTCAGAAAGATCCGGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCCCGCCAACGTAAATCAGCGCCCCAGCTGCAGCGTAGCCGGGCACCACCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAaga                                                                                                                     <  2:586843/239‑1 (MQ=255)
                                                                                                                            aGCGGAACGGCATCATCAGGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTCCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGa                                                                                                                  <  2:141818/238‑1 (MQ=255)
                                                                                                                             tcTGAACGGCAACACCACGGGGTTAAAAAACCCCGGAAAAGCGTCAGAAAGCTCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCCCGCCAACGGAAATCAGCGCACCAGCTGAAGCGTGGCCGGGCCCCACCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACt                                                                                                                <  1:177540/238‑1 (MQ=255)
                                                                                                                              cgtaaccgccaacctcccggcgaGAATAAACCCCGGAACAGATTCAGAAAGACCCCGCCAGTCCACGCGTGCCAGACTGAAAGCCACCACCACGCCACCGTAAATCGGCGACCCAGCTGCAGCGTGGCCGGGCCCAATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACt                                                                                                                <  1:590289/222‑1 (MQ=255)
                                                                                                                                tgaacggcctcccacggggggTAAAAAACGGCGAAACAGATTCAAAAAACCCCGGCCGTTCCACGCGTGCAAGACGGGAAGTCATCCGCCCCCCAACGTAAAACGGCGCACCAGCGGCAGCGTGGCCGGGCCCCACCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTc                                                                                                               <  1:27336/219‑1 (MQ=255)
                                                                                                                                                    gTAAAAAACCCCGGAACAGTTTCAGGAAGACCCTGCCAGTCCACGGGTGCCAGACGGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGgaag                                                                                           <  2:756073/238‑1 (MQ=255)
                                                                                                                                                      aataaACGCCGGAACAAATTAAGAAAGCCCCCACCAGTCCCCGCGTGCCAGGCTGGAAGTCCCCAGCCCCCCAAGGGAAACCACGCCACCAGCGGCAGCTTGGCCGGGCCCCACCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGaagaag                                                                                        <  1:221780/239‑1 (MQ=255)
                                                                                                                                                          acgcccGGAACAGATTCAGTAAGATCCGGCAAGTCCACGCGTGCCAGACTGGAAGTCATCAGCCCGCAAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCACCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTAcc                                                                                    <  2:23932/235‑1 (MQ=255)
                                                                                                                                                              ccGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACAGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGGTa                                                                                >  2:145355/1‑239 (MQ=255)
                                                                                                                                                               cGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCCGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATTTAAGCCCTAACGTGAGAAGTACCGATaa                                                                               >  2:258089/1‑239 (MQ=255)
                                                                                                                                                                       aTTCAGTAAGCTCCTGCAAGTTCACGCGTGCCAGCCTGGAAGTCCTCAGCACCCCAACGTAAATCAGCGCCCCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATaa                                                                                                                                                                                            <  1:105009/123‑1 (MQ=255)
                                                                                                                                                                       aTTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCCACGTAAATCAGCGCACCAGCTGCAGCGCAGCCTGGCCCCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAAACACCACGCCCGTCAAACCGGTACACACACCACCCGCAGCGCCCCACAATGACCCACTATAACCCGTCACGgaagaaacacacaaaaacaaaaag                                                                       >  1:381462/1‑221 (MQ=255)
                                                                                                                                                                        tccAGTAAGATCCGGCCAGTCCACGCTGGCCAGACTGGAAGTCAACAGCCCGCAAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGaa                                                                          <  2:289807/233‑1 (MQ=255)
                                                                                                                                                                           agagaagtCCTGCCATTTCCCCGGCGCCAACCTGGACGTGAGCAGCGCGCAAACGAAAAGCAGCGCACCGGCTGCGGCGTAGCATGCCACCATCCCCGACAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCGCTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTCCCGATAAACGAACCGGTCa                                                                   <  1:139392/232‑1 (MQ=255)
                                                                                                                                                                            gTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCAc                                                                  >  2:28921/1‑239 (MQ=255)
                                                                                                                                                                              aaGCTCCTGCCAGTTCACGGGGGCCAGCCTGGAAGTCATCAGCACGCCAACGTAAACGAGCGCCCCAGCTGCAGCGAAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCAc                                                                  <  2:784166/237‑1 (MQ=255)
                                                                                                                                                                               aaaacgtcccccattCACGCGTGGCAGACAGGAAAGCCACACCCCCGCAACCGTAAATCAGCCCACCACGGGGAGCGTACCCGGGCCCCACCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCCCTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCAc                                                                  <  1:618402/225‑1 (MQ=255)
                                                                                                                                                                                ataccccgcaaggtCACGCGTGCCAGACGGGAAGCCATCGACAGCCCAACGTAAACCAGCACTACAGCTGCAGCGAAGCCTGGCACCATCCCCGACAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGa                                                               <  2:162623/226‑1 (MQ=255)
                                                                                                                                                                                 aTCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGaa                                                                                                                                                                                  >  2:105009/1‑123 (MQ=255)
                                                                                                                                                                                 aTCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGgaag                                                                                           <  2:781965/209‑1 (MQ=255)
                                                                                                                                                                                 aTCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCAATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGCAGAACCCGTAAACGAACCGGTCACGGAag                                                             >  2:120465/1‑239 (MQ=255)
                                                                                                                                                                                  tCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGCTAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAAAGAACTCAAAATAAGCCGTAACGGAAGAAGTACCGCTAAAGGAACCGGTCACGGagaa                                                            >  2:187923/1‑236 (MQ=255)
                                                                                                                                                                                   ccgcccAGTTCACGCGTGCCAGACGGGAAGTCCTCAGCCCGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCCCCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAagag                                                           <  2:503903/235‑1 (MQ=255)
                                                                                                                                                                                     tGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGACCCGGCCACGGAAGAGAt                                                         >  2:344723/1‑239 (MQ=255)
                                                                                                                                                                                     tGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAAAGAt                                                         >  2:175317/1‑239 (MQ=255)
                                                                                                                                                                                     tGCCAGTTAACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCACCGCCAGCGGCGACATAAAGATAACCAGCAGGAACAGCAGCCCAACAACCCCTGCCGTCAGACAGGTAACACCGCCAACCTGTAAGCCTGCAGAGGCCTCAATATAAACCGCAACGGAAGAAGCAACGATAAACGAACCGGGCACGGaaacaca                                                         >  1:698140/1‑234 (MQ=255)
                                                                                                                                                                                       ccAGTTCCCGCGGGCCAGACTGGAAGTCATGAGCACGCCAACGTAAACCAGCGCACCAGGTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATAc                                                       <  2:485427/239‑1 (MQ=255)
                                                                                                                                                                                            ttCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCg                                                                                                                                         >  1:140170/1‑153 (MQ=255)
                                                                                                                                                                                            ttCACGCGTGCAAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCg                                                                                                                                         <  2:140170/153‑1 (MQ=255)
                                                                                                                                                                                             tcccggttGCCAGACTGGAAGTCACCGGCACGCAACCGTAAATCAGGGCACCAGCTGCAGCGTAGCCTGGCACCACCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTGTCg                                                  <  2:211445/232‑1 (MQ=255)
                                                                                                                                                                                             ccgcgcgTGCCAGTCTGGATGTCACCAGCCCCAAAACGTAAATCAGCGCACCAGCTGCAGGGAAGCGTGGCACAACCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTg                                                     <  2:393705/233‑1 (MQ=255)
                                                                                                                                                                                                      ccAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCCCGGAAGGGAAACTGTCgcacaacag                                          >  2:337856/1‑230 (MQ=255)
                                                                                                                                                                                                                      atcaGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACCACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGCGGACTCAATATAAACCGTACCGGAAGAAATACCGAAAAACAAACCGGTCACGGAACAGACACAGTAAACAAACAgaacacgcctcaagcg                           >  1:395073/1‑224 (MQ=255)
                                                                                                                                                                                                                          gCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGCggg                      >  2:596341/1‑239 (MQ=255)
                                                                                                                                                                                                                           cACGACAACGTAAATGAGCCCCCCAGCAGCAGGCTAGCCGGCACCCACCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGCGGGa                     <  2:12028/239‑1 (MQ=255)
                                                                                                                                                                                                                                cAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACCGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCACTCTAAGCCGTAACGGAAGAAGTACCGCTAAACGAACCGGTCACGGACGAGATACTGTCGCCATACAACGCCTGCTTCATCCCCGGGAAttt                 >  2:96015/1‑238 (MQ=255)
                                                                                                                                                                                                                                        aTCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGCTAAACGAACCGGTCACGGAAGAGATACTGTCGCCATACAGCGCCTGCTTCATGCGCGGGAATTTCCCCTtctc        >  2:885110/1‑239 (MQ=255)
                                                                                                                                                                                                                                        aTCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTGCCGACATACAGCGCCTCCTTCATGCGCGGGAATTTCCACTcct         >  2:689137/1‑238 (MQ=255)
                                                                                                                                                                                                                                            gcgcACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGCTAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCCGCTTCAAGCGCGGGAAATTCCCCTTCCCaac     >  1:272930/1‑236 (MQ=255)
                                                                                                                                                                                                                                             cgcACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGCGGGAATTTCCCCTTCTCATcc    >  2:514246/1‑238 (MQ=255)
                                                                                                                                                                                                                                             cgcACAAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCATCAGGAACAGCAGACCACCAACCACTGCCGACAGACCGGTACGACCGCCAACCGCAACGCCGGCAGAGGCCCCACTATAAGCCGGACAGGAAGAACTCCCGACACACGACCCGGCCACGGAAGGAACCCTTGCGCCATACAACGCCCGCTCCCTTCGCGGGCATTTCCCCCTCCCAAcc    >  1:772742/1‑238 (MQ=255)
                                                                                                                                                                                                                                              gcACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGCAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGCCACGGCAGAGATACTGTAGACAAACAACGCCTGCTTCATCGCGGGGAATTTCCCCTCCCCACCCGc  >  1:688555/1‑239 (MQ=255)
                                                                                                                                                                                                                                              gcACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAAAGATACTGTCGACATACAGCGCCTGCTTCATCCCCGGGAATTTCCCCTTCCCATcc  >  1:48077/1‑237 (MQ=255)
                                                                                                                                                                                                                                              gcACCAGTGCAAGCGTAGCCTGGCACCACCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGCGGGAATTTCCCCTTCTCATCCGc  <  2:641166/239‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       
ATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGCGGGAATTTCCCCTTCTCATCCGC  >  NC_000913/3895293‑3895765

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