Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 72,605 C→T G108E (GGG→GAG)  thiQ ← thiamine ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091372,6050CT97.8% 275.3 / NA 93G108E (GGG→GAG) thiQthiamine ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (31/59);  minor base G (0/2);  total (32/61)
Fisher's exact test for biased strand distribution p-value = 5.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCCACGCTGTGCGACACCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTACCAGACAACGCGCTAACGCCACTCGCTGTCGCTGACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCC  >  NC_000913/72355‑72855
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                          
ttCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGACTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCACCAGGTAACCGCGCCATTAAATTATCAATCCACATCTGGCGGGCGATAGCGTGCATTTTcc                                                                                                                                                                                                                                                            >  2:513750/1‑251 (MQ=255)
      ccTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGGGCTCCCCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCctgct                                                                                                                                                                                                                                                      <  2:92005/250‑1 (MQ=255)
      acTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCctgct                                                                                                                                                                                                                                                      <  1:272069/250‑1 (MQ=255)
           gtgACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGtgcc                                                                                                                                                                                                                                                 <  1:366056/249‑4 (MQ=255)
                  caccAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCTACGTCAGCATCTCCTGACGTAACGCCGGATCGAGCTCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCa                                                                                                                                                                                                                                          <  2:552869/248‑1 (MQ=255)
                   aTCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAg                                                                                                                                                                                                                                         <  1:422700/251‑1 (MQ=255)
                      aacaacGTCATTTTTGGCTGCTGGCAGCTCGTTCTCCCCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGttt                                                                                                                                                                                                                                      <  2:207465/246‑1 (MQ=255)
                        taacGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTc                                                                                                                                                                                                                                     >  1:140477/4‑250 (MQ=255)
                         accGTCATTTTGTGGTGTTGGCAGCTCGTGCTCACCAACGTCAGCATCGCCTGACGGAACGACGGATCGGGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGTTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCAGCTGGCGGGCGATAGCGTGCATTTGCTCCTGCTGTGCCGCGTTCAGTTTCa                                                                                                                                                                                                                                    <  2:288289/248‑1 (MQ=255)
                                        ggctgGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGGTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAAGCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGccc                                                                                                                                                                                                                    <  1:663638/250‑1 (MQ=255)
                                              cAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGAGCCATTAAATTATCAATCCCCATCTTGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGGTCAGTTTCAATCCCGGATTTAGACCCAGccc                                                                                                                                                                                                              >  2:177329/1‑251 (MQ=255)
                                               aGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCg                                                                                                                                                                                                             >  1:669622/1‑251 (MQ=255)
                                                   cGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATg                                                                                                                                                                                                          <  1:707109/250‑1 (MQ=255)
                                                    gggctccccagcgacaggaTCTCCAGACTTAACGCCGGATCGAGCGCAGAGACCGGTCCACCGAGCAATAAAACCGGCTGTTCGCGCACCAGACAACGCGTTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGGTTCAATCCCGGATTTAGCCCCAGCCCGATGtt                                                                                                                                                                                                        <  2:664746/233‑1 (MQ=255)
                                                    gTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGCGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTACCTTATCAATCCCCCTCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCCGTTTCAATCCCGGATTTAGCCCCAGCCCGATGtt                                                                                                                                                                                                        >  2:291200/1‑251 (MQ=255)
                                                    gTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGACCGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGtt                                                                                                                                                                                                        >  1:244066/1‑251 (MQ=255)
                                                    gTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGtt                                                                                                                                                                                                        >  1:134053/1‑251 (MQ=255)
                                                         cACCAGCGTCAGCATCTCCGTACGTACCGCCGGATCGGGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCGTACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGGGCATTTTCTCCGGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTcggc                                                                                                                                                                                                    <  2:374851/250‑4 (MQ=255)
                                                                  cagcacctcctaccGTAACGCCGGACCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCGGTTCGGGCACAAGACAACGCGCTAACGTCACGCGCTGTCGCGGACCGCCGGAAAGCTCGCCAGGTACCCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTc                                                                                                                                                                                           <  2:256999/237‑1 (MQ=255)
                                                                          ccTGACGTAACGCCGGTTCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTg                                                                                                                                                                                  >  1:524112/1‑251 (MQ=255)
                                                                          ccTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTg                                                                                                                                                                                  <  1:68298/251‑1 (MQ=255)
                                                                              aCGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACa                                                                                                                                                                              <  1:166643/251‑1 (MQ=255)
                                                                                 gaaCGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGCTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGGCGCTGACCGCCGGAAAGCTCGCCAGGTAACGGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAAACCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGgt                                                                                                                                                                           <  2:33074/250‑1 (MQ=255)
                                                                                   actccgcaTCGAGCGCAAAGAACGGGTCACCGGCCAATAAAATCGGCCGTCCGCGCACAAAAAACCGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCAGTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGgttg                                                                                                                                                                         <  2:246578/244‑1 (MQ=255)
                                                                                            cGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTg                                                                                                                                                                 <  1:453347/250‑1 (MQ=255)
                                                                                                 gcAGAGAACGGTTCAGCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAg                                                                                                                                                           <  1:677608/251‑1 (MQ=255)
                                                                                                      gaACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGGTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCg                                                                                                                                                      <  1:697199/251‑1 (MQ=255)
                                                                                                         cGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGac                                                                                                                                                    >  1:83218/1‑250 (MQ=255)
                                                                                                            ttCATCGAGCAATAAAATCGGCTGGTCCCGCACCAGACAACGCGCGAACGCCGCGCGTTGTCGCTGCCCGCCGGAAAGCTCGCCAGGAAACCGCGCCAGTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAAGCCCGGATTTAGCCCCAGCCCGATGTTCTGCTCGACCGTCAGGTGACTGAACAGGTTGTTTCCCTGAACCAGCATCGACAcc                                                                                                                                                 <  2:215659/250‑1 (MQ=255)
                                                                                                                 cGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACgg                                                                                                                                           <  1:15894/251‑1 (MQ=255)
                                                                                                                    gCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTATCGCTGACCGCCGGCCAGCTCGCCAGGTAACCACGCCATTAGATTATCAATCCCCATCTGGCGGGCGATATCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCCATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACGGGTTGTTTTCCTTCAACAGCCCCGCCCCCGGACGGCgg                                                                                                                                        >  2:103855/1‑250 (MQ=255)
                                                                                                                      aaTAAAATCGGCTGCACGCACACCAGACAACGCGCTAACGCCAAGCGCTGTCGCTGCACGCCGGAAAGCTCGCCAGGTATCCGCGCCATTAAATTATCAATCCCCATCTTGCGGGAGATAGAGTGCATTTTCTCCTGCTGGGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCAACCCAGAAGTTCTGCGCGACCGTCAGGTGGCTGCACAGGTATTTTACCTGTAACAGCATCggaaacggaaggcgtgc                                                                                                                                      >  2:606357/1‑235 (MQ=255)
                                                                                                                                  tGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGAGGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCTGGACggcggcgtcgg                                                                                                                          >  1:20449/1‑246 (MQ=255)
                                                                                                                                   gTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGAcgcgatcggaat                                                                                                                          >  2:150737/1‑240 (MQ=255)
                                                                                                                                   gTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGAcg                                                                                                                                    <  1:150737/240‑1 (MQ=255)
                                                                                                                                    ttCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACggcggcgtcgtgg                                                                                                                        >  1:116682/1‑244 (MQ=255)
                                                                                                                                        cgcaCCAGACACCGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACggcggcgtcgtggtgtg                                                                                                                    <  2:116682/248‑12 (MQ=255)
                                                                                                                                          caCCAGACAACGTGCTAGCGCCACTCGCTTTCGCTTACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAACTATCAATCCCCATCCTGCGGGCGATAGCGTGCATTTTCTCCTGCTGTTCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGAcggcgtcgtcgtggtgtgat                                                                                                                  <  2:446767/250‑16 (MQ=255)
                                                                                                                                                aCAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGGGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGCCGGCGTGCCGGCGGCGTCGTGGTGTGATCAACGc                                                                                                            >  2:161604/1‑251 (MQ=255)
                                                                                                                                                    cgtggtAACGCCACGAGCTGTCGTTACCCGCCGGAAAGCTCGCCAGGTAACCGGGCCATTTAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCGCCTTCTGTGCCGCGGTCAGTTTCAATCCCGGATTTAGCTCCAGCCCGATGTTCTGCGCGACCGTCAGTTGACTGAACAGGTTGTTTTCCTGAAACAGCATAGACCCCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGccggc                                                                                                        <  2:538745/246‑4 (MQ=255)
                                                                                                                                                    cgcgCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGccgtc                                                                                                        >  1:648462/1‑248 (MQ=255)
                                                                                                                                                    cgcgCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGccgtc                                                                                                        >  1:93357/1‑248 (MQ=255)
                                                                                                                                                     gcgcTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCAGTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAAAAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCg                                                                                                       <  1:78702/251‑1 (MQ=255)
                                                                                                                                                     cccgtAACGCCACGCGCTGGCGCTGACCGCCGGAAAACTCCGCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCg                                                                                                       <  1:146009/247‑1 (MQ=255)
                                                                                                                                                      cgcTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCg                                                                                                       >  1:454654/1‑250 (MQ=255)
                                                                                                                                                       gcTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGAt                                                                                                     <  1:630387/251‑1 (MQ=255)
                                                                                                                                                         tAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCg                                                                                                                                                >  2:8071/1‑206 (MQ=255)
                                                                                                                                                         tAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCg                                                                                                                                                <  1:8071/206‑1 (MQ=255)
                                                                                                                                                          aaCGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAg                                                                                                   <  1:389903/250‑1 (MQ=255)
                                                                                                                                                            cGCCACGCGCTGGCGCTGACCGCCGGAAAGCGCGCCAGGTAACCGCGCCATTAAATTATCAAGCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCa                                                                                                <  2:134053/251‑1 (MQ=255)
                                                                                                                                                             gCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGCCACCGGACGGCGTGACGGCGGCGTCGAGGTGTGATCAACGccgccggtagtccg                                                                                               >  2:486554/1‑239 (MQ=255)
                                                                                                                                                             gCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAg                                                                                               <  1:336952/251‑1 (MQ=255)
                                                                                                                                                             gCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAg                                                                                               <  1:139350/251‑1 (MQ=255)
                                                                                                                                                             ccgccgggctgtggctgACCGCCGGAAAGCTCGGCAGGGTACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGGGCATTTTCTCCTGCTGTGCCGCGTTCAGGTTCAATCCCGGATGTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACGGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAg                                                                                               <  2:628801/244‑1 (MQ=255)
                                                                                                                                                              ccACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGc                                                                                              <  1:399025/251‑1 (MQ=255)
                                                                                                                                                               catgctctgtcgctgCCCGCCGGAAAGCTCGCCATGTAACTGCGCCATTAAATGACCATTCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACCCCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCg                                                                                             <  2:203184/245‑1 (MQ=255)
                                                                                                                                                               cacgcgctgtcgctgACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCg                                                                                             <  1:522392/247‑1 (MQ=255)
                                                                                                                                                               cacgcgctgtcgctgACCGCCGGAAAGCGCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCAGCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCg                                                                                             <  1:316903/247‑1 (MQ=255)
                                                                                                                                                                   cgctgtcgctgACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAAcc                                                                                         <  1:177329/251‑1 (MQ=255)
                                                                                                                                                                          gctgACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCgctggct                                                                                  <  1:122261/251‑1 (MQ=255)
                                                                                                                                                                              aCCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCg                                                                              >  1:337658/1‑251 (MQ=255)
                                                                                                                                                                                cgccgGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTc                                                                            <  1:209334/251‑1 (MQ=255)
                                                                                                                                                                                 cccgGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCa                                                                           <  1:367126/250‑1 (MQ=255)
                                                                                                                                                                                         gCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCg                                                                   <  1:317465/251‑1 (MQ=255)
                                                                                                                                                                                              ccaggttaCCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGTGAACCGCTGGCTGGCTTCAGAAAACCTGCGAt                                                              <  1:606357/244‑1 (MQ=255)
                                                                                                                                                                                                gggTAACCGCGCCATTAAATTATCAACCCCCATCTGGCTGGCGATCGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGAGTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGGCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCa                                                            <  2:545725/250‑1 (MQ=255)
                                                                                                                                                                                                    acccGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGAtt                                                                                                                                                                                                                          <  2:349109/87‑1 (MQ=255)
                                                                                                                                                                                                    aaCCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGAtt                                                                                                                                                                                                                          >  1:349109/1‑89 (MQ=255)
                                                                                                                                                                                                    aaCCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGGGCATTTTCTCCTGCTGTGCCGCGTTCAGCTTCAATCCCGGATTTAGCCCAAGCCCGAAGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGGTCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCaagtt                                                        >  2:233736/1‑248 (MQ=255)
                                                                                                                                                                                                     aCCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACCGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTc                                                       >  2:274904/1‑251 (MQ=255)
                                                                                                                                                                                                     aCCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTc                                                       >  1:618288/1‑251 (MQ=255)
                                                                                                                                                                                                       cgcgCCATTAAATTATCAATCCCCATCTGGCGGGCGATTGCGTGCATTGTCGCATGCGGTTCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCCGCCCGATGTACTGCGCGACCGTCAGGCGACTGGACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCGGTCGTCAGCAACCCGGCGTTCAAGTTcag                                                     <  2:599803/251‑1 (MQ=255)
                                                                                                                                                                                                       cgcgCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTcag                                                     <  1:498405/251‑1 (MQ=255)
                                                                                                                                                                                                           ccATTAAATTCTCAATCCCCATCGGGCGGGCGATAGCGTGCATTTTCTCCTGCTTTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCACAGCCCGATGTTCTTCGCGACCTGCAGTTGACTGAACAGGTTGTTTTCCTGAAACAGCACCGACACCGGACGGCGTGCCGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTcagcatg                                                 <  2:215203/251‑3 (MQ=255)
                                                                                                                                                                                                            cATTAAATTATCAATCCCCATCTGGCGTGCGATAGCGTTCACTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGTTAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAggg                                                <  2:473237/251‑1 (MQ=255)
                                                                                                                                                                                                               tAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTa                                              <  1:505551/250‑1 (MQ=255)
                                                                                                                                                                                                                 aaTTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACtt                                           <  1:22573/251‑1 (MQ=255)
                                                                                                                                                                                                                  aTTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGCTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACttt                                          >  1:83352/1‑251 (MQ=255)
                                                                                                                                                                                                                   ttATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACtttt                                         >  1:653675/1‑251 (MQ=255)
                                                                                                                                                                                                                   ttATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAAAGCCGTCGATAGGCAGCGAACCGCTGGCTGGGTTCGGAAAAACAGCGAACGAGCTCAGCAGGGTACtttt                                         >  2:600580/1‑251 (MQ=255)
                                                                                                                                                                                                                    tATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGGTGTTTTCCTGAAACAGCATCGacac                                                                                                                                                  <  1:509143/145‑1 (MQ=255)
                                                                                                                                                                                                                       cAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTAcc                                      <  1:229188/250‑1 (MQ=255)
                                                                                                                                                                                                                           ccccTTCTGGCGGGCCATAGCGTGCATTTTCTCCTGCTGGGCCGCGTTCAGTTTCAATCCCGGATTGAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACGGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGAGAGCAAGTTCAGCAGGGTACTTTTACCGgcgc                                 <  2:244066/251‑1 (MQ=255)
                                                                                                                                                                                                                            cccATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGGCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGcgcc                                >  2:328623/1‑251 (MQ=255)
                                                                                                                                                                                                                               agcTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCt                             <  1:637950/249‑1 (MQ=255)
                                                                                                                                                                                                                                  tGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTgg                          >  1:647474/1‑251 (MQ=255)
                                                                                                                                                                                                                                  tGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAAAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTCCCTGCGCCGCTTgg                          >  2:664217/1‑251 (MQ=255)
                                                                                                                                                                                                                                   ggcgggcgATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTg                           <  1:87633/249‑1 (MQ=255)
                                                                                                                                                                                                                                          gATAGCGTGCAGTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGa                  <  1:103855/251‑1 (MQ=255)
                                                                                                                                                                                                                                                gTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAg                                                                                                   <  2:275797/164‑1 (MQ=255)
                                                                                                                                                                                                                                                gTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAg                                                                                                   >  1:275797/1‑164 (MQ=255)
                                                                                                                                                                                                                                                gTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGcc             >  1:164242/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    aTTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTg        <  1:401387/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    aTTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTg        <  1:161604/251‑1 (MQ=255)
                                                                                                                                                                                                                                                          cccTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTcgcc  <  2:647474/251‑1 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                          
TTCCACGCTGTGCGACACCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTACCAGACAACGCGCTAACGCCACTCGCTGTCGCTGACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCC  >  NC_000913/72355‑72855

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: