Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,832,800 C→A 7.1% A580S (GCG→TCG)  yicI ← putative alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,832,8000CA 7.1% 186.5 / 5.0 66A580S (GCG→TCG) yicIputative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base C (29/32);  new base A (2/3);  total (31/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 = 3.07e-01

GTTGATCGTTGTTGCCCAGCGCCAGTAGAGTGT‑TATCACGCACATAAACGGGCAGACTCAGGAAGCCGTGCTGCTGTTTATGCCAGCGACTACCGTCGAGTTCATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACC  >  NC_000913/3832509‑3833100
                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                            
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)
                                                                                                                                                                                                                                                                                                   cGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCa    >  2:610819/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                     cACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCAACGACCAGGCTTTCCGCCATTGATTCGTAGTTAACGTAACAACCCCCaac  >  1:227685/1‑299 (MQ=255)
                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                            
GTTGATCGTTGTTGCCCAGCGCCAGTAGAGTGT‑TATCACGCACATAAACGGGCAGACTCAGGAAGCCGTGCTGCTGTTTATGCCAGCGACTACCGTCGAGTTCATCGTTGTGCCACAGGTGTGTCCAGCGACCTTCCGGCAGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCAT‑GGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACC  >  NC_000913/3832509‑3833100

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: ‑