Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 198,802 T→C 5.8% A1436A (GCT→GCC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451198,8020TC 5.8% 284.6 / 4.4 84A1436A (GCT→GCCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (41/38);  new base C (3/2);  total (44/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

CGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGC  >  AE015451/198583‑199033
                                                                                                                                                                                                                           |                                                                                                                                                                                                                                       
cgcAGATGCAGGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTc                                                                                                                                                                                                                                        <  1:1248415/221‑1 (MQ=255)
             aaGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAAcc                                                                                                                                                                                                              <  2:772538/234‑1 (MQ=255)
             aaGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAAcc                                                                                                                                                                                                              <  2:1077862/234‑1 (MQ=255)
                  cGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAgcgc                                                                                                                                                                                                         >  2:664206/1‑234 (MQ=255)
                  cGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCAc                                                                                                                                                                                                       >  1:919041/1‑236 (MQ=255)
                            aacctgGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTc                                                                                                                                                                                              <  1:369648/235‑1 (MQ=255)
                                       gACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCgg                                                                                                                                                                                                                     <  2:1227206/201‑1 (MQ=255)
                                         ccccGGTCAGCACCGCGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAAcc                                                                                                                                                                                                              <  1:303758/206‑1 (MQ=255)
                                          cccGGTCAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACCCCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTa                                                                                                                                                                                                   <  2:1020631/216‑1 (MQ=255)
                                          cccGGTCAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAgg                                                                                                                                                                                  <  1:1110452/233‑1 (MQ=255)
                                          cccGGTCAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGc                                                                                                                                                                                <  1:248683/235‑1 (MQ=255)
                                               tCAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGAc                                                                                                                                                                          >  1:397972/1‑236 (MQ=255)
                                                      ccaGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGTCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCg                                                                                                                                                                               <  2:215333/222‑1 (MQ=255)
                                                      cGAGGTCACCGACGAGCCAGGGACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAgc                                                                                                                                                                                                           <  2:22135/196‑1 (MQ=255)
                                                          gTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAgg                                                                                                                                                                <  1:664211/235‑1 (MQ=255)
                                                                   gaGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCa                                                                                                                                                                                                                                       <  2:26210/155‑1 (MQ=255)
                                                                    agacAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACt                                                                                                                                                                                               <  2:366395/191‑1 (MQ=255)
                                                                       cAGGTCCGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAACCCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTa                                                                                                                                                                     <  2:320745/217‑1 (MQ=255)
                                                                         ggTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCa                                                                                                                                                 <  1:1181272/235‑1 (MQ=255)
                                                                          ttACGCCAGGCAACGAAGGCGATATCGCCCAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGAt                                                                                                                                                                             <  1:74250/205‑1 (MQ=255)
                                                                           tACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAGACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGc                                                                                                                                               >  1:424387/1‑235 (MQ=255)
                                                                              gCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCgg                                                                                                                                                                                                                     <  2:270529/162‑1 (MQ=255)
                                                                                     aaCGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGTTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACTCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGag                                                                                                                                                   >  1:312229/1‑220 (MQ=255)
                                                                                      aCGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGGACACTCACGACGCGCAggg                                                                                                                                                                                 >  2:271972/1‑190 (MQ=255)
                                                                                      aCGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCg                                                                                                                                                                                                     >  2:907031/1‑170 (MQ=255)
                                                                                       cGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAgg                                                                                                                                                                                                                         >  2:139685/1‑149 (MQ=255)
                                                                                          aGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTcgg                                                                                                                               >  1:939761/1‑236 (MQ=255)
                                                                                             cGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGt                                                                                                                                                                                                                                    >  1:341851/1‑132 (MQ=255)
                                                                                               atatCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCAGATCAGCCTGGGCATCAACTCGGCGGTg                                                                                                                          <  1:666548/235‑1 (MQ=255)
                                                                                                tatCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTg                                                                                                                          <  1:794507/235‑1 (MQ=255)
                                                                                                  tCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCacgac                                                                                                                                                                                         >  1:551825/1‑170 (MQ=255)
                                                                                                           tCACCCTCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGCCCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAggg                                                                                                                                                                                 <  2:224342/169‑1 (MQ=255)
                                                                                                               cctCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGaaa                                                                                                                                                                                                                <  2:321331/132‑1 (MQ=255)
                                                                                                                 tCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCa                                                                                                                                                                                                        <  2:529999/140‑1 (MQ=255)
                                                                                                                   aCGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGt                                                                                                       <  1:420962/235‑1 (MQ=255)
                                                                                                                             aGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGAt                                                                                                                                                                             >  2:36670/1‑155 (MQ=255)
                                                                                                                             aGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTc                                                                                            >  1:799220/1‑236 (MQ=255)
                                                                                                                                  tCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGtt                                                                                       <  1:1217254/236‑1 (MQ=255)
                                                                                                                                    ggTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCAcag                                                                                                                >  1:220364/1‑207 (MQ=255)
                                                                                                                                         ccGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTc                                                                                            <  2:919036/224‑1 (MQ=255)
                                                                                                                                          cGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGtg                                                                                                             <  1:120419/206‑1 (MQ=255)
                                                                                                                                            agaACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACg                                                                                                                      >  2:1251000/1‑195 (MQ=255)
                                                                                                                                            agaACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCt                                                                                             >  2:1246754/1‑220 (MQ=255)
                                                                                                                                                             ttCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCga                                                                                                   >  1:387896/1‑197 (MQ=255)
                                                                                                                                                                 cTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTg                                                                    <  2:723798/224‑1 (MQ=255)
                                                                                                                                                                         aTCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAAcc                                                                                              >  1:608940/1‑190 (MQ=255)
                                                                                                                                                                            aaCACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCAc                                                                                                                                                                                                                                 >  1:880149/1‑56 (MQ=255)
                                                                                                                                                                            aaCACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCt                                                   >  2:565745/1‑230 (MQ=255)
                                                                                                                                                                            aaCACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGCTGGGCGGTGctg                                                                            >  1:667149/1‑203 (MQ=255)
                                                                                                                                                                            aaCACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTAACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTcggcgg                                                                                                                            >  1:1141334/1‑157 (MQ=255)
                                                                                                                                                                            aaCACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCAc                                                                                                                                                                                                                                 <  2:880144/56‑1 (MQ=255)
                                                                                                                                                                             acacCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGACACCGACACCCTcgacga                                             <  1:1093595/235‑1 (MQ=255)
                                                                                                                                                                               acCGCCCTGGCTCACGACCTGGTCGTGGCCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATCACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACAccc                                                    >  1:678401/1‑226 (MQ=255)
                                                                                                                                                                               acCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCcgcg                                                                          <  2:928748/204‑1 (MQ=255)
                                                                                                                                                                               acCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCcgcg                                                                          >  1:928753/1‑204 (MQ=255)
                                                                                                                                                                               acCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAgg                                           <  1:422518/235‑1 (MQ=255)
                                                                                                                                                                                 cGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCAGCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCt                                                                                                                                                                      >  2:1266548/1‑110 (MQ=255)
                                                                                                                                                                                 cGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAggg                                                                                                                                                                                 <  2:656662/99‑1 (MQ=255)
                                                                                                                                                                                 cGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAggg                                                                                                                                                                                 >  1:656667/1‑99 (MQ=255)
                                                                                                                                                                                 cGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCgaga                                                                                                 <  1:168245/179‑1 (MQ=255)
                                                                                                                                                                                 cGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTcgacg                                              >  1:1206949/1‑230 (MQ=255)
                                                                                                                                                                                 cGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGCGGACGCCa                                                                                                                   >  1:306484/1‑161 (MQ=255)
                                                                                                                                                                                 cGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCAACTCGGCGGCGGACGCCACt                                                                                                                 >  1:390820/1‑163 (MQ=255)
                                                                                                                                                                                 cGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGACGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCATTCGAGAACCTTGCGCTGGGCGGTGCTGCGTCGGTCCACGTCACCGACACCCTcgac                                               >  2:1012040/1‑229 (MQ=11)
                                                                                                                                                                                    ccTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGtgct                                                                                                           <  2:608935/166‑1 (MQ=255)
                                                                                                                                                                                     cTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAGgtggtg                                      <  1:1193901/234‑1 (MQ=255)
                                                                                                                                                                                      tGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAgg                                                                                                                                                                <  1:1092289/111‑1 (MQ=255)
                                                                                                                                                                                      tGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAgg                                                                                                                                                                >  2:1092284/1‑111 (MQ=255)
                                                                                                                                                                                      tGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTcggcggc                                                                                                                           >  2:357921/1‑147 (MQ=255)
                                                                                                                                                                                      tGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAgg                                           >  1:686765/1‑228 (MQ=255)
                                                                                                                                                                                      tGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAgg                                           >  1:780574/1‑228 (MQ=255)
                                                                                                                                                                                      tGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGc                                    >  1:807599/1‑235 (MQ=255)
                                                                                                                                                                                              gACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCt                              >  2:1248430/1‑233 (MQ=255)
                                                                                                                                                                                                       gTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTcg                                                 >  1:1237203/1‑205 (MQ=11)
                                                                                                                                                                                                       gTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTcg                                                 <  2:1237198/205‑1 (MQ=11)
                                                                                                                                                                                                           cccTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCAc                                                            <  1:923949/190‑1 (MQ=255)
                                                                                                                                                                                                             cTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGa                                                         >  1:586685/1‑191 (MQ=255)
                                                                                                                                                                                                             cTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGGGCTGCCTTCGAGAACCTCg                                                                                           >  2:68333/1‑157 (MQ=255)
                                                                                                                                                                                                                 gCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCgaga                                                                                                 >  2:655989/1‑147 (MQ=255)
                                                                                                                                                                                                                 gCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCgaga                                                                                                 <  1:655994/147‑1 (MQ=255)
                                                                                                                                                                                                                  caacaaCGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAggg                                          >  1:283492/1‑200 (MQ=255)
                                                                                                                                                                                                                  caacaaCGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCAcc        >  1:1195775/1‑235 (MQ=255)
                                                                                                                                                                                                                   aacaacGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTc               >  2:1255051/1‑227 (MQ=255)
                                                                                                                                                                                                                        cGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAggc  >  1:725046/1‑235 (MQ=255)
                                                                                                                                                                                                                           |                                                                                                                                                                                                                                       
CGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAAGTGGATGTCACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGC  >  AE015451/198583‑199033

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