Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 205,450 C→T 8.3% G3652G (GGC→GGT PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451205,4500CT 8.3% 198.7 / 5.0 71G3652G (GGC→GGTPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (26/39);  new base T (2/4);  total (28/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.77e-01

GACACTCGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGTACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCAC  >  AE015451/205222‑205671
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                             
gACACTCGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGTACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCa                                                                                                                                                                                                                        <  1:530984/236‑1 (MQ=40)
gACACTCGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGTACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCa                                                                                                                                                                                                                        <  1:341693/236‑1 (MQ=39)
 acacTCGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGTACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAAc                                                                                                                                                                                                                          >  1:228707/1‑233 (MQ=38)
     tCGCCAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGTACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAgcgcg                                                                                                                                                                                                                   <  1:766197/236‑1 (MQ=37)
              cgccGGTCAGCACTACCGTCACCGACGAGCCTGGTACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAgc                                                                                                                                                                                                                      <  2:639848/222‑1 (MQ=11)
                  ggTCAGCACTACCGTCACCGACGAGCCTGGTACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAgcgcg                                                                                                                                                                                                                   <  2:602864/223‑1 (MQ=34)
                  ggTCAGCACTACCGTCACCGACGAGCCTGGTACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCCTACACCCa                                                                                                                                                                                                        >  2:496165/1‑234 (MQ=18)
                      aGCACTACCGTCACCGACGAGCCTGGTACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTg                                                                                                                                                                                                   >  1:831700/1‑235 (MQ=34)
                           tACCGTCACCGACGAGCCTGGTACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGACACCAGCGCGcctg                                                                                                                                                                                                               >  1:304911/1‑216 (MQ=21)
                           tACCGTCACCGACGAGCCTGGTACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTgcg                                                                                                                                                                                                 >  1:292142/1‑232 (MQ=34)
                                tCACCGACGAGCCTGGTACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGc                                                                                                                                                                                           <  2:831697/233‑1 (MQ=38)
                                     gacgaGCCTGGTACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGcc                                                                                                                                                                                                                 >  1:229441/1‑206 (MQ=34)
                                               gTACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACg                                                                                                                                                                           <  2:498536/234‑1 (MQ=39)
                                                tACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACg                                                                                                                                                                           <  2:323646/233‑1 (MQ=39)
                                                       ggCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCa                                                                                                                                                                 <  1:811767/236‑1 (MQ=255)
                                                        gCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTgc                                                                                                                                                                                                  >  2:396123/1‑202 (MQ=35)
                                                           aCCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGACCAGATc                                                                                                                                                             <  1:679037/236‑1 (MQ=255)
                                                             ccGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGc                                                                                                                                                                                           <  2:936678/204‑1 (MQ=18)
                                                                 gcggcAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCa                                                                                                                                                            <  2:365619/231‑1 (MQ=38)
                                                                         aaCGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGa                                                                                                                                                                                                                             >  2:312613/1‑158 (MQ=14)
                                                                           cGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAgcgcg                                                                                                                                                                                                                   <  1:265581/166‑1 (MQ=35)
                                                                             aaGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCacgtc                                                                                                                                            <  1:757764/235‑5 (MQ=9)
                                                                              aGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGcc                                                                                                                                                                  <  2:518359/212‑1 (MQ=35)
                                                                                tcGACCTGGTCACGGTCACCATCACCGCCGACCAGCCTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCa                                                                                                                                                                                                                        <  2:32247/155‑1 (MQ=255)
                                                                                        gTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAgcgcg                                                                                                                                                                                                                   >  1:317245/1‑153 (MQ=38)
                                                                                        gTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAgcgcg                                                                                                                                                                                                                   <  2:317242/153‑1 (MQ=38)
                                                                                          cAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGAt                                                                                                                                                                                        <  2:803405/178‑1 (MQ=12)
                                                                                          cAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGAt                                                                                                                                                                                        >  1:803408/1‑178 (MQ=37)
                                                                                                 aCCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACg                                                                                                                                                                           <  1:171273/184‑1 (MQ=37)
                                                                                                      cACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGAt                                                                                                                                                                                        >  2:757418/1‑166 (MQ=16)
                                                                                                      cACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAgcgcg                                                                                                                                                                                                                   <  1:191696/139‑1 (MQ=37)
                                                                                                       aCCGCCGCCCAGACTCCGGTAGCCGAGAACGTCAAACCGCCCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCt                                                                                                                                                                      <  2:101366/183‑1 (MQ=255)
                                                                                                       aCCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCTGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCtt                                                                                                                 <  1:511548/236‑1 (MQ=38)
                                                                                                        ccgccgACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACt                                                                                                                                             <  2:928162/207‑1 (MQ=35)
                                                                                                         ggccgACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAg                                                                                                                                                                <  2:478872/186‑1 (MQ=21)
                                                                                                          gccgACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGc                                                                                                                                                                                                                               >  2:39618/1‑123 (MQ=33)
                                                                                                            cgACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTa                                                                                                                                                                                <  2:900688/168‑1 (MQ=37)
                                                                                                            cgACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTa                                                                                                                                                                                >  1:900691/1‑168 (MQ=37)
                                                                                                            cgACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACCGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAg                                                                                                                                                                                                                       <  2:519855/129‑1 (MQ=17)
                                                                                                              aCCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCAGGCTgc                                                                                                                                                                                                  >  1:161927/1‑148 (MQ=17)
                                                                                                                cAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGcctg                                                                                                                                                                                        >  1:202566/1‑153 (MQ=14)
                                                                                                                cAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACg                                                                                                                                                                                    >  2:235227/1‑160 (MQ=35)
                                                                                                                            gCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGggcgg                                                                                            <  1:71364/236‑1 (MQ=38)
                                                                                                                                        aagccGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAg                                                                                                                                                                                                                       <  2:282042/98‑1 (MQ=17)
                                                                                                                                        aaaCCGACCTTACCTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATg                                                                                                                                 <  2:635073/187‑1 (MQ=37)
                                                                                                                                         aaCCGCCCTTCCCTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCGGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCa                                                                                                                                                            <  2:98320/159‑1 (MQ=255)
                                                                                                                                               ccTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGt                                                                          <  1:481093/235‑1 (MQ=255)
                                                                                                                                                 ttCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGAc                                                                       >  1:786711/1‑236 (MQ=40)
                                                                                                                                                       gTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCCCCGTACACCCACGCTGCGCAGGGc                                                                                                                                                                                           <  1:517440/114‑1 (MQ=17)
                                                                                                                                                                aaCCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGAc                                                                                                                    >  2:532886/1‑176 (MQ=33)
                                                                                                                                                                     gCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACcg                                                            <  1:312616/227‑1 (MQ=39)
                                                                                                                                                                              tCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAg                                                      >  2:437469/1‑224 (MQ=40)
                                                                                                                                                                               cACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTc                                                                              >  2:485973/1‑199 (MQ=40)
                                                                                                                                                                                          cGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCCACGCTTCGGTTCaa                                                                            >  1:217426/1‑190 (MQ=255)
                                                                                                                                                                                            tGACCCTGAGCAACACCGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACcg                                                            <  1:634815/204‑1 (MQ=255)
                                                                                                                                                                                              aCCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGt                           <  1:583978/235‑1 (MQ=255)
                                                                                                                                                                                                 cTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTACAACTGGGTGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCg                         >  1:523686/1‑234 (MQ=255)
                                                                                                                                                                                                      caacaaTGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTc                                                                                   >  2:440974/1‑171 (MQ=255)
                                                                                                                                                                                                      caacaaCGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGGTCAAGTGACCGACACCACCGACGAAGGGGGGGc                                               >  1:165177/1‑207 (MQ=255)
                                                                                                                                                                                                           gcGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTTCAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGGCGGTCACCGAAgg              <  1:575854/234‑1 (MQ=255)
                                                                                                                                                                                                            cGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGaa                >  1:411536/1‑232 (MQ=40)
                                                                                                                                                                                                               ccAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCTCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCACCTCGGCTGTGGATGCCACTGGTGCCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAggcggc          <  1:420315/235‑1 (MQ=255)
                                                                                                                                                                                                               ccAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAggcggc          <  1:916866/235‑1 (MQ=40)
                                                                                                                                                                                                                 aGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACCCCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAggcggc          <  1:255974/233‑1 (MQ=38)
                                                                                                                                                                                                                    tCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCg                                                                     >  2:858236/1‑171 (MQ=40)
                                                                                                                                                                                                                       ccATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCAc  <  1:532889/235‑1 (MQ=39)
                                                                                                                                                                                                                           cAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAggcggc          <  1:846136/223‑1 (MQ=39)
                                                                                                                                                                                                                             aGGCCGGTAAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGAcc                                     <  2:556868/194‑1 (MQ=255)
                                                                                                                                                                                                                                ccGGTGAAACCAGCGCGCCGTACACCCACGCTGCTCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGGAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCAc  >  2:478112/1‑226 (MQ=255)
                                                                                                                                                                                                                                ccGGTGAAACCAGCGCGCCGTACACCCACGCTGCTCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCAc  <  1:478115/226‑1 (MQ=255)
                                                                                                                                                                                                                                ccGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGCCGAAg                                                      >  2:141275/1‑174 (MQ=255)
                                                                                                                                                                                                                                ccGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCgaga      >  1:53148/1‑222 (MQ=39)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                             
GACACTCGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGTACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCAC  >  AE015451/205222‑205671

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