Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 216,949 T→C 6.3% D7485D (GAT→GAC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451216,9490TC 6.3% 239.5 / 4.8 77D7485D (GAT→GACPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (40/32);  new base C (3/2);  total (43/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.00e-01

CGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGACACTGGACAAGACCCCGGTCAGCACCGCCGTGA  >  AE015451/216729‑217181
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                        
cgcAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATcaa                                                                                                                                                                                                                            <  1:749273/235‑1 (MQ=34)
     ggCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAAcc                                                                                                                                                                                                                       <  1:432250/235‑1 (MQ=34)
     ggCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAAcc                                                                                                                                                                                                                       <  1:526261/235‑1 (MQ=37)
           gACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGc                                                                                                                                                                                                                <  1:1268281/236‑1 (MQ=40)
             cGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGc                                                                                                                                                                                                                >  2:1080884/1‑234 (MQ=38)
                  aCAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCCACCCCTTCGGTGACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAAcc                                                                                                                                                                                                                       >  1:850521/1‑222 (MQ=2)
                        aTGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATc                                                                                                                                                                                                                              <  2:361506/209‑1 (MQ=32)
                           cTGGCCAGATCAGCCGCGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCGTCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCAt                                                                                                                                                                                                              <  2:890815/222‑1 (MQ=255)
                            tGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTg                                                                                                                                                                                                >  1:825671/1‑235 (MQ=255)
                            tGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTg                                                                                                                                                                                                >  1:1047111/1‑235 (MQ=255)
                              gCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCa                                                                                                                                                                                                                                 >  1:700492/1‑200 (MQ=32)
                                      aGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCGAGATCACTTACACCATCACCCTGACCAACAAAGACGGCCTGc                                                                                                                                                                                                                                   <  1:395392/190‑1 (MQ=9)
                                                tCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGATGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACGTACACCATCACCCTGGCCAACAAAGATGGTCTGCCAATCAACAACAACAGt                                                                                                                                                                                                                  >  2:686099/1‑197 (MQ=255)
                                                  agcTCGGCGGTGGACGACCCTGGTGCGACCTTCGAGAACCTCGAGCGGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAg                                                                                                                                                                                              <  1:384456/213‑1 (MQ=255)
                                                    cTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTAGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATc                                                                                                                                                                                                                              <  2:694145/181‑1 (MQ=21)
                                                    cTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATc                                                                                                                                                                                                                              >  1:694150/1‑181 (MQ=34)
                                                    cTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCg                                                                                                                                                                       <  1:1017746/236‑1 (MQ=255)
                                                       ggcggTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATc                                                                                                                                                                                                                              <  2:182591/178‑1 (MQ=33)
                                                         cggTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGc                                                                                                                                                                                                                                   >  1:376508/1‑171 (MQ=31)
                                                         cggTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTg                                                                                                                                                                     <  2:655915/233‑1 (MQ=255)
                                                          ggTGGACGCCACTGGTGCGACCTTCGAGAACCTCGTGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGGCCCCTTCGCTAACGGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAAcc                                                                                                                                                                                                                       <  2:693673/182‑1 (MQ=255)
                                                          ggTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGACGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTCACCTTCACGCtc                                                                                                                                                                                                >  2:65136/1‑204 (MQ=255)
                                                                      tGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTg                                                                                                                                                                                                            >  2:734533/1‑181 (MQ=255)
                                                                                ttCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGCGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGcct                                                                                                                                                                                                              >  1:390284/1‑167 (MQ=255)
                                                                                ttCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATg                                                                                                                                                                                                                                         >  1:1010626/1‑142 (MQ=32)
                                                                                ttCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATg                                                                                                                                                                                                                                         <  2:1010621/142‑1 (MQ=32)
                                                                                ttCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACgg                                                                                                                                                          >  2:1216632/1‑221 (MQ=255)
                                                                                ttCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACt                                                                                                                                             >  2:1199769/1‑234 (MQ=255)
                                                                                gccGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGa                                                                                                                                                                                                           <  2:25391/170‑1 (MQ=255)
                                                                                    agaACCGCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGc                                                                                                                                                                                                                                   <  2:72823/144‑1 (MQ=17)
                                                                                      aaCCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCg                                                                                                                                      <  1:1080889/235‑1 (MQ=255)
                                                                                           cGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGt                                                                                                                                                                                                                  >  1:37862/1‑154 (MQ=255)
                                                                                           cGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGt                                                                                                                                                                                                                  <  2:37861/154‑1 (MQ=255)
                                                                                           cGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTg                                                                                                                                >  1:140475/1‑236 (MQ=255)
                                                                                           cGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTg                                                                                                                                <  1:486267/236‑1 (MQ=255)
                                                                                                      gTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCCCCg                                                                                                                                                     >  1:569615/1‑204 (MQ=255)
                                                                                                           gcgAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTAcc                                                                                                                                       >  1:56804/1‑213 (MQ=255)
                                                                                                            cgAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTc                                                                                                                  <  2:492176/233‑1 (MQ=255)
                                                                                                               aGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAACATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTc                                                                                                            <  1:744847/236‑1 (MQ=255)
                                                                                                                ggTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACg                                                                                                              <  2:448252/233‑1 (MQ=255)
                                                                                                                    gATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTAc                                                                                                       <  1:1199774/236‑1 (MQ=255)
                                                                                                                        tCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCg                                                                                                                           <  1:275588/212‑1 (MQ=39)
                                                                                                                            cGACACCACCGACGATGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGACGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGAc                                                                                                <  1:712161/235‑1 (MQ=255)
                                                                                                                              acacCACCGACGTAGTGGTTGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGCTCACCTACACCATCACCCTGCTCAACAAAGATGGTCTGCCCATCAACAACCACAGTGCATTGAcc                                                                                                                                                                                                         >  2:106123/1‑128 (MQ=14)
                                                                                                                              acacCACCGACGTAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGATCAACAAAGATGGTCTGCCAATCAACAACACCAGTGCATTGACCTTCACGCTGAGCGAccg                                                                                                                                                                                        >  2:106122/1‑143 (MQ=255)
                                                                                                                                     cgacgaAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGc                                                                                                                                           >  1:572881/1‑183 (MQ=39)
                                                                                                                                     cgacgaAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGc                                                                                                                                           <  2:572876/183‑1 (MQ=40)
                                                                                                                                        cgaggtggtgGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCAATCAACAACCACACGGCATTGAACTTCACGCTGACCGACGGCAAAACCGTCATCACCGTGCCga                                                                                                                                                                 >  2:505001/5‑157 (MQ=11)
                                                                                                                                        cgaAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCCACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCAc                                                                                                                                                       >  2:41800/1‑168 (MQ=255)
                                                                                                                                        cgaAGTGGTGGCCAAACTGCCCGCCACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCg                                                                                                                                                                       <  2:285908/152‑1 (MQ=255)
                                                                                                                                               gtggCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGa                                                                                                                         >  2:175813/1‑191 (MQ=38)
                                                                                                                                                tggCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAatcgatc                                                                              <  2:691532/233‑1 (MQ=39)
                                                                                                                                                 ggCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGAcc                                                                                               >  2:855304/1‑215 (MQ=255)
                                                                                                                                                   ccAAGCTGACCGCAACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGATGGTCTGCCGATCAACAACcacgc                                                                                                                                                                                                                   >  1:791135/1‑95 (MQ=11)
                                                                                                                                                         tGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTggc                                                                  >  1:188351/1‑236 (MQ=39)
                                                                                                                                                          gACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCt                                                                                                                                                                                                                                     <  1:929375/72‑1 (MQ=35)
                                                                                                                                                                    ccTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGATGGTCTGCAGATCAACAACcac                                                                                                                                                                                                                     >  2:693930/1‑78 (MQ=11)
                                                                                                                                                                    ccTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGt                                                                                                                   >  2:316769/1‑176 (MQ=38)
                                                                                                                                                                        cGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCATGCCGGCCAACgg                                                                                                                                                          >  1:34430/1‑133 (MQ=37)
                                                                                                                                                                           tAACCGAAGGCGGCGAGATCACCAACACCATCACCCTGACCAACAAAGGTGGTCTGCCAATCAACAACCACAGTACAGTGACCTGCACGCTGAGCGACGGCTAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCg                                                                                                                                      <  1:879667/150‑1 (MQ=11)
                                                                                                                                                                                gAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCATGTGATCAAat                                                                                   >  2:1247732/1‑196 (MQ=38)
                                                                                                                                                                                     cggcgAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCTGTACCa                                                                                                     >  1:18285/1‑173 (MQ=25)
                                                                                                                                                                                      ggcgAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGt                                                                                                                                     >  2:233959/1‑140 (MQ=34)
                                                                                                                                                                                                 cTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGa                                                                                                                                                                                                           <  2:432201/59‑1 (MQ=35)
                                                                                                                                                                                                        aTCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGACACTGGACAAGAcccc                   <  1:762026/236‑1 (MQ=255)
                                                                                                                                                                                                            cccTGACCAACAAAGTTGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCACCGGCCCCGTGGGCACTGCTACCGTCCCTGCCCCGGACAACGTCGACGt                                                                                                             >  2:442986/1‑142 (MQ=11)
                                                                                                                                                                                                            cccTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCCTTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGAACCTGTGATCAAATCGATCGCAACCgct                                                                     >  2:308323/1‑180 (MQ=21)
                                                                                                                                                                                                            cccTGACCAACAAAGATGGCCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCAcc                                                                                                                                                                        <  1:445344/83‑1 (MQ=18)
                                                                                                                                                                                                              cTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGACACTGGACAAGACCCCGGTCa              >  1:898587/1‑235 (MQ=255)
                                                                                                                                                                                                                  ccAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGCCCCTGTGATCAAatcgat                                                                               >  2:178691/1‑166 (MQ=18)
                                                                                                                                                                                                                    aacaaAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGGCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGAGGTTGGCAAGTTCGAAc                                         >  1:119050/1‑202 (MQ=255)
                                                                                                                                                                                                                    aacaaAGATGGTCTGCCAATCAACAACCACAGTGCATTCACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCCACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGt                                                                                                                   >  2:53895/1‑128 (MQ=17)
                                                                                                                                                                                                                    aacaaAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGTGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGACACTGGACAAGACCCCGGTCa              >  2:566241/1‑229 (MQ=255)
                                                                                                                                                                                                                    aacaaAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGACACTGGACAAGACCCCGGTCa              <  1:566246/229‑1 (MQ=255)
                                                                                                                                                                                                                        aaGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTGGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTAACACTGGACAAGACCCCGGTCAGCAccgcc      >  1:1205472/1‑233 (MQ=255)
                                                                                                                                                                                                                         aGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGa                                                          >  2:1178343/1‑180 (MQ=40)
                                                                                                                                                                                                                         aGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGACACTGGACAAGACCCCGGTCAGCACCGCCGTGa  >  1:278111/1‑236 (MQ=255)
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                        
CGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGACACTGGACAAGACCCCGGTCAGCACCGCCGTGA  >  AE015451/216729‑217181

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: