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.1% 362.9 / NA 107K392K (AAG→AAApdeKc‑di‑GMP phosphodiesterase PdeK
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (56/49);  minor base G (0/2);  total (56/51)
Fisher's exact test for biased strand distribution p-value = 2.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 107 total aligned reads displayed.

TAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCA  >  NC_000913/3684199‑3684687
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
tAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCt                                                                                                                                                                                                                                                >  1:559772/1‑251 (MQ=255)
       cAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCTCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGc                                                                                                                                                                                                                                         <  1:368498/251‑1 (MQ=255)
            cacaGCACGCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCg                                                                                                                                                                                                                                    <  2:219909/251‑1 (MQ=255)
               aGCACTCAATGCGCTCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAACGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGt                                                                                                                                                                                                                                  >  2:39951/1‑250 (MQ=255)
               aGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGt                                                                                                                                                                                                                                  >  1:221441/1‑250 (MQ=255)
                  ccTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAATCTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCGGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCTTTTATGGCGAGCGGTAAAt                                                                                                                                                                                                                              <  2:42747/250‑1 (MQ=255)
                     cAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGc                                                                                                                                                                                                                            <  1:140612/250‑1 (MQ=255)
                     cAATGCGATCGATTAAGCCATCCGGCAGGGCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGc                                                                                                                                                                                                                            <  1:215148/250‑1 (MQ=255)
                      aaTGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCg                                                                                                                                                                                                                          <  1:220210/251‑1 (MQ=255)
                          cgatcgatTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGc                                                                                                                                                                                                                         >  2:447777/1‑248 (MQ=255)
                              cgatTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAAt                                                                                                                                                                                                                   <  1:456130/250‑1 (MQ=255)
                              cgatTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGAGCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATa                                                                                                                                                                                                                  <  2:347691/251‑1 (MQ=255)
                                atTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAgc                                                                                                                                                                                                                <  2:556279/251‑1 (MQ=255)
                                     gCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGAc                                                                                                                                                                                                           <  1:546687/251‑1 (MQ=255)
                                     gCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGAc                                                                                                                                                                                                           <  1:330768/251‑1 (MQ=255)
                                      acATCAGGCAGGGCCCAACTGCCATCCGGGTGCTGGATACGCAGTAACACTTCCGCACTGATCAGTTTACAGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGATTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTCTGGCGAGCGGTAAATGCCGCAGAACTAGcgcgcct                                                                                                                                                                                                          >  2:155234/2‑248 (MQ=255)
                                           cggaaGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTGACACTTACGAACTGACCAGTTTACAGCTAGTCATCTCGACCTGTGGCTGTCACCAAATAGAAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACGGCTTCTGGGCGGCTTCCATCTGCTGCGGATCCAAGAAATGAATCTGATTTTTGCCTTTATGGCGCGCGGTCAATGCCGCAGACATAGCGCGCCTGTaa                                                                                                                                                                                                      >  1:33348/5‑250 (MQ=255)
                                            ggCAGGTCCCACCTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCACAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTaa                                                                                                                                                                                                      >  1:482452/1‑249 (MQ=255)
                                               aGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACCCTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTg                                                                                                                                                                                                 <  1:426348/251‑1 (MQ=255)
                                               aGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGcc                                                                                                                                                                                                                           >  2:44429/1‑225 (MQ=255)
                                               aGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGcc                                                                                                                                                                                                                           <  1:44429/225‑1 (MQ=255)
                                                 gTCCCAACGGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTGGTGGCTTTATGGCGAGCTGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGtt                                                                                                                                                                                               <  2:389209/251‑1 (MQ=255)
                                                  tCCCAACTGCCTTCCGGTTGCTGGATAAGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTAGACCTGTGGCGGTAACAAAATAGCAAACTGATGATTTTCCAGTGAATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCAATATGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAACAGCCGGCGAACTAACGCGACTGTAAATCTGttg                                                                                                                                                                                              >  2:28346/1‑250 (MQ=255)
                                                         tGCCATACGTTTGGTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGATGGTCATATCGCCCAGTGGCTGTACCCAAATAGCAAACTGATGATTTTCAAGTGCATTAAGGATATCGCTATCTTCCGTCAACAGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGCTTTTTGCCTTTATGGCGAGCGGTCAATGCAGCAGAAATAGCACGACTGTAAAACTGTTctgcgg                                                                                                                                                                                         >  2:389256/1‑249 (MQ=255)
                                                          gCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTg                                                                                                                                                                                       <  1:100889/250‑1 (MQ=255)
                                                           ccATCCGGTTGCTGGATACGCAGGAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTAGACCTGTGGCGGTAACCAAATAGCCAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATAAAAGAACTGCATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGAGAGACGGTATCGCTGTTCGGCGGTg                                                                                                                                                                                       >  2:499325/1‑249 (MQ=255)
                                                             aTCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGAt                                                                                                                                                                                   <  2:226004/251‑1 (MQ=255)
                                                               ccGGTTGCTGGCTACGCCGTAACACTTCCGCACTGACCCGTTTCCCGCTGGTCATCTCGCCCTGTGGCTGTAACCAAATCGAAAACTGATGATTTTCCCGTGCCTTCAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGATGCGGATCAAAGAACTGCATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCAGCAGAAATACCGCGACTGTTAAGCTGTTCGGCGGTGAGATcc                                                                                                                                                                                 >  2:82929/1‑250 (MQ=255)
                                                                  gTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCg                                                                                                                                                                              <  2:221441/251‑1 (MQ=255)
                                                                   ttGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAACTAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCgg                                                                                                                                                                             >  1:530464/1‑250 (MQ=255)
                                                                        ggATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAAc                                                                                                                                                                        <  1:26085/251‑1 (MQ=255)
                                                                            aCGCAGTACCACGTCCGCACTGTCCAGTTTACCGCTTGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTTCATTAAGGATATCTCTCTCTTCCGTCAACGGCTTCTGGGCGGCTTCCATCTGCGGCGGATCAAAGATCTGAATCGGATTTTTGCCTTTATGGCGGGCGGTAAATGCCGCATAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATc                                                                                                                                                                     <  1:378199/250‑1 (MQ=255)
                                                                            aCGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATcg                                                                                                                                                                    >  1:392471/1‑251 (MQ=255)
                                                                               cAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATcgcca                                                                                                                                                                 <  1:503570/251‑1 (MQ=255)
                                                                                aGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATcgccac                                                                                                                                                                <  1:385432/251‑1 (MQ=255)
                                                                                 gTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTTGGCGACTTCCCTCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAACTAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAAAACATcgccacg                                                                                                                                                               >  1:516895/1‑251 (MQ=255)
                                                                                   aaCACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGTCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGAGTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATcgccacgcc                                                                                                                                                             <  2:238874/251‑1 (MQ=255)
                                                                                              aCTGACCAGTTTACCGCTGGTCCTCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGt                                                                                                                                                  >  1:187959/1‑251 (MQ=255)
                                                                                                 gACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGcg                                                                                                                                               >  1:518842/1‑251 (MQ=255)
                                                                                                  aCCAGTTTACCGCTGGTCGTCTCGACCTGTGGCTGTAACCAAGTAGCAAACTGATGCTTTTCCAGTGCATTAAGGATATCGCCCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCCGAAATAGCGCGCCTGTAAAGCTGTTCGGCGGTGAGATAGCCGGAGCACATTGCAACGCCAATGCTACAGAGcgg                                                                                                                                              >  2:541841/1‑251 (MQ=255)
                                                                                                  aCCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGcgg                                                                                                                                              >  1:101152/1‑251 (MQ=255)
                                                                                                  aCCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGcgg                                                                                                                                              <  1:31866/251‑1 (MQ=255)
                                                                                                     aGTTTACCGCTGGTCATCTCGCCCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGTGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGcggacg                                                                                                                                           <  2:384967/251‑1 (MQ=255)
                                                                                                     aGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGcggacg                                                                                                                                           >  1:443410/1‑251 (MQ=255)
                                                                                                          aCCGCGGGTCATCACGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTGCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACAGCGCCACGCCAATGCTACAGTGCGGACGGAGtt                                                                                                                                      <  1:481522/251‑1 (MQ=255)
                                                                                                            cGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGCAGAACGTCGCCACGCCAATGCTACAGTGCGGGCGGAGTTg                                                                                                                                     >  2:330438/1‑250 (MQ=255)
                                                                                                             gCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGaa                                                                                                                                   <  1:488175/251‑1 (MQ=255)
                                                                                                             ccgatttaTCTCGACCTGTGGCGGGAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCGCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGAGTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGAGCGCCGTAGAACATCACCACGCCAAGGCTACAGTGCGGACGGAGTTGa                                                                                                                                    <  2:102066/243‑1 (MQ=255)
                                                                                                                 gTCATCTCGACCTGTGGCTGTAACCAAATAGAAAACTGATGATTTTCCAGTTCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCCAAGAACTGAATCTGATATTTGCCTTTATGGCGAGCGGTAAATGCCGCAGCAATAGCTCGACTGTAAAGCTGTTAGGCGGTGAGTTCGCCGCAGAACATCGCCACGACAATGCTACAGTGCGGACGGAGTTGATTACg                                                                                                                               >  1:490159/1‑251 (MQ=255)
                                                                                                                  tCATCTCGACCTGTGGCTGTAACCAAATAGCAACCTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGCAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCACTGCTACAGTGCGGACGGAGTTGAATACGt                                                                                                                              >  2:213512/1‑251 (MQ=255)
                                                                                                                    aTCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTCGAACATCGCCCCGCCAATGCTACAGTGCGGACGGAGTTGacgatcgg                                                                                                                             >  1:102066/1‑243 (MQ=255)
                                                                                                                       tcGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAAt                                                                                                                         >  1:190446/1‑251 (MQ=255)
                                                                                                                         gACCTGTGGCTGTAACCAAATAGCAAAATGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGCGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCg                                                                                                                       >  2:517396/1‑251 (MQ=255)
                                                                                                                           ccTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGc                                                                                                                     >  1:119797/1‑251 (MQ=255)
                                                                                                                              gtgGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTGGAATACGTTCAATCGGCa                                                                                                                    <  2:467538/249‑1 (MQ=255)
                                                                                                                              gtgGCTGTAACCAAAGAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTTAATACGTTCAATCGGCAgg                                                                                                                  <  1:512343/251‑1 (MQ=255)
                                                                                                                               tgGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGgc                                                                                                                 >  1:294204/1‑251 (MQ=255)
                                                                                                                                ggatGTAACCAAATAGCAAACTGATGCTTTTGCAGTGCATTAATGATATCGATCTCTTCCGTCTTCCGCTTCTGGGTGGATTCCATCTGCTGCGGATCAAAGAACTTAAGCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCTCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGGGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGAGGACGGTGTTGTAGACGTTCAATCGGCATgcg                                                                                                                <  2:161499/248‑1 (MQ=255)
                                                                                                                                ggCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGgcg                                                                                                                >  2:572720/1‑251 (MQ=255)
                                                                                                                                       aCCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCgctcgctc                                                                                                         <  2:294204/251‑1 (MQ=255)
                                                                                                                                       aCCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGGGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTGCAATCGGCAGGCgctcgctc                                                                                                         <  1:501654/251‑1 (MQ=255)
                                                                                                                                        ccAACTAGCAACCTGTTGATTTTCCAGTGCATTAAGGATATCGCTATCTTCAGACAACCGCTTCTGGGCCGCTTCCCTCAGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTCTATGGCGCGCGGTAAATGCCGCAGAAATAGCGCGACTGTCAAGCTGTTCGGCGGTGAGATCGCCGTAGACCATCGCCACGCCACTGCTACAGTGCGGACGGCGTTGAATACGTTCAATCGGCAGgcgaccgctcc                                                                                                        >  2:582736/1‑243 (MQ=255)
                                                                                                                                           aaTAGCAAACTGATTATTTTCCAGTGCATGAACGATATCGCTCTCTTCCGTCAACCGCTTCTTGTCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTGTGCCTTTATTGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTATAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGa                                                                                                     <  1:515663/251‑1 (MQ=255)
                                                                                                                                              aacAAACTGATGATTTTCCAGTGAATTAAGGCAATCGCTCTCTACCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGCACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAGAGCTGTTCGGCGGTGAGATCGCAGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGATGAATACGTTCAAACGGCAGGCGCACGCTCATGAt                                                                                                    >  2:336652/3‑249 (MQ=255)
                                                                                                                                                cAAACTGTTGATTTCCCAGCGCATTATGGATATCGCTCTCTTCCTTCAACCGCTTCTGGGCGTCTTCCATCTGCTGCGGATCAAAGAACTGAATCGGATTGTTGCCGTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTCGAACATCGCCACGCCAATGCTACATGGCGTACGGAGTTGAATACGGTCAATCGGCAGGCGCTCGCTCATGATAGTg                                                                                                <  1:499325/251‑1 (MQ=255)
                                                                                                                                                cAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTg                                                                                                <  2:570223/251‑1 (MQ=255)
                                                                                                                                                 aaaCTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTg                                                                                                >  1:494980/1‑250 (MQ=255)
                                                                                                                                                      gatgatTTTACAGTGCATTCAGGATATCGCTCTATTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACt                                                                                          >  1:449653/1‑251 (MQ=255)
                                                                                                                                                      gatgatTTGCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACt                                                                                          >  2:343108/1‑251 (MQ=255)
                                                                                                                                                       atgatTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGCGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACtt                                                                                         >  1:106370/1‑251 (MQ=255)
                                                                                                                                                                cAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCAATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGACGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGCTCGTGAGCACTTGCTGACCt                                                                                 >  2:300673/1‑250 (MQ=255)
                                                                                                                                                                          aggattgcgctcgctTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTTCCGTTAGGGCTAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTGCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGc                                                                      <  2:558268/244‑1 (MQ=255)
                                                                                                                                                                             atatCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTg                                                                   >  1:53146/1‑251 (MQ=255)
                                                                                                                                                                                tCGCTCTCTTCCGTAAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGCATCTGCTTTTTGCCTTTATGGCGAGCGGTAACTGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCCTGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGcc                                                                 >  1:91295/1‑250 (MQ=255)
                                                                                                                                                                                    tctctTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGAGATTGCGTGCCACGGt                                                            >  1:255817/1‑251 (MQ=255)
                                                                                                                                                                                       ctTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTcc                                                         >  1:345444/1‑251 (MQ=255)
                                                                                                                                                                                            gTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCACTGCTACAGTACGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTac                                                    >  1:347525/1‑251 (MQ=255)
                                                                                                                                                                                              cAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGAACTAAGGTGATTGCGTGCCACGGTTCCTGTacac                                                  >  1:199207/1‑251 (MQ=255)
                                                                                                                                                                                               aaCCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACAcc                                                 <  1:561128/251‑1 (MQ=255)
                                                                                                                                                                                                  cgcatCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGAGTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATATTGAGCACTTGCTGACCTAAGGGGATTGCGTTCCACGGTTCCTGTACACCGtt                                              <  2:263884/247‑1 (MQ=255)
                                                                                                                                                                                                      tCTGGTCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTTATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCGGTTCGGCGGTGGGATCGCCGTAGAACATCGCCACGCTAATTCTACAGTGCGGACGGAGTTGAATACTTTCAATCGGCATGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTgg                                            <  1:92453/249‑1 (MQ=255)
                                                                                                                                                                                                      tCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGc                                           <  1:364314/250‑1 (MQ=255)
                                                                                                                                                                                                         gggCGGCTTCCTTCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCGCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGACACGGTTCCTGTACACCGTTGTCaa                                         <  2:137256/249‑1 (MQ=255)
                                                                                                                                                                                                         gggCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATg                                       <  1:193148/251‑1 (MQ=255)
                                                                                                                                                                                                         gggCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTGATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGAAAAGCGGTCCGGCGGTGAGAGCGCCGTAGAACATCGTCACGCCAATGCTACAGTGCGGACGGAGGTGAATACGTGCAATCGGCAGGCGCGCGCGCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCGGTCCACCGGTggtaa                                         <  2:164219/249‑4 (MQ=255)
                                                                                                                                                                                                            cggcTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGAc                                     >  1:86167/1‑250 (MQ=255)
                                                                                                                                                                                                              gcTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTGTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGc                                  <  2:220557/251‑1 (MQ=255)
                                                                                                                                                                                                                ttCCATCTGCTGCGGATCAAAGAACTGAATCTGCTTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCa                                 >  1:480457/1‑250 (MQ=255)
                                                                                                                                                                                                                   cATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCAAGGTTCCTGTACACCGTTGGCAATGACAGCAAAGt                             >  1:172244/1‑251 (MQ=255)
                                                                                                                                                                                                                     tCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCAACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATCGTGAGCACTTGCTGACCTAAGATGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCaaattc                            >  1:566162/1‑247 (MQ=255)
                                                                                                                                                                                                                     tCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACCGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTc                            >  1:778/1‑250 (MQ=255)
                                                                                                                                                                                                                     tCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACCGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTc                            <  2:778/250‑1 (MQ=255)
                                                                                                                                                                                                                        gctgcGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGGTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCGCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTc                            <  2:269589/247‑1 (MQ=255)
                                                                                                                                                                                                                                  caaGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGACATAGCGAGACTGTAATGCTGTTCGGAGGTGAGATCGACGTAGAAAATAGCCACGCCAATGCTACAGCGCGGGCGGGGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACATAAGGGGATTGCGAGCCACGGTTCCTGTACACCGTTGGCAAATAAAGCAAAGTCATGACCACTAATCt              >  2:477659/2‑251 (MQ=255)
                                                                                                                                                                                                                                  aaaGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGCGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCt              >  2:551547/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGAGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCGAATCTGCGCGAgg      >  1:197077/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gAATCTGATTTTTGCCGTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAg       <  2:479429/250‑1 (MQ=255)
                                                                                                                                                                                                                                             tCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATc   >  1:465905/1‑251 (MQ=255)
                                                                                                                                                                                                                                              cTGATTTTTGCCTTTGTGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTACAGCTGTTCGGCGGAGAGATCGCCGTAGACCATCGCCACGCCAATGCTAAAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGAGAGTGAGCACTTGCTGGCCTAAGGGGATTGCGGGCCACGGTTCCTGAACCCCGTTGGCCATGACAGCAAAGGCATAACCGCTAATCTGCGCGAGGagc   >  2:212598/1‑248 (MQ=255)
                                                                                                                                                                                                                                              cTGATTTTTGCATTTATGGCGAGAGGTAAATGCCGCAGAAATAGCGAGACTGTAAAGCTGTTCGGCGGTGAGATCGACGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATAAGTTCAATCGGCAGGCGCTCGCTCATGATAGCGAGCACTTGCTGACCTAAGGTGGTTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTAATAAACGCTAATATGCGACAGGATct  >  2:223202/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
TAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCA  >  NC_000913/3684199‑3684687

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