Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,684,443 C→T K392K (AAG→AAA pdeK ← c‑di‑GMP phosphodiesterase PdeK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,684,4430CT98.9% 292.5 / NA 87K392K (AAG→AAApdeKc‑di‑GMP phosphodiesterase PdeK
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (43/43);  total (44/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 = 9.98e-01

AACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCAT  >  NC_000913/3684200‑3684688
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
aaCCATCAGCCCACAGCACTCAATGAGCTCGCTTAAGCCCTCCGGCAGGCCCCCACAGCCATCCGGTTGCTGGCTACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATCTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGAGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCtt                                                                                                                                                                                                                                                >  2:298314/1‑251 (MQ=255)
 aCCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCTGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTTACCAGTTTACCGCTGGTCATCTCGACCTGTTGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGTCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCttt                                                                                                                                                                                                                                               <  1:350434/251‑1 (MQ=255)
    aTCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACCCTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATg                                                                                                                                                                                                                                            <  1:478468/251‑1 (MQ=255)
    aTCAGCCCACAGCACACACTGCACTCGATTACGCAATCCGGCGGGTCCCAACTGCCCTAAGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAAAAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGCACTGATACTGATTTTTGCCTTTATg                                                                                                                                                                                                                                            >  2:372677/1‑251 (MQ=255)
       aGCCCACAGCACTCAATGCGATCGATTAAGCCCTCCGGCAGGTCACAACTGCCATCCGGTTGCTGGATACGCAGTCACACTTCAGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCg                                                                                                                                                                                                                                         >  1:503583/1‑251 (MQ=255)
       aGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCg                                                                                                                                                                                                                                         <  1:310054/251‑1 (MQ=255)
           cacaGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCg                                                                                                                                                                                                                                     <  1:24040/251‑1 (MQ=255)
               gCACTCAATGCGATCGATTAAGCCATCCGGCAGGGCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACTCTTCCGCAATGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAAGCTGACTTTTGCCTTTATGGCGAGCTGTaa                                                                                                                                                                                                                                 <  2:227174/251‑1 (MQ=255)
                        gCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAg                                                                                                                                                                                                                        <  2:317656/251‑1 (MQ=255)
                          gatcgttTCAGCCATCCGGCCGGTACCAGCTGCCCTCCGGTTTCGCGATATGACGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGATCTGTGGCTGTATCCAAATAGCAAACTGATGATTTTACAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGAATTTTGCCTTTTTGGCGAGCGGTAAATGCTGCAGa                                                                                                                                                                                                                       >  2:264336/1‑250 (MQ=255)
                          gatcgatTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGaa                                                                                                                                                                                                                      <  2:514379/251‑1 (MQ=255)
                            tcgatTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGCTTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGaaa                                                                                                                                                                                                                     >  2:377020/1‑250 (MQ=255)
                              gatTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGATCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGaatt                                                                                                                                                                                                                    >  2:309480/1‑247 (MQ=255)
                               atTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGCGCGGTAATTGCCGCAGAAATAgc                                                                                                                                                                                                                 >  2:296490/1‑251 (MQ=255)
                                 tAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAgcg                                                                                                                                                                                                                <  1:180810/250‑1 (MQ=255)
                                 tAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAgcg                                                                                                                                                                                                                <  1:363137/250‑1 (MQ=255)
                                     ccATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACt                                                                                                                                                                                                           >  1:78334/1‑251 (MQ=255)
                                      cATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTg                                                                                                                                                                                                          >  2:350031/1‑251 (MQ=255)
                                           ggCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTaa                                                                                                                                                                                                       >  1:561399/1‑249 (MQ=255)
                                            gCAGTTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACCCTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGTCTGTAACCAAATAGCAACCTGATGATTTTCCAGTGCATTAATGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCTGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCTGTAAATCCCGCAGAAATAGCGCGACTGTAAAGc                                                                                                                                                                                                    <  1:582120/251‑1 (MQ=255)
                                               ggTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACCCTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGt                                                                                                                                                                                                 <  2:200669/251‑1 (MQ=255)
                                                       cTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTcggcggg                                                                                                                                                                                         >  2:107563/1‑250 (MQ=255)
                                                               cGGTTGCTGGATAAGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTACAGTGCATTAAGGCTATCGCTATCTTCCGTCAACCGCTTCTGGGCGGCTTCAATCAGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATAGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGc                                                                                                                                                                                 >  1:571144/1‑251 (MQ=255)
                                                                ggTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGcc                                                                                                                                                                                <  1:262205/251‑1 (MQ=255)
                                                                      tGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGGAATAGCGCGACTGTAACGCTGCTAGGCGGTGAGATCGCCGTAGa                                                                                                                                                                           >  1:392636/1‑250 (MQ=255)
                                                                        gATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTCCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATTGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCTGTGAGATCTCCGTAGAACa                                                                                                                                                                        <  2:499618/251‑1 (MQ=255)
                                                                           aCGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATcg                                                                                                                                                                     <  2:416968/251‑1 (MQ=255)
                                                                           aCGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATc                                                                                                                                                                      <  2:513049/250‑1 (MQ=255)
                                                                             gCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGCAATAGCGCGACTGTAAAGCTGTTc                                                                                                                                                                                               <  1:187634/223‑1 (MQ=255)
                                                                             gCAGTAACACTTCCGCAATGACCAGTTTACCGCTGGTCATCTCGACCTATGTCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGACAACCGCTTCTGGGCGGCTTCCATCTGCTGCGCATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGCGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTc                                                                                                                                                                                               >  2:187634/1‑223 (MQ=255)
                                                                                  aaCACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCGTCTGGGCGGCTTCCATCTGCTGCGGATCAAGGAACGGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGGGCGACTGTAAAGCTGGTCGGCGGTGAGATCGCCGTAGAACATcgccacgcc                                                                                                                                                              <  1:34729/251‑1 (MQ=255)
                                                                                       ttCCGCAATGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAACGCTGTTCGGCGGGGAGATCGCCGTAGCACATCGCCACGCCAATGc                                                                                                                                                         >  2:567824/1‑251 (MQ=255)
                                                                                        tCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCt                                                                                                                                                        >  1:148805/1‑251 (MQ=255)
                                                                                         ccGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGGGCATTACGGATATCGCTCTCTTCCGTCAACCGCTTCTGGTCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATAGCGCGCGGAAAATGCCCCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGATGAGATCGCCGTAGAActccccaacgccaatactt                                                                                                                                                       >  1:432615/1‑233 (MQ=255)
                                                                                          cGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTAc                                                                                                                                                      >  1:100491/1‑251 (MQ=255)
                                                                                                    aGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGcggacg                                                                                                                                            <  2:138169/251‑1 (MQ=255)
                                                                                                     gTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCTCGCCAATGCTACAGTGcggacgg                                                                                                                                           <  1:122180/251‑1 (MQ=255)
                                                                                                                gTCATCTCGACCGGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATGAATGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAGTAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTGGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACg                                                                                                                                <  1:82189/251‑1 (MQ=255)
                                                                                                                  cATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGtt                                                                                                                              >  2:227896/1‑251 (MQ=255)
                                                                                                                       cGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCaa                                                                                                                           <  2:235025/249‑1 (MQ=255)
                                                                                                                       cGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATc                                                                                                                         <  2:511229/251‑1 (MQ=255)
                                                                                                                        gACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGc                                                                                                                                                                                                                             <  2:30599/150‑1 (MQ=255)
                                                                                                                        gACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGc                                                                                                                                                                                                                             >  1:30599/1‑150 (MQ=255)
                                                                                                                        gACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCg                                                                                                                        <  1:458634/251‑1 (MQ=255)
                                                                                                                            ggtgGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCACCCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAg                                                                                                                    <  2:40155/250‑1 (MQ=255)
                                                                                                                               ggCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGgc                                                                                                                  <  1:158376/250‑1 (MQ=255)
                                                                                                                               ggCTGTAACCAAATAGCAAACGGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGGCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGgcg                                                                                                                 <  2:184533/251‑1 (MQ=255)
                                                                                                                                      aCCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCACTCGGCAGGCgctcgctc                                                                                                          >  1:495349/1‑251 (MQ=255)
                                                                                                                                           aTAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTCCAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGAt                                                                                                     >  1:525525/1‑251 (MQ=255)
                                                                                                                                           aTAGCAAACTGATGATTTTACAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCGCGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGAt                                                                                                     >  1:461768/1‑251 (MQ=255)
                                                                                                                                             aGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGCACATCGGCACGACAATGCTACAGTGAGGACGGAGTTGAATACGTTCAATAGACAGGCACTCGCTCATGAt                                                                                                     >  1:258789/1‑249 (MQ=255)
                                                                                                                                             aGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATTCGTTCAATCGGCAGGCGCTCGCTCATGATAg                                                                                                   <  1:197508/251‑1 (MQ=255)
                                                                                                                                             aGCAAACTGATGATTTTCCAGGGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCTCTCGCTCATGATAg                                                                                                   <  1:469527/251‑1 (MQ=255)
                                                                                                                                                 aaCTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTg                                                                                                 >  1:438624/1‑249 (MQ=255)
                                                                                                                                                   cTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGc                                                                                              >  1:104430/1‑250 (MQ=255)
                                                                                                                                                     gatgatcttcaaccgCATTAACGATATCGCTCTCTTCCGTCACCCGCTTCTGGGCGGCTTCACTCTGCTGCGGATCACAGCACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTACCTGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTCGAACATCGCCACGCCAATGCTACAGTGCGGAAGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACt                                                                                           >  2:276487/15‑251 (MQ=255)
                                                                                                                                                     gatgatTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCATCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATCGTGAGCACt                                                                                           >  1:45507/1‑251 (MQ=255)
                                                                                                                                                         atTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCt                                                                                       <  1:377020/251‑1 (MQ=255)
                                                                                                                                                          ttttCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTg                                                                                      >  1:110855/1‑251 (MQ=255)
                                                                                                                                                             tCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGAcc                                                                                   >  2:420067/1‑251 (MQ=255)
                                                                                                                                                                  tGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAgg                                                                              >  2:495101/1‑251 (MQ=255)
                                                                                                                                                                   gCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAgg                                                                              <  1:243380/250‑1 (MQ=255)
                                                                                                                                                                    cATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTCACCTAAGGTg                                                                            >  2:172683/1‑251 (MQ=255)
                                                                                                                                                                     aTTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGa                                                                           <  1:13465/251‑1 (MQ=255)
                                                                                                                                                                         aGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTgg                                                                       >  2:422273/1‑250 (MQ=255)
                                                                                                                                                                          ggATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCg                                                                      >  1:23142/1‑251 (MQ=255)
                                                                                                                                                                           gATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGt                                                                     >  1:386870/1‑251 (MQ=255)
                                                                                                                                                                                  ctctctTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACgg                                                              >  1:450905/1‑251 (MQ=255)
                                                                                                                                                                                  ctctctTCCGTAAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGCAATGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTACAGCTGTTCGGCGGTGAGATCGCCGTAGCACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGCGAGCACTTGCTGCCCTAAGGTGATTGCGTGCCACgg                                                              >  2:251932/1‑251 (MQ=255)
                                                                                                                                                                                      ctTCCGTCAAACGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAACGAACTAACTCTGATTTTTGCCTTTATGGCGGGCGGTAAATGCCGCAGAACAAGCGCGACTGTAACGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGTTACCGGGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATAATAGTGAGAATTTGCTGACATAAGGGGATTGCGGGCCACGGTTcc                                                          >  2:83051/1‑251 (MQ=255)
                                                                                                                                                                                         ccGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTg                                                        <  1:350031/250‑1 (MQ=255)
                                                                                                                                                                                                       tGGGCGGCTTCCATCTGCTGCGGCTCACCGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGGCGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTAACCCGTTGGCAAt                                         >  2:60759/1‑251 (MQ=255)
                                                                                                                                                                                                        gggCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCa                                           <  1:399609/248‑1 (MQ=255)
                                                                                                                                                                                                        gggCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATg                                        <  1:161693/251‑1 (MQ=255)
                                                                                                                                                                                                          gcggcTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGAc                                      <  1:358449/251‑1 (MQ=255)
                                                                                                                                                                                                           cggcTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACa                                     <  1:29089/251‑1 (MQ=255)
                                                                                                                                                                                                           cggcTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACa                                     <  1:39702/251‑1 (MQ=255)
                                                                                                                                                                                                                   aTCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTg                                                                                         >  1:116901/1‑191 (MQ=255)
                                                                                                                                                                                                                   aTCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTg                                                                                         <  2:116901/191‑1 (MQ=255)
                                                                                                                                                                                                                         tgcGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCCGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATaaac                       >  2:285628/1‑249 (MQ=255)
                                                                                                                                                                                                                             gATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCt                    >  1:218899/1‑250 (MQ=255)
                                                                                                                                                                                                                                  aaGAACTGAATCTGATTTTTGCCTTTATGGCGGGCGGTAAATGCCGCAAAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTg              <  2:448966/251‑1 (MQ=255)
                                                                                                                                                                                                                                     aaCTGAATCTGCTTTTTGCCTTTATGGCGCGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTgcgc           >  1:419673/1‑251 (MQ=255)
                                                                                                                                                                                                                                         gAATCTGATTTTTGCCTTTCTGCCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATAGCCGTAGAACATCGCCACGCCAATGCTAAAGTGCGGAGGGAGTTGGATACGTTCAATCGGCAGGCGCTAGCCCATGATAGTGAGCACTTGCTGACCTCAGGTGATTGCGTGCCACGGTTCCTGTACACCCTTGGCAATGAAAGCAAAGGCATAACCGCTAATCTGCGCGAgg       >  1:226588/1‑251 (MQ=255)
                                                                                                                                                                                                                                         gAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAAGGACAGCAAAGTCATAACCGCTAATCTGCGCGAg        <  1:181949/250‑1 (MQ=255)
                                                                                                                                                                                                                                         gAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCGGTTCGGCGGTGAGATCGCCGTAGAACAGCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAg        <  2:229732/250‑1 (MQ=255)
                                                                                                                                                                                                                                              ggATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGGTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCAt  <  1:449532/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
AACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCAT  >  NC_000913/3684200‑3684688

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