Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 213,079 C→T 8.4% G6195G (GGC→GGT PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451213,0790CT 8.4% 234.0 / 5.1 80G6195G (GGC→GGTPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (38/35);  new base T (2/5);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGC  >  AE015451/212849‑213311
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                        
ttCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACCCCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGc                                                                                                                                                                                                                                     <  1:232979/236‑1 (MQ=37)
ttCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACt                                                                                                                                                                                                                                       >  2:974506/1‑234 (MQ=39)
ttCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGc                                                                                                                                                                                                                                     <  1:418662/236‑1 (MQ=37)
    agaACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGc                                                                                                                                                                                                                                          <  1:1008640/227‑1 (MQ=38)
         cTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCg                                                                                                                                                                                                                                <  1:525304/232‑1 (MQ=38)
                    cGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACa                                                                                                                                                                                                                 >  1:170825/1‑236 (MQ=40)
                         ccGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTc                                                                                                                                                                                                            <  1:974509/236‑1 (MQ=255)
                            cgAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTAc                                                                                                                                                                                                         >  1:325236/1‑236 (MQ=255)
                              aaGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACg                                                                                                                                                                                                        <  1:373299/235‑1 (MQ=255)
                                     aTGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGc                                                                                                                                                                                                                                     <  1:1073667/199‑1 (MQ=25)
                                          cccGACACCACCGACGAAGTCGTGGCCCAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACg                                                                                                                                                                                                        <  1:251066/222‑1 (MQ=255)
                                              acacCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACt                                                                                                                                                                                                                           >  2:368656/1‑200 (MQ=14)
                                               caccacCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTg                                                                                                                                                                                      >  1:749401/1‑236 (MQ=255)
                                                accaccGACGAAGACGTGGCAAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGa                                                                                                                                                                                           <  1:308851/230‑1 (MQ=255)
                                                     cgacgaAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTg                                                                                                                                                                                 <  1:919200/235‑1 (MQ=255)
                                                          aaGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATg                                                                                                                                                                              >  1:726372/1‑233 (MQ=255)
                                                               gTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGAt                                                                                                                                                                         >  2:931965/1‑233 (MQ=17)
                                                               gTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGc                                                                                                                                                                      <  1:641854/236‑1 (MQ=255)
                                                                tGGCCAAGCTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACt                                                                                                                                                                                                                                       <  1:225623/170‑1 (MQ=14)
                                                                     aaGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCCCCGTCCCTCCCCCGGACAATGt                                                                                                                                                                                                             >  1:1074472/1‑191 (MQ=255)
                                                                          gACCGCCACCCCGTCGGTCACCGAGGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAAt                                                                                                                                                                                                               >  2:290415/1‑184 (MQ=255)
                                                                           cccGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAg                                                                                                                                                           <  1:205057/234‑1 (MQ=40)
                                                                           aCCGCCACCCCGTCGGTCACCGAAGACGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCg                                                                                                                                                                                                                                <  1:751711/166‑1 (MQ=25)
                                                                             cGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACCCCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCg                                                                                                                                                                                                                                <  2:214680/164‑1 (MQ=25)
                                                                             cGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGGCCCTGTCGTGATGTc                                                                                                                                                                            <  1:334014/216‑1 (MQ=39)
                                                                                cACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTg                                                                                                                                                                                                                          >  2:961158/1‑167 (MQ=39)
                                                                                cACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTg                                                                                                                                                                                                                          <  1:961161/167‑1 (MQ=39)
                                                                                            cACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACACTGTCCACGTCg                                                                                                                                                                                                     >  2:213788/1‑176 (MQ=255)
                                                                                                gAAGGCGGCGAGATTCCCTACACCAGCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAAt                                                                                                                                                                                                               <  2:301411/162‑1 (MQ=37)
                                                                                                       gcgAGATTACCTACACCATCACGCTGACGAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGCCCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTg                                                                                                                              <  1:1047413/236‑1 (MQ=255)
                                                                                                          agaTTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGAc                                                                                                                            <  1:724553/235‑1 (MQ=38)
                                                                                                          agaTTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGa                                                                                                                             >  2:539936/1‑234 (MQ=40)
                                                                                                            aTTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTg                                                                                                                              >  2:652053/1‑231 (MQ=40)
                                                                                                            aTTACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGc                                                                                                                                   <  1:1085684/226‑1 (MQ=39)
                                                                                                                    cacCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGa                                                                                                                 <  1:840682/236‑1 (MQ=255)
                                                                                                                       cATTACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCg                                                                                                                                                                                                                                <  1:773949/122‑1 (MQ=18)
                                                                                                                            cccTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCaa                                                                                                                                                                                              >  1:563753/3‑151 (MQ=2)
                                                                                                                            cccTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGCCACCGTGGGCACTGCCACCGTCTCTGCCCCGGACAATGTCTACGTCGGTACCat                                                                                                                                                                                              <  2:563750/149‑2 (MQ=11)
                                                                                                                            cGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCaa                                                                                                                                 >  1:1031828/1‑212 (MQ=39)
                                                                                                                            cGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGACACTGCCCCGGACAATGTCTACGTCGGTACCAACGCCCCTg                                                                                                                                                                                      >  2:98562/1‑159 (MQ=25)
                                                                                                                             gCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAATGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGCCAAGACGCCGGTc                                                                                                         >  1:170659/1‑235 (MQ=255)
                                                                                                                             gCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGAcac                                                                                                                          >  2:641805/1‑216 (MQ=38)
                                                                                                                                    aacaaAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAg                                                                                                                                          >  2:934403/1‑195 (MQ=40)
                                                                                                                                    aacaaAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGc                                                                                                      >  2:167164/1‑231 (MQ=255)
                                                                                                                                         aGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGc                                                                                                                                                                                                                                     >  1:574989/1‑99 (MQ=33)
                                                                                                                                         aGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCg                                                                                                                                                                                                                                <  1:980528/104‑1 (MQ=33)
                                                                                                                                              gCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGAcac                                                                                                                          >  2:834393/1‑199 (MQ=38)
                                                                                                                                                    cGATTAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAg                                                                                                                                                           <  1:432331/162‑1 (MQ=255)
                                                                                                                                                    cGATTAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAg                                                                                                                                                           >  2:432328/1‑162 (MQ=255)
                                                                                                                                                    cGATCAACATCCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCt                                                                                                                                                                                       <  1:1048423/134‑1 (MQ=21)
                                                                                                                                                            aacCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTAcccc                                                                         >  1:407248/1‑236 (MQ=16)
                                                                                                                                                                    tGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTggcc                                                                                                                     >  2:314572/1‑182 (MQ=255)
                                                                                                                                                                    tGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACg                                                                 <  1:421167/236‑1 (MQ=255)
                                                                                                                                                                              ttCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGATCGAGCAACTGACGCTgg                                                                                                                       >  2:215798/1‑172 (MQ=255)
                                                                                                                                                                                cACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCg                                                                                                                                      <  2:147680/155‑1 (MQ=255)
                                                                                                                                                                                 aCGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCg                                                                                                >  2:862823/1‑192 (MQ=255)
                                                                                                                                                                                 aCGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGt                                                    >  1:144544/1‑236 (MQ=255)
                                                                                                                                                                                   gCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGc                                                            <  2:638234/226‑1 (MQ=255)
                                                                                                                                                                                   gCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGc                                                            >  1:638237/1‑226 (MQ=255)
                                                                                                                                                                                   gCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGCTGGCAAGCTCGAGCACCTGACGCTGGa                                                                                                                      >  2:309859/1‑168 (MQ=255)
                                                                                                                                                                                          gACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGACGGCGCGGACGTTGGCAAGTTCGACCAACTGACTCTCGCCAAgccccc                                                                                                             >  2:94472/1‑165 (MQ=11)
                                                                                                                                                                                           aCGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACt                                                        <  2:941463/222‑1 (MQ=255)
                                                                                                                                                                                                aaaaCCGTCATCACCGTGCCGGCTAATGGCACCGTGGGCACTGCCACTGTCACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTgg                                                     <  2:405591/220‑1 (MQ=255)
                                                                                                                                                                                                aaaaCCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCa                                      >  1:912810/1‑235 (MQ=255)
                                                                                                                                                                                                  aaCCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCg                                                                                 <  2:325233/190‑1 (MQ=255)
                                                                                                                                                                                                    ccGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCa                                                                                                                                                                                               <  1:950603/78‑1 (MQ=39)
                                                                                                                                                                                                       tcatcaCCGTGCCGGCTAATGGCACCGTGGGTACTGCCACCGTCACTGCCCCGGATAACGTCTACGTCGGTAGCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTgg                                                                                                                       <  1:967286/147‑1 (MQ=11)
                                                                                                                                                                                                             ccGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACTTCGGTCAcc                                                                                          <  2:191405/170‑1 (MQ=255)
                                                                                                                                                                                                                tGCCGGCAAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCt                       <  1:13372/234‑1 (MQ=255)
                                                                                                                                                                                                                  ccGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGCACCCCAGGCAACGAAGGCGACTTGGTCAAGGGCACCATCACGGCCGACCAGACCTCGGt                   >  1:128319/1‑236 (MQ=255)
                                                                                                                                                                                                                   cGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACACTGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCg                                                                                 >  2:162599/1‑173 (MQ=255)
                                                                                                                                                                                                                     gCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAAcc                                                                                                                                                                  <  1:533955/90‑1 (MQ=255)
                                                                                                                                                                                                                        aaCGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCCACGACCCTGTCGTGCTGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCACGCTCGAGCAACTGACGCTGGACaa                                                                                                                   >  2:206486/1‑134 (MQ=255)
                                                                                                                                                                                                                        aaCGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCg              >  1:901422/1‑235 (MQ=255)
                                                                                                                                                                                                                               cggTGGGCACTGCCACCGTCCCTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGc                                                                                   <  1:107521/157‑1 (MQ=255)
                                                                                                                                                                                                                               ccGTGGGTACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGACCGCAACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACGc                                                                                                                          >  2:145960/1‑120 (MQ=255)
                                                                                                                                                                                                                                ggTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGc                                                            >  2:162721/2‑181 (MQ=255)
                                                                                                                                                                                                                                 gTGGGTACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTAGCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTc      <  2:454528/234‑1 (MQ=255)
                                                                                                                                                                                                                                  tGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGACAGCACGTCGGTGACcgacga                                                                                     >  2:497021/1‑154 (MQ=255)
                                                                                                                                                                                                                                   gggCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGCGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGc  >  1:1028760/1‑236 (MQ=255)
                                                                                                                                                                                                                                    ggCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGAcc                        <  1:167166/213‑1 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                        
TTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGC  >  AE015451/212849‑213311

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