Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 195,732 A→C 16.4% N413T (AAC→ACC)  PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451195,7320AC16.4% 212.7 / 32.8 93N413T (AAC→ACC) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base A (47/29);  new base C (7/9);  total (54/39)
Fisher's exact test for biased strand distribution p-value = 2.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.03e-01

CGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCAC  >  AE015451/195523‑195936
                                                                                                                                                                                                                 |                                                                                                                                                                                                            
cGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACg                                                                                                                                                                                                   >  2:672143/1‑221 (MQ=12)
   ttCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTg                                                                                                                                                                                 >  1:39820/1‑236 (MQ=16)
    tCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGa                                                                                                                                                                                <  1:996755/236‑1 (MQ=255)
                   aCCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTg                                                                                                                                                                                                          >  1:737002/1‑195 (MQ=255)
                                    aggcggTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTcc                                                                                                                                                >  1:795760/2‑236 (MQ=255)
                                            aTTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGccg                                                                                                                                                                           <  1:991877/201‑3 (MQ=2)
                                              tgcGAGCAACTGACCCCGACTCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACCATTAcc                                                                                                                                                                                     <  1:661016/187‑1 (MQ=255)
                                                            ccggAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGccg                                                                                                                                                                           <  1:124073/182‑3 (MQ=2)
                                                               gaacccGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGa                                                                                                                                                                                >  1:631477/5‑177 (MQ=255)
                                                               gAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGa                                                                                                                                                                                <  2:631475/177‑1 (MQ=255)
                                                                    cGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                           <  1:543885/177‑1 (MQ=40)
                                                                          cGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAAt                                                                                                                                  <  2:599781/212‑1 (MQ=25)
                                                                           ggCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGtt                                                                                                                                                                                  <  2:1238793/163‑1 (MQ=12)
                                                                               cAGCCCACGACCACCGACTCGGTCGACGCCACCACCGCGACCCTCACGGCGAGCCCGTCGGTCGCCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACgg                                                                                                                                                                                                  <  2:457636/143‑1 (MQ=14)
                                                                                          cAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAAc                                                                                                                                                                                                    >  2:1236842/1‑130 (MQ=255)
                                                                                          cAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAAc                                                                                                                                                                                                    <  1:1236845/130‑1 (MQ=255)
                                                                                                    gTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGcc                                                                                  >  1:491548/1‑234 (MQ=255)
                                                                                                      cgatgccaccaccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                           >  1:983427/6‑143 (MQ=255)
                                                                                                      cgatgccaccaccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                           <  2:983425/138‑1 (MQ=255)
                                                                                                      cGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGCAACCGTGAACCTGTCCAACGGc                                                                                                                                                                                                 >  1:273702/1‑121 (MQ=255)
                                                                                                      cGACACCACTACCGCGACCCTGACGACGAGCCCGTCGGTCACCGAAGGAGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCa                                                                                 >  1:839929/1‑233 (MQ=255)
                                                                                                      cGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGcc                                                                                                                                                                                                >  1:519775/1‑122 (MQ=25)
                                                                                                                  cgcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACgg                                                                                                                 >  1:1253799/1‑189 (MQ=37)
                                                                                                                  cgcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTc                                                                   >  1:1409269/1‑235 (MQ=255)
                                                                                                                         cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAACCGACTAcctc                                                                                                                                                                                   >  2:493501/1‑114 (MQ=255)
                                                                                                                         cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGa                                                                                                                                >  1:55301/1‑167 (MQ=38)
                                                                                                                         cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGc                                                            >  1:787387/1‑235 (MQ=255)
                                                                                                                               gCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGAcc                                                                                                                                           >  1:389806/1‑150 (MQ=39)
                                                                                                                                    cccGTCGGTCACCGAAGGCGGCGTGATCACTTACACCTGTACCCTAAACACATCTGCCCACCCACCAGTAACCGTGACACTGTCCAACCCCCAAACGATTACCGTAGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACggtt                                                                                                               >  1:705607/1‑171 (MQ=11)
                                                                                                                                      cGTCGGTGACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                           <  2:797562/111‑1 (MQ=255)
                                                                                                                                         cGGTCACCGAAGGGGGCGTGATCACTTACCCCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCa                                                                                                                                    <  2:313529/147‑1 (MQ=255)
                                                                                                                                           gtgaCCGAAGGTGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGt                                                                                                                             <  2:175697/149‑1 (MQ=16)
                                                                                                                                           gtgaCCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACa                                                              <  1:64729/212‑1 (MQ=9)
                                                                                                                                            tCACCGAAGGTGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGt                                                                                                                                                        <  2:772713/124‑1 (MQ=255)
                                                                                                                                            tCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAGCGATCACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCa                                                                                                      <  2:378212/174‑1 (MQ=12)
                                                                                                                                                         cgTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAAc                                                                                                                                                                                            <  2:272921/75‑1 (MQ=255)
                                                                                                                                                           tGATCACTTACACCGTGATCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                           <  1:1479188/90‑1 (MQ=255)
                                                                                                                                                           tGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGGCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGAcc                                                                                                          >  1:464821/1‑155 (MQ=255)
                                                                                                                                                           tGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGCCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTc                                                                                                              >  1:223860/1‑151 (MQ=255)
                                                                                                                                                            gATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATc                         >  1:739644/1‑235 (MQ=255)
                                                                                                                                                             aTCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGac                                                                >  1:1139372/1‑194 (MQ=37)
                                                                                                                                                                        cGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCg                                                                                                                                                                            <  2:1041621/76‑1 (MQ=255)
                                                                                                                                                                        cGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGTCCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCg                                              <  1:220796/202‑1 (MQ=255)
                                                                                                                                                                           gACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGa                                                                 >  2:487983/1‑180 (MQ=9)
                                                                                                                                                                           gACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGtt                                                                                                                                                                                  <  1:265213/67‑1 (MQ=255)
                                                                                                                                                                           gACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGtt                                                                                                                                                                                  >  2:265211/1‑67 (MQ=255)
                                                                                                                                                                             cccTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCgagacc                                                                                     >  2:40758/1‑156 (MQ=38)
                                                                                                                                                                               cTGAGCAATCCTGCCCAGACGCCGGTAACAGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGCAAAACCCCGGGCAGCGTCGATTTCAAGACCCCGGCCAATGACGTCTACAACAACGGCTCg                                                                                                             >  2:723497/1‑132 (MQ=11)
                                                                                                                                                                               cTGAGCAACGCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGAcaccac     <  1:1406306/236‑1 (MQ=255)
                                                                                                                                                                               cTGAGCAACCCTGCCCAGACGCCGGTGACCGTGAACCTGTCCAACGGCCAAACCATTAc                                                                                                                                                                                      >  2:1082816/1‑59 (MQ=11)
                                                                                                                                                                               cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACgg                                                                                                                                                                                                  >  1:976053/1‑47 (MQ=255)
                                                                                                                                                                               cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAAc                                                                                                                                                                                            <  2:1039354/53‑1 (MQ=255)
                                                                                                                                                                               cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAg                                                                                                                                              >  1:1007246/1‑99 (MQ=255)
                                                                                                                                                                               cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGaatc                                                >  2:488177/1‑191 (MQ=2)
                                                                                                                                                                               cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCCGCGTCCCCAt                                                                                            >  2:465546/1‑149 (MQ=255)
                                                                                                                                                                               cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCGCGACCGTCAGCGTCACCATCgagaac                                                                                     >  2:420081/1‑154 (MQ=255)
                                                                                                                                                                                  aGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGAcaccaccac  <  1:487985/236‑1 (MQ=255)
                                                                                                                                                                                  aGCAACCCTGCCCAGACGCCGGTCACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAg                                                                                                                                              <  2:1444158/96‑1 (MQ=12)
                                                                                                                                                                                  aGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAGACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGt                                                                                                                                                        <  2:1422964/86‑1 (MQ=255)
                                                                                                                                                                                  aGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAAc                                                                                                                                                                                                    <  1:933246/42‑1 (MQ=255)
                                                                                                                                                                                  aGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAAc                                                                                                                                                                                                    >  2:933244/1‑42 (MQ=255)
                                                                                                                                                                                  aGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTgg                                                                             <  2:510574/161‑1 (MQ=40)
                                                                                                                                                                                   gCAACCCTGCCCAGACGCCGGTAACTGTGAACCTGTCCAACg                                                                                                                                                                                                   >  1:1194313/1‑42 (MQ=255)
                                                                                                                                                                                   gCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                           <  2:334991/66‑1 (MQ=255)
                                                                                                                                                                                   gCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCCTc                                                                                           >  2:391955/1‑146 (MQ=255)
                                                                                                                                                                                   gCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGAcaccaccac  >  1:1404849/1‑235 (MQ=255)
                                                                                                                                                                                     aaCCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAAc                                                                                                                                                                                                    >  2:433815/1‑39 (MQ=255)
                                                                                                                                                                                         cTGCCAAGACGCCGGTAACCGTGATCCTGACCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGACGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGt              <  1:385177/217‑1 (MQ=255)
                                                                                                                                                                                          tGCCCAGACGCCGGTAACCGTGAACCTGTCCAACgg                                                                                                                                                                                                  <  2:338921/36‑1 (MQ=40)
                                                                                                                                                                                          tGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGtt                                                                                                                                                                                  >  2:1362123/1‑52 (MQ=255)
                                                                                                                                                                                          tGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGa                                                                                                                                                                                <  2:1444565/54‑1 (MQ=255)
                                                                                                                                                                                          tGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGa                                                                                                                                                                                >  1:1444568/1‑54 (MQ=255)
                                                                                                                                                                                          tGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGAcc                               <  2:674289/199‑1 (MQ=255)
                                                                                                                                                                                          tGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGAcc                               >  1:674291/1‑199 (MQ=255)
                                                                                                                                                                                          tGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACAGGCGGCAACTTCGaa                                                                >  1:415499/1‑166 (MQ=25)
                                                                                                                                                                                          tGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGGCTACAACAACGGCTCGAACGTCAGCGTCACCATCga                                                                                         >  1:5260/1‑141 (MQ=255)
                                                                                                                                                                                          tGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTccc                                                                                                                                               >  1:654066/1‑86 (MQ=255)
                                                                                                                                                                                           gCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGTAACTTCGac                                                                >  2:1017523/1‑164 (MQ=255)
                                                                                                                                                                                               aGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTg                                                                                                                                                                                 <  1:1132635/48‑1 (MQ=255)
                                                                                                                                                                                               aGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAg                                                                                                                                              >  1:193682/1‑83 (MQ=38)
                                                                                                                                                                                                 ccGCCGGTAACCGTGCCCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGAc                                <  1:234641/190‑1 (MQ=255)
                                                                                                                                                                                                  cGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGac          <  1:1225827/212‑1 (MQ=255)
                                                                                                                                                                                                    ccGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTCAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCAcaa                           >  2:1013454/1‑191 (MQ=255)
                                                                                                                                                                                                    ccGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGGAATTTc                                                                   >  1:202174/1‑153 (MQ=255)
                                                                                                                                                                                                    ccGGTAACCGTGAACCTGTCCAACGGCCAAACCATTACCgg                                                                                                                                                                                   >  2:326101/1‑40 (MQ=255)
                                                                                                                                                                                                       agacccGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCg                                        <  1:142014/173‑1 (MQ=255)
                                                                                                                                                                                                         aaCCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTATAACCCAGGGCAGCGTCGATTTCCAGTCCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTggcggc                                                                         >  2:927946/1‑142 (MQ=255)
                                                                                                                                                                                                         aaCCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGAcc                                                                                                                                           >  1:641557/1‑76 (MQ=40)
                                                                                                                                                                                                         aaCCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACgga                                                                                                                >  1:294919/1‑102 (MQ=25)
                                                                                                                                                                                                         aaCCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACTGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAATTTCGAACAACTGACCCCGAACCCGACGCCGGCCCAGACCACGATCAACGACt                  >  2:1173711/1‑197 (MQ=255)
                                                                                                                                                                                                         aaCCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCCACaacaac                                                                                                                   >  2:1365617/1‑100 (MQ=255)
                                                                                                                                                                                                               gAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCt                                                                                                               >  1:771458/1‑98 (MQ=35)
                                                                                                                                                                                                               gAACCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCg                <  1:456095/193‑1 (MQ=18)
                                                                                                                                                                                                                 |                                                                                                                                                                                                            
CGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCAC  >  AE015451/195523‑195936

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