Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 864,741 Δ1 bp 5.3% coding (389/750 nt) moeB ← molybdopterin synthase sulfurylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913864,7410T. 5.3% 257.9 / 6.3 76coding (389/750 nt)moeBmolybdopterin synthase sulfurylase
Reads supporting (aligned to +/- strand):  ref base T (29/41);  new base . (2/2);  total (33/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

GACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAATTACGCCGATCAACGGTGCCATTACGCCTGCTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAG  >  NC_000913/864458‑865027
                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                              
gACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAATTACGCCGATCAACGGTGCCATTACGCCTGCTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTa                                                                                                                                                                                                                                                                                           <  1:74613/289‑1 (MQ=255)
gACGATGTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAATTACGCCGATCAACGGTGCCATTACGCCTGCTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTa                                                                                                                                                                                                                                                                                           <  1:150238/289‑1 (MQ=255)
       ttCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAATTACGCCGATCCACGGTCCCATTACGCCTGCTTCCACGCAGGTTAATGCCTTTTCACCAAACAAACGGCGAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTCACGACGGTGCTTTGACCTTCCATCCGAATTGCCGCCCCGGAACCCTGCGGTACTTTCGCGGCCAAACAGCCTGCGTTCGGTTGATTACGTACCGCAACGTTATACGTCAAGTCGa                                                                                                                                                                                                                                                                                    >  1:12376/1‑289 (MQ=255)
                     tttGCCATAACCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAATTACGCCGATCAACGGTGCCATTACGCCTGCTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCgtg                                                                                                                                                                                                                                                                       <  1:110715/288‑1 (MQ=255)
                          cATAACCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAATTACGCCGATCAACGGTGCCATTACGCCTGCTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAg                                                                                                                                                                                                                                                                  >  1:198851/1‑288 (MQ=255)
                          cATAACCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAATTACGCCGATCAACGGTGCCATTACGCCTGCTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTACCGTACAGTCGAGCACCAGATCGTGTTCAGca                                                                                                                                                                                                                                                                >  2:505071/1‑289 (MQ=255)
                          cATAACCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAATTACGCCGATCAACGGTGCCATTACGCCTGCTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTACCGTACAGTCGAGCACCAGATCGTGTTCAGca                                                                                                                                                                                                                                                                >  2:298148/1‑289 (MQ=255)
                                 tGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAATTACGCCGATCAACGGTGCCATTACGCCTGCTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGcaatcaat                                                                                                                                                                                                                                                          <  1:167251/289‑1 (MQ=255)
                                   ccAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAATTACGCCGATCAACGGTGCCATTACGCCTGCTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAAtg                                                                                                                                                                                                                                                         >  1:427727/1‑288 (MQ=255)
                                                cGCTTCCATTGCCTGCAACGAACCAATTACGCCGATCAACGGTGCCATTACGCCTGCTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCg                                                                                                                                                                                                                                           <  1:91664/289‑1 (MQ=255)
                                                       aTTGCCTGCAACGAACCAATTACGCCGATCAACGGTGACATTACGCCTGCTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCcag                                                                                                                                                                                                                                   <  2:425767/289‑1 (MQ=255)
                                                         tGCCTGCAACGAACCAATTACGCCGATCAACGGTGCCATTACGCCTGCTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCcagc                                                                                                                                                                                                                                  <  2:11446/289‑1 (MQ=255)
                                                             tGCAACGAACCAATTACGCCGATCAACGGTGCCATTACGCCTGCTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCa                                                                                                                                                                                                                                                                                 <  2:323709/238‑1 (MQ=255)
                                                             tGCAACGAACCAATTACGCCGATCAACGGTGCCATTACGCCTGCTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCa                                                                                                                                                                                                                                                                                 >  1:323709/1‑238 (MQ=255)
                                                             tGCAACGAACCAATTACGCCGATCAACGGTGCCATTACGCCTGCTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGt                                                                                                                                                                                                                               >  1:160492/1‑288 (MQ=255)
                                                               cAACGAACCAATTACGCCGATCCACGGTGCCATTACGCCTGCTTGCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGCTTACGTACCGCAACGTTAGCCGTACCGTCGAGCACCAGATCGTGTTACGCAATCAATGCTGCAAGTTCTGCGTCCTCCAGCCGtcc                                                                                                                                                                                                                             >  1:124346/1‑286 (MQ=255)
                                                               cAACGAACCAATTACGCCGATCAACGGTGCCATTACGCCTGCTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCa                                                                                                                                                                                                                            >  1:475891/1‑289 (MQ=255)
                                                                aaCGAACCAATTACGCCGATCAACGGTGCCATTACGCCTGCTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCa                                                                                                                                                                                                                            <  2:427727/288‑1 (MQ=255)
                                                                     aCCAATTACGCCGATCAACGGTGCCATTACGCCTGCTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTg                                                                                                                                                                                                                         >  1:286764/1‑286 (MQ=255)
                                                                               ccGATCAACGGTGCCATTACGCCTGCTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATcg                                                                                                                                                                                                            <  2:354548/289‑1 (MQ=255)
                                                                                  aTCAACGGTGCCATTACGCCTGCTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGa                                                                                                                                                                                                         >  2:43614/1‑289 (MQ=255)
                                                                                     aaCGGTGCCATTACGCCTGCTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGatat                                                                                                                                                                                                      <  1:89057/289‑1 (MQ=255)
                                                                                                    ccTGCTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCCGTTGATTACGTACCGCAACGTTACCCGTACATTCGAGCACCAGATCGCGTTCAGCAATCAATGCTGCCAGTTCTGCGTCCTCCAGCAGTGCATTGACCGGCGTAATCGCGTTGTGTGGGTTGATCcgcg                                                                                                                                                                                       >  1:48687/1‑289 (MQ=255)
                                                                                                    ccTGCTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTa                                                                                                                                                                                                                                                                                           <  2:279265/189‑1 (MQ=255)
                                                                                                    ccTGCTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTa                                                                                                                                                                                                                                                                                           >  1:279265/1‑189 (MQ=255)
                                                                                                    ccTGCTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTc                                                                                                                                                                                                                                                                                      <  2:464086/194‑1 (MQ=255)
                                                                                                    ccTGCTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTc                                                                                                                                                                                                                                                                                      >  1:464086/1‑194 (MQ=255)
                                                                                                    ccTGCTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTACCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGt                                                                                                                                                                                      <  1:505071/289‑1 (MQ=255)
                                                                                                    ccTGCTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTACCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGt                                                                                                                                                                                      <  1:298148/289‑1 (MQ=255)
                                                                                                      tGCTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTc                                                                                                                                                                                     >  2:33742/1‑289 (MQ=255)
                                                                                                        cTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGCTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATcg                                                                                                                                                                                                            <  2:202871/264‑1 (MQ=255)
                                                                                                        cTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGCTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATcg                                                                                                                                                                                                            >  1:202871/1‑264 (MQ=255)
                                                                                                                 aGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCAc                                                                                                                                                                           <  1:43614/288‑1 (MQ=255)
                                                                                                                         gCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGAtt                                                                                                                                                                  <  2:198851/289‑1 (MQ=255)
                                                                                                                             tttCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAATGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCa                                                                                                                                                               >  1:140845/1‑288 (MQ=255)
                                                                                                                             tttCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTg                                                                                                                                                                                                                                             <  2:45817/210‑1 (MQ=255)
                                                                                                                             tttCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTg                                                                                                                                                                                                                                             >  1:45817/1‑210 (MQ=255)
                                                                                                                             tttCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCa                                                                                                                                                               >  2:48367/1‑288 (MQ=255)
                                                                                                                                 ccaaaacaaacGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTc                                                                                                                                                          <  1:235190/286‑1 (MQ=255)
                                                                                                                                 aCCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTc                                                                                                                                                          <  2:160492/289‑1 (MQ=255)
                                                                                                                                        aaacGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGcc                                                                                                                                                   <  1:506460/289‑1 (MQ=255)
                                                                                                                                            ggCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGAc                                                                                                                                               <  1:35573/289‑1 (MQ=255)
                                                                                                                                                 aGGCAGCGATAGGACGGTTCACCGTCCTTATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTgg                                                                                                                                          <  1:97647/289‑1 (MQ=255)
                                                                                                                                                 aGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTgg                                                                                                                                          <  1:48367/289‑1 (MQ=255)
                                                                                                                                                    cgtcGATAGGACGGTTCACCGACCTGATAAGTAAAGACGGTGATTTGACCTTCCACACGAATGGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCa                                                                                                                                        <  2:12376/285‑1 (MQ=255)
                                                                                                                                                    cAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCa                                                                                                                                        <  1:33742/288‑1 (MQ=255)
                                                                                                                                                       cGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTc                                                                                                                                                                                                                                                                                      >  2:264778/1‑143 (MQ=255)
                                                                                                                                                       cGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTc                                                                                                                                                                                                                                                                                      <  1:264778/143‑1 (MQ=255)
                                                                                                                                                        gATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAATGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCa                                                                                                                                     <  2:140845/287‑1 (MQ=255)
                                                                                                                                                              aCGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCa                                                                                                                              <  2:475891/288‑1 (MQ=255)
                                                                                                                                                              aCGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCa                                                                                                                              <  2:286764/288‑1 (MQ=255)
                                                                                                                                                                   tCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTc                                                                                                                         <  1:126199/288‑1 (MQ=255)
                                                                                                                                                                     aCCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCa                                                                                                                                     <  2:195086/274‑1 (MQ=255)
                                                                                                                                                                     aCCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCa                                                                                                                                     >  1:195086/1‑274 (MQ=255)
                                                                                                                                                                                        aCGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCg                                                                                                                                                                                                                  >  1:447391/1‑178 (MQ=255)
                                                                                                                                                                                        aCGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCg                                                                                                                                                                                                                  <  2:447391/178‑1 (MQ=255)
                                                                                                                                                                                             gATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAgg                                                                                                                                                                                  <  2:473827/205‑1 (MQ=255)
                                                                                                                                                                                             gATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAgg                                                                                                                                                                                  >  1:473827/1‑205 (MQ=255)
                                                                                                                                                                                                        ccATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCCTCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGaact                                                                                    >  1:195109/1‑286 (MQ=255)
                                                                                                                                                                                                        ccATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGCGCGAAACCGTGTCGAAGTc                                                                                   >  2:459528/1‑289 (MQ=255)
                                                                                                                                                                                                        ccATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTc                                                                                   >  2:236661/1‑289 (MQ=255)
                                                                                                                                                                                                                 tggccggtacGGAAACCAGGCGTACCGTCCCGGCAAAAAGGCCTGAGTTCAGTGGATTACGTACCTCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGt                                                                          <  1:459528/280‑1 (MQ=255)
                                                                                                                                                                                                                                            tCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGtt                                                                                                                  >  1:244944/1‑222 (MQ=255)
                                                                                                                                                                                                                                            tCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGtt                                                                                                                  <  2:244944/222‑1 (MQ=255)
                                                                                                                                                                                                                                             cgcgGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCcagc                                                                                                                                                                                                                                  >  1:274302/1‑109 (MQ=255)
                                                                                                                                                                                                                                             cgcgGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCcagc                                                                                                                                                                                                                                  <  2:274302/109‑1 (MQ=255)
                                                                                                                                                                                                                                                          ccTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGCTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATc                                                                                                                                      >  2:51632/1‑188 (MQ=255)
                                                                                                                                                                                                                                                          ccTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGCTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATc                                                                                                                                      <  1:51632/188‑1 (MQ=255)
                                                                                                                                                                                                                                                                 tCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGcc                          >  2:431021/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                    gTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTgg                                                                                                                                          <  2:40750/174‑1 (MQ=255)
                                                                                                                                                                                                                                                                    gTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTgg                                                                                                                                          >  1:40750/1‑174 (MQ=255)
                                                                                                                                                                                                                                                                     ttGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCGTGCCAGATACTGCGAGGCTGCACAGCCGAgg                      >  2:84342/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                      tGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAgg                      <  1:431021/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                ccGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGcc           <  1:248068/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      cGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTa      >  1:360913/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                         tATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAg  <  2:247105/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                              
GACGATTTTCCCGCTGGCCGGTTTGCCATAACCTGCCAGCATTTTGATCGCTTCCATTGCCTGCAACGAACCAATTACGCCGATCAACGGTGCCATTACGCCTGCTTCCACGCAGGTTAATGCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTTCGCGGCAAAACAGCCTGCGTTCAGTTGATTACGTACCGCAACGTTATCCGTACAGTCGAGCACCAGATCGTGTTCAGCAATCAATGCTGCAAGTTCTGCGTCATCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAG  >  NC_000913/864458‑865027

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