Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,832,797 T→G 6.5% N581H (AAC→CAC)  yicI ← putative alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,832,7970TG 6.5% 183.4 / 3.9 64N581H (AAC→CAC) yicIputative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base T (25/33);  new base G (2/2);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.22e-01

TAATCGGGACGTTGATCGTTGTTGCCCAGCGCCAGTAGAGTGT‑TATCACGCACATAAACGGGCAGACTCAGGAAGCCGTGCTGCTGTTTATGCCAGCGACTACCGTCGAGTTCATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCC  >  NC_000913/3832499‑3833097
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                            
tAATCGGGACGTTGATCGTTGTTGCCCAGCGCCAGTAGAGTGT‑TATCACGCACATAAACGGGCAGACTCAGGAAGCCGTGCTGCTGTTTATGCCAGCGACTACCGTCGAGTTCATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCg                                                                                                                                                                                                                                                                                                            <  2:111529/301‑1 (MQ=255)
          tttGATCGTTGTTGCCCAGCGCCAGTAGATTGT‑TATCACGCACATAAACGGGCAGACTCAGGAAGCCGTGCTGCTGTTTATGCCAGCGACTACCGTCGAGTTCATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTGCGAACgcgcggg                                                                                                                                                                                                                                                                                                  <  2:169130/300‑4 (MQ=255)
            gggttgttgttGTCCAGCGGCAGTA‑TGTGTGTATCACGCAGATATATGGGCAGACATAGGAAGGCGGGATGCTGTGTATGGCAGCGACTTCCGGCGAGTTCGTCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACCGCACATCGCCCGCTGCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTGCGAAGgcgcggggt                                                                                                                                                                                                                                                                                                <  1:317419/296‑6 (MQ=255)
                        cccAGCGCCAGTAGAGTGT‑TATCACGCACATAAACGGGCAGACTCAGGAAGCCGTGCTGCTGTTTATGCCAGCGACTACCGTCGAGTTCATCGTTGTTCCCCAGGTGTGTCCAGCGACCTTCCGGCATGTAGAACTGCACATCGCCCGCTTCAGTGCACACCGCCGCAACCACCACGCTTTCGCATAACATGTATTGACGGTCAAGGTAATCACAAGCCGGCTCATCCGTGAACTCCATCATCAT‑GTCCCGCATCATTCGCGTACACCCCCCGTCCGCAccccaaacaccccaacacagaa                                                                                                                                                                                                                                                                                    >  1:283645/1‑280 (MQ=255)
                                  gTAGAGTGT‑TATCACGCACATAAACGGGCAGACTCAGGAAGCCGTGCTGCTGTTTATGCCAGCGACTACCGTCGAGTTCATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCCTCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGCTCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCCCAGCGTCACGATACAAAAAAGGcatcat                                                                                                                                                                                                                                                                          >  1:150847/1‑301 (MQ=255)
                                         gt‑TATCACGCACATAAACGGGCAGACTCAGGAAGCCGTGCTGCTGTTTATGCCAGCGACTACCGTCGAGTTCATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCACCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCaa                                                                                                                                                                                                                                                                   >  2:224483/1‑300 (MQ=255)
                                               cccGCACATAAACGGGCAGACTCAGGAAGCCGTGCTGCTGTTTATGCCAGCGACTACCGTCGAGTTCATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAg                                                                                                                                                                                                                                                              <  1:83194/299‑1 (MQ=255)
                                                      aTAAACGGGCAGACTCAGGAAGCCTTGCTGCTGTTTATGCCAGCGACTACCGTCGAGTTCATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTg                                                                                                                                                                                                                                                       <  2:280201/301‑1 (MQ=255)
                                                      aTAAACGGGCAGACTCAGGAAGCCGTGCTGCTGTTTATGCCAGCGACGACCGTCGAGTTCATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTg                                                                                                                                                                                                                                                       <  2:437972/301‑1 (MQ=255)
                                                       tAAACGGGCAGACTCAGGAAGCCGTGCTGCTGTTTATGCCAGCGACTACCGTCGAGTTCATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCCTCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCGTGCGGCATTTCAGTTGCATga                                                                                                                                                                                                                                                      >  2:350027/1‑301 (MQ=255)
                                                             cttctttccctccactacgcttttccgttttGCCAGCGACTACCGGCGAGTTCATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAgc                                                                                                                                                                                                                                                <  1:59810/271‑1 (MQ=255)
                                                             ctctttccctacacgacgctcttccgagcttGCCAGCGACTACCGTCGAGTTCATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAgc                                                                                                                                                                                                                                                <  1:103820/271‑1 (MQ=255)
                                                                 gACTCAGGAAGCCGTGCTGCTGTTTATGCCAGCGACTACCGTCGAGTTCATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCCTACCCCGCACGTTCGCACGCGCAGCCTCACGCAACAACAACAGCACCATCAGGCAATTCAATTTCGTTCAGAAACCcac                                                                                                                                                                                                                                            >  1:606136/1‑301 (MQ=255)
                                                                    tCAGGAAGCCGTGCTGCTGTTTATGCCAGCGACTACCGTCGAGTTCATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGcaccac                                                                                                                                                                                                                                         <  2:376879/301‑1 (MQ=255)
                                                                       ggAAGCCGTGCTGCTGTTTATGCCAGCGACTACCGTCGAGTTCATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACACTCACGTTGTCGCCTAACATGTATTGCCGGTCAAGGTAATCACAAGCCGGATCGTCCGTGAACTCCATCAACAT‑CGCCCGACTCATCCGCGTACCCCGCGCGTACGCACGCGCAGCTTCACGATCCATATACGGCCACATGCGGCATTTCACTTACGAgaaaacgcaacccacata                                                                                                                                                                                                                                      >  1:38812/1‑286 (MQ=255)
                                                                         aaGTCGTGCTGTTTTTTTTGCCAGCGACTACCGTTGGGTTCATCGTTGTGCTACAGGTGTGTCCAGCGACCTGCCGGCAGGTAGAACTGCACATCACCCGCTTCAGTGAACACCGGCGCAACCATCACGTGGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCAGCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATcac                                                                                                                                                                                                                                    <  1:224240/301‑1 (MQ=255)
                                                                            cgtggctggtGTTTATGCCAGCGGCTGCCGTCGGGTTCTGCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTGGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAgg                                                                                                                                                                                                                                 <  1:566998/297‑1 (MQ=255)
                                                                              gTGCTGCTGTTTATGCCAGCGACTACCGTCGAGTTCATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAggcc                                                                                                                                                                                                                               <  2:632304/301‑3 (MQ=255)
                                                                                    ttgTTTATGCCAGCGACTACCGTCGAGTTCATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACtcatca                                                                                                                                                                                                                         <  2:307763/300‑1 (MQ=255)
                                                                                       tttATGCCAGCGACTACCGTCGAGTTCATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCg                                                                                                                                                                                                                      <  2:150998/301‑1 (MQ=255)
                                                                                         ccgtgcgcttttccgatctgAGTTCATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTa                                                                                                                                                                                                                    <  2:371776/282‑1 (MQ=255)
                                                                                           tGCCAGCGACTACCGTCGAGTTCATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAgcagatcggaagagcacacgtctgaactccaa                                                                                                                                                                                                                  >  2:103944/1‑271 (MQ=255)
                                                                                            gCCAGCGACTACCGTCGAGTTCATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCACTTCAGTTGCGTGAAGAAGCCCAACACAACACAGGaaccaccaaccaagcc                                                                                                                                                                                                                 >  1:106829/1‑286 (MQ=255)
                                                                                              actcgacggtagTCGAGTTCATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGccc                                                                                                                                                                                                               >  2:621139/12‑301 (MQ=255)
                                                                                                       ccGTCGAGTTCATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGAACGTCCGGGCACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGNCGCTTCACGATACAGATACGGCATCATCACGCCTCTCAGTAGCGTGACAAAGCGCAACACACCACAGAACCCCACACCGTACACCCAcaccccac                                                                                                                                                                                                      >  1:385832/1‑294 (MQ=255)
                                                                                                            gAGTTCATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACCGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCAACGTaagaacggaagagcgccgtg                                                                                                                                                                                                 >  1:371280/1‑282 (MQ=255)
                                                                                                              gTTCATCGTGGTGCCACAGGTGTGTCCAGCTACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGGGAAAACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGGCAAGGTAATCACAAGCCGGATCGTCCGGGAACGCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGTAGAAGCGCACCACATTACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGa                                                                                                                                                                                               <  1:349564/301‑1 (MQ=255)
                                                                                                                 cATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGCATCAACATCCTAGGCCCACGGcaaccgaaaaacatt                                                                                                                                                                                            >  1:263442/1‑288 (MQ=255)
                                                                                                                    cGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATCCGGCATCATGCGGCATTTCAGTTGCGTGACAAAGCGCACCACATCACAAGACTCATCAACCTCAGCCCACGGCACACGATaacacttact                                                                                                                                                                                         >  1:37294/1‑293 (MQ=255)
                                                                                                                           ccACAGGTGTTTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCgtgt                                                                                                                                                                                  <  2:263710/301‑1 (MQ=255)
                                                                                                                            cacaGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCCCGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCCCATCACAGGACCCATCATCGTAGGCCCACGGCACACGATAACAttaccaaccatgtt                                                                                                                                                                                 >  1:52149/1‑289 (MQ=255)
                                                                                                                                     gtCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCCCCACATCACAGGACTCACCATCGTAAGCCCACGGCACACCATAAGATTTACTACCGTGTAAACAGCtct                                                                                                                                                                        >  1:447788/1‑299 (MQ=255)
                                                                                                                                             aCCTTCTGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGGATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGgag                                                                                                                                                                <  2:435845/301‑1 (MQ=255)
                                                                                                                                                   cGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCCCAAGACTCATCCTCGTAAGCCCACGGCACCAGATAAGATTTGCTAACGTGTAAACGGCAAAGGATggcaaacaaa                                                                                                                                                          >  1:170639/1‑293 (MQ=255)
                                                                                                                                                    ggCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATATAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCCCACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGcaaaa                                                                                                                                                         >  2:550172/1‑300 (MQ=255)
                                                                                                                                                                     aTCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACAATAAGATTTGCTACCGTGTACAAGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAACGt                                                                                                                                        >  2:161658/1‑301 (MQ=255)
                                                                                                                                                                        gCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCCCGGCACACGATAAGATTTGCTACAGTGTAAACGGCTATGGATGGAAAACAAACCAAACGCGCACAAGctctat                                                                                                                                     >  1:95470/1‑296 (MQ=255)
                                                                                                                                                                          ccgtgtCAGGGAATACGGGCGCAACCATCAGGTGGTCGCCTTACATGTATTTACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACGCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGAGACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTa                                                                                                                                   <  1:592760/296‑1 (MQ=255)
                                                                                                                                                                              ttCATTTAACACCGGCGCAACCATCATGGTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCTTACCCCGCTCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCTTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAAcgctggtgcgcgtttggtttgc                                                                                                                               <  2:38864/301‑20 (MQ=255)
                                                                                                                                                                               tCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACCGGACTCCTCATCCTAGGCCCACGGCACCCCATAAGATTTGCTACCGTTTAAACAGCTATGCCTGGACAGCCAAACAAACGCGACCCACAGATTGTAAACGt                                                                                                                              >  1:448318/1‑301 (MQ=255)
                                                                                                                                                                                         ccGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAAAGCGAACCAGCGTTTGTAAAAGACCGCCGGacg                                                                                                                    >  1:278354/1‑299 (MQ=255)
                                                                                                                                                                                          cGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGAGCGCCGGAGCg                                                                                                                   >  2:201955/1‑301 (MQ=255)
                                                                                                                                                                                                  cctgcatgtggTCGCCGAACAGGTATGGTCGGTAAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCAGCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGAACGCGCAGCTTCACGATATAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTc                                                                                                           <  1:549499/291‑1 (MQ=255)
                                                                                                                                                                                                   cATTACGTTGTCGCTTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCa                                                                                                          <  2:106953/301‑1 (MQ=255)
                                                                                                                                                                                                   cATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCAATGGCAGAAGAGCAAACCAAACGCGCACCAGCATTTGTAAACGAGAGCAGAAGCGGTATTTACa                                                                                                          >  2:242466/1‑301 (MQ=255)
                                                                                                                                                                                                     acACGTTGTCGCTTAAAATGTATTGACGGTAAAGGTAATCACAAGCCGGTTCGGCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCaaa                                                                                                        <  1:620281/300‑1 (MQ=255)
                                                                                                                                                                                                      taCGTTGTCGCCTAACATTTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAg                                                                                                       <  2:170836/300‑1 (MQ=255)
                                                                                                                                                                                                            gTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCCTGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCCTTTGTAAACGTCCGCCGGAGACGCATTTTCAAAccccccc                                                                                                 >  1:399585/1‑300 (MQ=255)
                                                                                                                                                                                                              cGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCCAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGAGGTAATTTCAAACCCACCCat                                                                                               >  1:266236/1‑301 (MQ=255)
                                                                                                                                                                                                                   atcttGTTTTGACGGTCAAGTGAATTCCAAGCCGGATCGTCCGGGAACTCCATCA‑CATGGGCCCGCATCATCGGCGGACGCCGCGCGTTCGCACGCGCAGCGTCACGATACAGATACGGCATCATGCGTCATTGCAGTTGCGTGAAGAAGCGCACCACAGCACACGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGt                                                                                          <  1:313594/297‑1 (MQ=255)
                                                                                                                                                                                                                     tttGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTgg                                                                                        <  2:386328/299‑1 (MQ=255)
                                                                                                                                                                                                                          ttgtaCGGTCAAGTTAATCCTAATCCGGATCTTCCGTGAACGCCATCATCAT‑GGCCCGCAGCATCGGCGGACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGAGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCC‑ACGATATCggcggctttg                                                                                  <  1:416919/297‑10 (MQ=255)
                                                                                                                                                                                                                              ccGGTCAAGTTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAg                                                                               <  2:606897/300‑1 (MQ=255)
                                                                                                                                                                                                                                      ggTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAAcc                                                                       <  2:400099/301‑1 (MQ=255)
                                                                                                                                                                                                                                      ggTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACAACACAGGACTCATCATCGTAGGCCCACGGCCCACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGCGCAAACCAAACGCGCACCCGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCCCCACAAGCCAAAccc                                                                       >  1:592428/1‑301 (MQ=255)
                                                                                                                                                                                                                                       gTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCt                                                                      <  2:266504/301‑1 (MQ=255)
                                                                                                                                                                                                                                               aaGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGcc                                                              <  2:52217/301‑1 (MQ=255)
                                                                                                                                                                                                                                                      gtCGTCCGGGAACTCCTTTCTCAT‑GGCCCGCCTCCTCGGCGTACCCCTCGCGTTCGCACGCGCAGCTTCACGATACAGATACCGCATTATGCGTCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGAc                                                       <  1:242199/300‑1 (MQ=255)
                                                                                                                                                                                                                                                         gTCCGGGAACTCCATCATCAT‑GGCACGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGATACCGTGTAAACGGCTATGGCTGGAGAGAAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCAGGAGCGGTATTTTAAAAGCCGCCGCGATCGTGGCTCCAGAAGCCAAAAACTGCAAGGACAATAaaaaaa                                                    >  2:412950/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                 attCCATCATCAT‑GGTCCGAATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTAACGATCCAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGAACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACcgcg                                            <  2:31493/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                    ccATCATCCT‑GGCCTGGTTCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAg                                         <  1:201712/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 cccGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCAt                             <  2:448352/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        gacgtgtgctcttccgatcctCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCg                      <  2:412174/281‑1 (MQ=255)
                                                                                                                                                                                                                                                                                         gtatcgtgaagctgcgcgTGCGAACGCGCGGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCga                     >  2:473754/14‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                   gctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAACCCTGAAAGGCCAATAGACAAACCACCCCCCAGGCATTCCCCCAATGaaccgaacataccgtt           >  1:512374/4‑286 (MQ=255)
                                                                                                                                                                                                                                                                                                            tCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAACGGCCAATAGACAAACCACAGCGCAGGCTTTCCGCCATAGATTCgagatcgaaacaccgacgagt  >  1:411611/1‑281 (MQ=255)
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                            
TAATCGGGACGTTGATCGTTGTTGCCCAGCGCCAGTAGAGTGT‑TATCACGCACATAAACGGGCAGACTCAGGAAGCCGTGCTGCTGTTTATGCCAGCGACTACCGTCGAGTTCATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCC  >  NC_000913/3832499‑3833097

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑