Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 203,674 T→C 7.3% D3060D (GAT→GAC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451203,6740TC 7.3% 293.2 / 6.4 91D3060D (GAT→GACPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (48/36);  new base C (2/5);  total (50/41)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACTACCGTCACC  >  AE015451/203454‑203908
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                          
cgcAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATcaa                                                                                                                                                                                                                              >  1:686894/1‑235 (MQ=39)
   aGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACATCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGTTCACTTACACCATCACGCTGACCAACAAAGATGGTCTGCCGCTCaacaac                                                                                                                                                                                                                          >  1:645269/1‑236 (MQ=255)
       tgaggaCGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCGCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAAcc                                                                                                                                                                                                                         <  1:162937/233‑1 (MQ=255)
             cGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGc                                                                                                                                                                                                                  <  1:28866/234‑1 (MQ=38)
             cGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCCATAACCACAGTGCAt                                                                                                                                                                                                                >  1:162839/1‑236 (MQ=255)
                 tACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGAc                                                                                                                                                                                                            <  1:739781/236‑1 (MQ=40)
                  aCAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGAc                                                                                                                                                                                                            >  1:885102/1‑235 (MQ=40)
                  aCAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGAc                                                                                                                                                                                                            >  1:1249343/1‑235 (MQ=40)
                    aaCGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCtt                                                                                                                                                                                                         >  1:923409/1‑236 (MQ=39)
                               ccAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGa                                                                                                                                                                                                 <  2:885099/233‑1 (MQ=40)
                                      aGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCa                                                                                                                                                                                        <  1:68122/235‑1 (MQ=38)
                                                tCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGc                                                                                                                                                                                                                  <  2:896059/199‑1 (MQ=255)
                                                tCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGc                                                                                                                                                                                                                  >  1:896062/1‑199 (MQ=255)
                                                tCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGAc                                                                                                                                                                                            >  1:695036/1‑221 (MQ=40)
                                                tCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGGCTGCCGATCAACAAcc                                                                                                                                                                                                                         >  1:358293/1‑192 (MQ=38)
                                                     tCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCt                                                                                                                                                                                                                                       >  1:823646/1‑173 (MQ=38)
                                                     tCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTAACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCg                                                                                                                                                                         >  1:1187719/1‑235 (MQ=37)
                                                     tCGGCGGTGGACGCCACTGGTGCAACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTg                                                                                                                                                                                                  >  2:1154029/1‑210 (MQ=17)
                                                         cggTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCg                                                                                                                                                                    >  1:1137429/1‑236 (MQ=37)
                                                             ggACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCa                                                                                                                                                                <  1:254978/236‑1 (MQ=37)
                                                                   cACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCt                                                                                                                                                                                                          >  1:874617/1‑188 (MQ=39)
                                                                         tGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGc                                                                                                                                                                                                                                     >  1:1023789/1‑155 (MQ=39)
                                                                               cTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCaa                                                                                                                                                                                       <  1:47700/195‑1 (MQ=38)
                                                                                ttCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCAcc                                                                                                                                                                          >  2:898847/1‑207 (MQ=38)
                                                                                ttCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCggtctcgt                                                                                                                                                                  >  1:288075/1‑208 (MQ=255)
                                                                                accGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGACGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCg                                                                                                                                                                                                                                   <  2:819647/148‑1 (MQ=255)
                                                                                   gagaACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGCCCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCAcc                                                                                                                                                                          <  1:1027056/204‑1 (MQ=21)
                                                                                   gagaACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGc                                                                                                                                                                                         <  2:566056/189‑1 (MQ=39)
                                                                                      aaCCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGt                                                                                                                                       >  1:1278171/1‑236 (MQ=40)
                                                                                             aGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGc                                                                                                                                             >  1:97844/1‑223 (MQ=37)
                                                                                              gCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGccc                                                                                                                               >  1:907099/1‑236 (MQ=2)
                                                                                                 gggCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAAcc                                                                                                                                                                                                                         <  2:301975/143‑1 (MQ=39)
                                                                                                     ggTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGAc                                                                                                                          >  1:702533/1‑234 (MQ=40)
                                                                                                          cgcgAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCaaa                                                                                                                                                                                      >  2:286792/1‑169 (MQ=38)
                                                                                                          cgcgAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTc                                                                                                                    <  1:798009/235‑1 (MQ=16)
                                                                                                              aaGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACg                                                                                                                <  1:723885/235‑1 (MQ=255)
                                                                                                                 gTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTa                                                                                                                  <  1:1242025/230‑1 (MQ=255)
                                                                                                                  tGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTc                                                                                                              <  2:1278168/233‑1 (MQ=255)
                                                                                                                        tccccGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCTCTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATcaa                                                                                                                                                                                                                              <  1:281064/112‑1 (MQ=255)
                                                                                                                        tCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAAc                                                                                                                                                                                    <  2:1274991/157‑1 (MQ=255)
                                                                                                                        tCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAAc                                                                                                                                                                                    >  1:1274994/1‑157 (MQ=39)
                                                                                                                         ccccGCCACCACCGACGAAGTCGTGGCCACGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGc                                                                                                                                                                                         <  2:288074/149‑1 (MQ=255)
                                                                                                                         cACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCAcc                                                                                                                                                        <  2:644295/184‑1 (MQ=39)
                                                                                                                         cACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCAcc                                                                                                                                                        >  1:644298/1‑184 (MQ=39)
                                                                                                                             gACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGa                                                                                                                                     >  1:258704/1‑199 (MQ=255)
                                                                                                                              acacCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACg                                                                                                                                                                                           >  2:1283499/1‑144 (MQ=39)
                                                                                                                                  cgccGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGc                                                                                                                                                                                                                                     <  2:361352/96‑1 (MQ=255)
                                                                                                                                       acgaAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCg                                                                                                                                        <  1:309079/186‑1 (MQ=40)
                                                                                                                                       acgaAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCg                                                                                                                                        >  2:309078/1‑186 (MQ=40)
                                                                                                                                       acgaAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGc                                                                                                                                                                                                                  <  2:455996/112‑1 (MQ=12)
                                                                                                                                             tCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAatcgatc                                                                                <  1:391725/236‑1 (MQ=37)
                                                                                                                                               gTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTa                                                                                                                  <  2:60602/200‑1 (MQ=40)
                                                                                                                                               gTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGGGATCAAagcgatcg                                                                               >  1:146199/1‑235 (MQ=255)
                                                                                                                                                          gACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATcaaca                                                                                                                                                                                                                            >  1:630968/1‑83 (MQ=255)
                                                                                                                                                          gACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAg                                                                   >  1:205160/1‑236 (MQ=37)
                                                                                                                                                          gACCGCCACCCCGTCGATCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATcaaca                                                                                                                                                                                                                            <  2:630965/83‑1 (MQ=255)
                                                                                                                                                              gcgaCCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATc                                                                                         >  2:1017216/4‑210 (MQ=25)
                                                                                                                                                                cACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCg                                                              <  1:788925/235‑1 (MQ=37)
                                                                                                                                                                    ccGTCGGTCACCGAAGGTGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCggac                                                                                                                                                                      >  2:288491/1‑124 (MQ=25)
                                                                                                                                                                             aCCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCaa                                                                                                                                                               >  2:904557/1‑125 (MQ=35)
                                                                                                                                                                                  aGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGc                                           <  1:898850/236‑1 (MQ=38)
                                                                                                                                                                                      ggcgAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGt                                                                                                           >  1:621151/1‑168 (MQ=38)
                                                                                                                                                                                       gcgAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGCAACTg                                      <  1:1121313/236‑1 (MQ=37)
                                                                                                                                                                                        cgAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGAACCAACGACCCTGCGACCAACTCGATCGCAACCGTTGAAGGt                                                                 >  1:295359/1‑208 (MQ=255)
                                                                                                                                                                                        cgAGAGCACTTACACAATCACCCTGACCAACAAAGATGGTCTGCCTATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGTCAAAACCGTCATCACCGTGCCGGACAACGGCACCGTGGGCa                                                                                                                                                 <  1:80906/128‑1 (MQ=11)
                                                                                                                                                                                                    cacCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAg                          <  1:1040084/235‑1 (MQ=14)
                                                                                                                                                                                                    cacCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAg                          <  1:522224/235‑1 (MQ=14)
                                                                                                                                                                                                     acCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAg                          <  2:1117331/234‑1 (MQ=14)
                                                                                                                                                                                                       cATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACg                       <  1:1275008/235‑1 (MQ=37)
                                                                                                                                                                                                        attaCCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACg                       <  1:880615/231‑1 (MQ=18)
                                                                                                                                                                                                            cccTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGtcatca                                                                                                                                                                            >  2:582590/1‑81 (MQ=35)
                                                                                                                                                                                                            cccTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGcccc                                                                                                                              >  1:861862/1‑127 (MQ=38)
                                                                                                                                                                                                            cccTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGcccc                                                                                                                              <  2:861859/127‑1 (MQ=39)
                                                                                                                                                                                                            cccTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCACCGACCCTGTGATCAAatcgatcg                                                                               >  2:228746/1‑174 (MQ=25)
                                                                                                                                                                                                            cccTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGtt                                                                       >  1:240678/1‑182 (MQ=35)
                                                                                                                                                                                                            cccTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTACGtt                                                >  2:278919/1‑205 (MQ=25)
                                                                                                                                                                                                            cccTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTc                 >  1:708896/1‑236 (MQ=37)
                                                                                                                                                                                                            cccTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGGACCAACGACCCtgt                                                                                             >  2:509980/1‑160 (MQ=37)
                                                                                                                                                                                                             ccTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCACCGCCCCt                                                                                               >  2:201537/1‑157 (MQ=25)
                                                                                                                                                                                                             ccTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGCACAACGTCTACGTCGGTACCAACGACCCtgc                                                                                             >  2:10059/1‑158 (MQ=25)
                                                                                                                                                                                                             ccTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCAAGATGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCa                <  1:685254/236‑1 (MQ=37)
                                                                                                                                                                                                              cTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGCTCGCAACCGCTGACGGGCCGGAc                                                           >  2:330088/1‑192 (MQ=255)
                                                                                                                                                                                                               tGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGc              <  1:228746/236‑1 (MQ=35)
                                                                                                                                                                                                                    aacaaAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCg                                                                                                                             >  1:1140190/1‑120 (MQ=39)
                                                                                                                                                                                                                    aacaaAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCg                                                                                                                             <  2:1140187/120‑1 (MQ=39)
                                                                                                                                                                                                                      caaAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGc                                                                                                                                 <  1:1122941/114‑1 (MQ=40)
                                                                                                                                                                                                                      caaAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGc                                                                                                                                 >  2:1122938/1‑114 (MQ=40)
                                                                                                                                                                                                                        aaGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAg                          >  2:477265/1‑215 (MQ=37)
                                                                                                                                                                                                                         aGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCg                                              >  1:1139389/1‑194 (MQ=37)
                                                                                                                                                                                                                          gATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACTACCGTCAc   >  1:1087510/1‑236 (MQ=34)
                                                                                                                                                                                                                           aTGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACTACCGTCAcc  <  1:592482/236‑1 (MQ=33)
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                          
CGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACTACCGTCACC  >  AE015451/203454‑203908

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