Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 961,699 A→T 7.3% L239Q (CTG→CAG)  selD ← Selenide, water dikinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451961,6990AT 7.3% 237.3 / 6.0 68L239Q (CTG→CAG) selDSelenide, water dikinase
Reads supporting (aligned to +/- strand):  ref base A (36/26);  new base T (3/2);  total (39/29)
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

CCTCCGGGGCAACAGCTACCAGCAAGCCACCGCTGGTCTGCGGGTCGCACAGCAGGTGCTTCTGGTCGTCGCTCAGGGCGCTGATCTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAGTGGTGAGGATGCCGATGCCCAAGGGCTTGGTCAGGTACAGCTGGCAGCCTGCGCTGGCGGTG  >  AE015451/961481‑961920
                                                                                                                                                                                                                          |                                                                                                                                                                                                                             
ccTCCGGGGCATCAGCTACCAGCAAGCCACCGCTGGTCTGCGGGTCGCACAGCAGGTGCTTCTGGTCGTCGCTCAGGGCGCTGATCTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCGCGACCAGGTGGcc                                                                                                                                                                                                               >  1:1052508/1‑235 (MQ=255)
                        aGCCACCGCTGGTCTGCGGGTCGCACAGCAGGTGCTTCTGGTCGTCGCTCAGGGCGCTGATCTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCg                                                                                                                                                                                      <  1:328230/236‑1 (MQ=255)
                                                       gTGCTTCTGGTCGTCGCTCAGGGCGCTGATCTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGt                                                                                                                                                                                    >  1:253534/1‑207 (MQ=255)
                                                             cTGGTCGTCGCTCAGGGCGCTGATCTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCgg                                                                                                                                                  >  1:688898/1‑235 (MQ=255)
                                                               ggTCGTCGCTCAGGGCGCTGATCTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTg                                                                                                                                               >  1:12629/1‑236 (MQ=255)
                                                                    tcgCTCAGGGCGCTGATCTTGGGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGaa                                                                                                                                                      <  1:1008838/224‑1 (MQ=255)
                                                                      gCTCAGGGCGCTGATCTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGcc                                                                                                                                             <  2:564143/231‑1 (MQ=255)
                                                                        tCAGGGCGCTGATCTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGcc                                                                                                                                             <  1:630138/229‑1 (MQ=255)
                                                                        tCAGGGCGCTGATCTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGcc                                                                                                                                             >  2:630133/1‑229 (MQ=255)
                                                                            ggCGCTGATCTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAg                                                                                                                                                                                                                              <  2:574257/144‑1 (MQ=255)
                                                                              ggcTGATCTTGTGCCCGTAGCTGTCGTAGTTACGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGc                                                                                                                                                                                                                >  2:1183417/2‑156 (MQ=255)
                                                                                     cTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGcc                                                                                                                                             <  1:141376/216‑1 (MQ=255)
                                                                                        gtgCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGa                                                                                                                                                                                           <  1:1038795/167‑1 (MQ=255)
                                                                                        gtgCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGcc                                                                                                                                             <  1:931227/213‑1 (MQ=255)
                                                                                              gTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTcgc                                                                                                                >  1:1262498/1‑236 (MQ=255)
                                                                                                aGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAGGCCgggg                                                                                                                                                                                            >  2:213806/1‑156 (MQ=255)
                                                                                                aGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTg                                                                                                                                               <  1:912405/203‑1 (MQ=255)
                                                                                                aGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTg                                                                                                                                               >  2:912400/1‑203 (MQ=255)
                                                                                                aGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTg                                                                                                                             >  2:10261/1‑221 (MQ=255)
                                                                                                   tGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGaacc                                                                                                                                                    >  1:133352/1‑193 (MQ=255)
                                                                                                        tAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAgg                                                                                                        <  1:1244044/234‑1 (MQ=255)
                                                                                                        tAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGc                                                                                                       >  1:681136/1‑235 (MQ=255)
                                                                                                          tttGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGa                                                                                                                                                                                           <  2:420096/148‑1 (MQ=255)
                                                                                                            ggcgcAGGGTTCCACAGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCg                                                                                                                                                                                      <  1:333467/151‑1 (MQ=255)
                                                                                                               gcAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCgtgt                                                                                                                                     <  1:145272/198‑1 (MQ=255)
                                                                                                                cAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTcc                                                                                               >  1:1187374/1‑235 (MQ=255)
                                                                                                                  gggTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTg                                                                                             >  1:1071982/1‑235 (MQ=255)
                                                                                                                      tCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGc                                                                                                                                                 >  2:848402/1‑179 (MQ=255)
                                                                                                                       ccACCGGGGATGCAGCGTTCGGCCAGGTAATGGTCGACACTAGGAAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGAcc                                                                                                                                                                                                                       <  2:406055/108‑1 (MQ=255)
                                                                                                                                  gCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCgt                                                                                                                                       >  1:1281710/1‑177 (MQ=255)
                                                                                                                                   cAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAg                                                                                                                    <  1:524099/195‑1 (MQ=255)
                                                                                                                                            gCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGCGCAGGCCGAAGCCggg                                                                                                                                                                                             >  2:192716/1‑112 (MQ=255)
                                                                                                                                               aGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTAGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAgt                                                               <  1:959838/236‑1 (MQ=255)
                                                                                                                                                ggTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCC‑‑‑‑‑‑ACC‑‑‑TGGTCGAGCTGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCgtgt                                                                                                                                     <  1:869889/156‑1 (MQ=255)
                                                                                                                                                   aaTGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGGGGGTGCACAGCCAGTcg                                                                                                                 >  1:353697/1‑182 (MQ=255)
                                                                                                                                                   aaTGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTc                                                                                                                  >  2:1076775/1‑181 (MQ=255)
                                                                                                                                                     tGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCg                                                                    <  1:732776/225‑1 (MQ=255)
                                                                                                                                                          gACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCATGACGAAGCCGGTGACGTCGGTCATGGCCTTGACGCCGt                                                                                                                                                                   >  2:702252/1‑125 (MQ=255)
                                                                                                                                                            cacTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGGTTGCCGAAGCgg                                                                                                                                                  >  2:869813/1‑140 (MQ=255)
                                                                                                                                                             acTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGCTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGcac                                                                                          <  1:43222/195‑1 (MQ=255)
                                                                                                                                                              cTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGGGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCCTCTTTTCGGCa                                                                 >  1:136859/1‑219 (MQ=255)
                                                                                                                                                               tAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAGTGGTGAGgatgccgat                                                >  1:919602/1‑235 (MQ=255)
                                                                                                                                                               tAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGGGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGcc                                                                                                                                                                                                               <  2:25723/76‑1 (MQ=255)
                                                                                                                                                                  gcaggcGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGGGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGcc                                                                 >  1:979098/1‑214 (MQ=255)
                                                                                                                                                                      gcGTGGCATGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCgtg                                                                                                                                      <  2:506510/142‑1 (MQ=255)
                                                                                                                                                                       cGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAGTGGTGAGGATGCCGATGCCCAAgg                                        <  2:681131/235‑1 (MQ=255)
                                                                                                                                                                         tGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCTGCTCGGCCACCTGGTCGAGCAGGCCGAAGCCGGTGACGTCGGt                                                                                                                                                                                    >  2:979401/1‑93 (MQ=25)
                                                                                                                                                                         tGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCTGCTCGGCCACCTGGTCGAGCAGGCCGAAGCCGGTGACGTCGGt                                                                                                                                                                                    <  1:979406/93‑1 (MQ=37)
                                                                                                                                                                         tGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCa                                                                                                                     >  1:341639/1‑156 (MQ=255)
                                                                                                                                                                         tGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAGTGGTGAGGATGCCGATGCCCAAGGGCt                                     >  1:591532/1‑236 (MQ=255)
                                                                                                                                                                         tGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGGGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGc                                                                              >  2:1050069/1‑195 (MQ=255)
                                                                                                                                                                          ggCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAGTGGTGAGGATGCCGATGCCCAAGGGCtt                                    >  1:474194/1‑236 (MQ=255)
                                                                                                                                                                               ggcggcGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTc                                                                                                <  2:50928/171‑1 (MQ=255)
                                                                                                                                                                                   gcgTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCCCACGCAGCTTGGCCTTCttt                                                                       >  2:1000921/1‑192 (MQ=255)
                                                                                                                                                                                   gcgTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAGTGGTGAGGATGCCGATGCCCAAGGGc                                      >  1:407664/1‑225 (MQ=255)
                                                                                                                                                                                         ttCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAGTGGTGAGGATGCCGATGCCCAAGGGCTTGGTCAGGTACAGCt                      <  1:10261/235‑1 (MQ=255)
                                                                                                                                                                                                 cgcgGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGGCCTGGCCTTGACCCCGTCCCGTTTGCCGAAGCGGCTGCCGGGCggg                                                                                                                                      >  2:179289/1‑113 (MQ=255)
                                                                                                                                                                                                     gTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAGTGGTGa                                                          >  2:1173606/1‑187 (MQ=255)
                                                                                                                                                                                                           ccGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTc                                                                     >  1:989088/1‑170 (MQ=255)
                                                                                                                                                                                                            cGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAGTGGTGAGGATGCCGATGCCCAAgg                                        <  1:437385/198‑1 (MQ=255)
                                                                                                                                                                                                               tGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAGTGGTGAGGATGCCGATGCCCAAGGGCTTGGTCAGGTACAGCTGGCAGCCTGCGCTGGCGgtg  <  2:591527/233‑1 (MQ=255)
                                                                                                                                                                                                               cggcTTCGGCCTGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACCCGCCGCTTGGCCTTCTTTTCGGCAGTGGTGAg                                                         >  1:278203/4‑178 (MQ=255)
                                                                                                                                                                                                                 ccTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAGTGGTGCGGATGCCGATGCCCAAGGGCTTGGTCAGGTACAGCTgg                    >  1:304172/1‑213 (MQ=255)
                                                                                                                                                                                                                 ccTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAGTGGTGAGGATGCCGATGCCCAAgg                                        >  2:173068/1‑193 (MQ=255)
                                                                                                                                                                                                                 ccTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAGTGGTGAGGATGCCGATGCCCAAgg                                        <  1:173071/193‑1 (MQ=255)
                                                                                                                                                                                                                      ggCCTGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAGTGGTGAGGATGCCGATGCCCAAGGGCTTGGTCAGGTACAGCTGGCAGCCTGCGCTg        >  1:1025203/1‑220 (MQ=255)
                                                                                                                                                                                                                       gCCTGCTCGTCCACCTGGTTGGGCAGGCCGAGGCCGGTGACGGCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGcc                                                                                                                                             <  1:1134332/86‑1 (MQ=255)
                                                                                                                                                                                                                        ccAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGCAGCGGCTGCCGCGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACCCGCAg                                                                                    >  2:52049/1‑142 (MQ=255)
                                                                                                                                                                                                                          |                                                                                                                                                                                                                             
CCTCCGGGGCAACAGCTACCAGCAAGCCACCGCTGGTCTGCGGGTCGCACAGCAGGTGCTTCTGGTCGTCGCTCAGGGCGCTGATCTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAGTGGTGAGGATGCCGATGCCCAAGGGCTTGGTCAGGTACAGCTGGCAGCCTGCGCTGGCGGTG  >  AE015451/961481‑961920

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