Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 210,337 C→T 6.2% G5281G (GGC→GGT PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451210,3370CT 6.2% 244.7 / 2.8 79G5281G (GGC→GGTPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (38/36);  new base T (3/2);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

CCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGAC  >  AE015451/210106‑210571
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                          
ccACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGc                                                                                                                                                                                                                                        <  1:1287126/236‑1 (MQ=255)
ccACGCTGCGCAGGGTGATGACGTTTACAACGACGCGGGCCAGATCAGCCTGGGCATCACCTCAGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGc                                                                                                                                                                                                                                        <  1:807970/236‑1 (MQ=255)
           aGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCaacaac                                                                                                                                                                                                                             >  1:462843/1‑236 (MQ=255)
             ggTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGTACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGCTCAACAAcc                                                                                                                                                                                                                            <  2:677318/235‑1 (MQ=255)
             ggTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACca                                                                                                                                                                                                                           >  1:674562/1‑236 (MQ=255)
             ggTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACca                                                                                                                                                                                                                           >  1:268046/1‑236 (MQ=255)
              gTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAAcc                                                                                                                                                                                                                            <  2:1078959/234‑1 (MQ=255)
                    aCGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGc                                                                                                                                                                                                                     <  1:984261/235‑1 (MQ=255)
                          aCAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGAc                                                                                                                                                                                                               >  1:1255485/1‑235 (MQ=255)
                             cagacgCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATc                                                                                                                                                                                                                                   <  1:623192/210‑1 (MQ=255)
                               gacgCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCAc                                                                                                                                                                                                         <  1:709982/236‑1 (MQ=255)
                                              aGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGc                                                                                                                                                                                            <  2:1072449/234‑1 (MQ=255)
                                                   gggCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCa                                                                                                                                                                                           >  1:68664/1‑230 (MQ=255)
                                                       aTCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGtcat                                                                                                                                                                                 >  1:765838/1‑236 (MQ=255)
                                                        tCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGtcat                                                                                                                                                                                 >  1:992518/1‑235 (MQ=255)
                                                          aCCGCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGAcc                                                                                                                                                                                                              <  1:230286/204‑1 (MQ=255)
                                                            cTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGa                                                                                                                                                                                                <  1:397506/216‑1 (MQ=255)
                                                             tCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGt                                                                                                                                                                           >  2:886847/1‑236 (MQ=19)
                                                             tCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTg                                                                                                                                                                                                                      >  1:1096873/1‑193 (MQ=255)
                                                               ggcggTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGACTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGTCCGTCATCAcc                                                                                                                                                                             <  1:1206565/232‑1 (MQ=255)
                                                                gcggTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGcc                                                                                                                                                                        <  1:129651/236‑1 (MQ=255)
                                                                 cggTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGATCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCg                                                                                                                                                                       >  1:845671/1‑236 (MQ=255)
                                                                      gACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCa                                                                                                                                                                   <  1:1152629/235‑1 (MQ=255)
                                                                         gCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGGGCGGCTTCGGGTCAAGTGACCGACACCACCGACGAAGTGGGGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCAGc                                                                                                                                                                                                                                        >  2:397501/1‑163 (MQ=255)
                                                                                 tGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGTGGCGAGATCACTTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGTGCCGGCAAAACCGtcag                                                                                                                                                                                 >  1:355995/1‑209 (MQ=255)
                                                                                   cGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGc                                                                                                                                                                                            <  2:471245/197‑1 (MQ=255)
                                                                                        ttCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATc                                                                                                                                                                                                                                   >  2:189126/1‑153 (MQ=255)
                                                                                        ttCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGt                                                                                                                                                         <  1:8465/227‑1 (MQ=255)
                                                                                            agaACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGc                                                                                                                                                                     <  2:52685/211‑1 (MQ=255)
                                                                                             gaACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCg                                                                                                                                           >  2:996949/1‑236 (MQ=16)
                                                                                              aaCCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCCCCGAAGGCGGCGAGATCACCTACACCCTCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCa                                                                                                                                                                                                          <  1:393690/172‑1 (MQ=255)
                                                                                                ccTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCTTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCt                                                                                                                                                                                                                                          <  1:261249/138‑1 (MQ=255)
                                                                                                 cTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTg                                                                                                                                                                          <  2:189458/201‑1 (MQ=255)
                                                                                                  tGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACt                                                                                                                                      <  1:534898/236‑1 (MQ=255)
                                                                                                   ggAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCgg                                                                                                                                                                      >  1:389847/1‑203 (MQ=255)
                                                                                                   ggAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGc                                                                                                                                                                         >  2:777971/1‑200 (MQ=255)
                                                                                                     aGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCAACACGCTGACCAACAAAGACGGc                                                                                                                                                                                                                                            >  2:1052690/1‑131 (MQ=255)
                                                                                                       cTGGGCGGTGCGGCTTCGGTTCAAGAGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTAGGTCACCGAAGGCGGGGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCggt                                                                                                                                                                                                                                    >  2:798339/1‑135 (MQ=255)
                                                                                                          ggCGGTGCCGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCg                                                                                                                                <  1:7955/234‑1 (MQ=255)
                                                                                                             ggTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATcaacaa                                                                                                                                                                                                                              <  1:1158476/137‑1 (MQ=255)
                                                                                                             ggTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATcaacaa                                                                                                                                                                                                                              >  2:1158471/1‑137 (MQ=255)
                                                                                                             ggTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGCCATCACCg                                                                                                                                                                            >  2:242521/1‑187 (MQ=255)
                                                                                                              gTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGATCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTGACTGCCCCGGACAAc                                                                                                                          >  1:12450/1‑236 (MQ=37)
                                                                                                                    cTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATc                                                                                                                                                                                                                                   <  1:1013865/125‑1 (MQ=40)
                                                                                                                     ttCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAGGCTGACTGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATcaacaa                                                                                                                                                                                                                              <  1:179361/129‑1 (MQ=255)
                                                                                                                                tGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACg                                                                                                                                                                                              <  2:187117/150‑1 (MQ=37)
                                                                                                                                       caccacCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGAccc                                                                                                   >  2:1069936/1‑234 (MQ=255)
                                                                                                                                          cctcgacgaAGTGGTGGCCAAACTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGc                                                                                                                                                                                                                                        <  2:627246/95‑1 (MQ=11)
                                                                                                                                             cgacgaAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACg                                                                                                                                                                >  1:470653/1‑167 (MQ=38)
                                                                                                                                              gacgaAGTGGTGGCCAAACTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGc                                                                                                                                                                                                                     <  2:522589/113‑1 (MQ=255)
                                                                                                                                                    gtggtgGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGc                                                                                                                                                <  2:379581/176‑1 (MQ=38)
                                                                                                                                                    gtggtgGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGa                                                                                                                                        >  2:1271879/1‑184 (MQ=255)
                                                                                                                                                      ggtggCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGa                                                                                     <  2:992513/233‑1 (MQ=255)
                                                                                                                                                             aaGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCg                                                                                                                                           >  2:526559/1‑172 (MQ=40)
                                                                                                                                                                     tGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCCCCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTa                                                                                                                     >  2:122288/1‑186 (MQ=255)
                                                                                                                                                                            ccGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACg                                                                                                       >  1:70346/1‑193 (MQ=40)
                                                                                                                                                                            ccGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACg                                                                                                       <  2:70345/193‑1 (MQ=40)
                                                                                                                                                                            ccGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCaa                                                                                                                                                                  >  2:1251179/1‑134 (MQ=18)
                                                                                                                                                                             cGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTc                                                           >  1:773319/1‑236 (MQ=255)
                                                                                                                                                                                cGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTg                                                                                                                                                                                                     >  2:138683/1‑95 (MQ=18)
                                                                                                                                                                                     aCCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATcaacaa                                                                                                                                                                                                                              <  2:65587/65‑1 (MQ=40)
                                                                                                                                                                                       cGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCATAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGtt                                                   <  2:674741/234‑1 (MQ=255)
                                                                                                                                                                                                       aCCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAAcc                                                                                                                                                                                                                            <  2:68025/49‑1 (MQ=38)
                                                                                                                                                                                                        ccTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTgg                                  <  2:1215776/234‑1 (MQ=255)
                                                                                                                                                                                                               cATCACCCTGAACAACAACGGCGGCCTGCCGAGCAACAACCACAGTCCATTGAGCTTCACGCTGAGCTACGGCAATACCGTCATCACCGTGCCGGCCAACGGCAcc                                                                                                                                                           >  2:899609/1‑106 (MQ=11)
                                                                                                                                                                                                                aTCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGc                                              <  2:560440/214‑1 (MQ=255)
                                                                                                                                                                                                                  cACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGt                                                                           >  1:686810/1‑183 (MQ=40)
                                                                                                                                                                                                                    cccTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCAcc                                                                                                                                                           <  1:734547/101‑1 (MQ=35)
                                                                                                                                                                                                                    cccTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCg                                                                                                                                           >  1:146965/1‑117 (MQ=38)
                                                                                                                                                                                                                    cccTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCg                                                                                                                                           <  2:146962/117‑1 (MQ=38)
                                                                                                                                                                                                                    cccTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGt                     >  1:664271/1‑235 (MQ=255)
                                                                                                                                                                                                                    cccTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCCTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCAcc                                                                                                                                                           >  2:734542/1‑101 (MQ=25)
                                                                                                                                                                                                                      cTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAg                  >  1:893287/1‑236 (MQ=255)
                                                                                                                                                                                                                      cTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACg                                                                                                       >  1:604545/1‑151 (MQ=25)
                                                                                                                                                                                                                       tGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGCGGGCACTGCAACCGTGACTGCCCCGGACAACGCCTACGTCGGCAGCACCGCCCCTGTCgcg                                                                                            >  2:201954/1‑159 (MQ=255)
                                                                                                                                                                                                                          ccAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTc                                                                                                                 >  2:621909/1‑137 (MQ=37)
                                                                                                                                                                                                                            aacaaAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATGACCGTGCCGGCCAACgg                                                                                                                                                               >  2:307483/1‑89 (MQ=18)
                                                                                                                                                                                                                            aacaaAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTgg                                                                        >  2:488633/1‑176 (MQ=255)
                                                                                                                                                                                                                                      gCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGCCCCTGTCGTGATGt                                                                                        >  1:41139/1‑150 (MQ=255)
                                                                                                                                                                                                                                      gCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACcgac  <  1:501689/236‑1 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                          
CCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGAC  >  AE015451/210106‑210571

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: