Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 202,117 T→C 11.0% S2541S (AGT→AGC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451202,1170TC11.0% 184.0 / 6.0 62S2541S (AGT→AGCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (31/23);  new base C (4/3);  total (35/27)
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.96e-01

AGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGA  >  AE015451/201882‑202346
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                     
aGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGt                                                                                                                                                                                                                                       <  1:612380/236‑1 (MQ=255)
                    gACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGAcgc                                                                                                                                                                                                                   <  1:872835/236‑1 (MQ=255)
                                 tCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGa                                                                                                                                                                                                       >  1:387525/1‑235 (MQ=255)
                                        tACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTgcgc                                                                                                                                                                                                                                   >  1:180295/1‑200 (MQ=255)
                                                cTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCacgac                                                                                                                                                                                                                     <  2:1074696/206‑1 (MQ=255)
                                                cTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCacgac                                                                                                                                                                                                                     >  1:1074699/1‑206 (MQ=255)
                                                     gAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCacgacg                                                                                                                                                                                                                    >  2:806626/1‑202 (MQ=40)
                                                     gAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCac                                                                                                                                                                                                                        >  2:1067328/1‑198 (MQ=255)
                                                        ccAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCGACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGAcgc                                                                                                                                                                                                                   >  2:803343/1‑200 (MQ=255)
                                                                  cAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCg                                                                                                                                                                                                           <  1:1010631/198‑1 (MQ=40)
                                                                   aGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTg                                                                                                                                                                                                                                      <  2:931189/170‑1 (MQ=255)
                                                                       aaCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCa                                                                                                                                                                 <  1:443028/235‑1 (MQ=255)
                                                                        aCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGATCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGAGAt                                                                                                                                                                                                         <  2:558340/194‑1 (MQ=255)
                                                                         cGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTg                                                                                                                                                                                                                                      >  2:815774/1‑164 (MQ=255)
                                                                            aGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTcc                                                                                                                                                            <  1:999568/235‑1 (MQ=255)
                                                                                gACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCtg                                                                                                                                                        <  1:1083090/235‑1 (MQ=255)
                                                                                 aCATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCtgtg                                                                                                                                                      >  1:685698/1‑236 (MQ=255)
                                                                                                 cATCACAGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACTATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTc                                                                                                                                         <  2:1036342/233‑1 (MQ=255)
                                                                                                    cACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCt                                                                                                                                   >  1:412754/1‑236 (MQ=255)
                                                                                                                              agaACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGACgcgc                                                                                                                                                                                                                 >  2:172985/1‑132 (MQ=255)
                                                                                                                               gaACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTg                                                                                                                        <  2:788869/220‑1 (MQ=255)
                                                                                                                                  cGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCtt                                                                                                     <  1:631455/236‑1 (MQ=255)
                                                                                                                                      aaGCCAACCGTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCGACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACg                                                                                                        <  1:803346/229‑1 (MQ=255)
                                                                                                                                        gCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTc                                                                                               <  1:612268/236‑1 (MQ=255)
                                                                                                                                                 cACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAgg                                                                                                                                                                                            <  2:716044/134‑1 (MQ=255)
                                                                                                                                                     gTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACg                                                                                                        <  2:576832/214‑1 (MQ=255)
                                                                                                                                                        cACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCAAGTCCGCTGTGGacgg                                                                                                                                                 >  2:924496/1‑167 (MQ=11)
                                                                                                                                                         aCGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTc                                                                              >  1:282280/1‑236 (MQ=255)
                                                                                                                                                         aCGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCGGTGGac                                                                                                                                                   >  2:363517/1‑166 (MQ=255)
                                                                                                                                                          cGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGac                                                                                    <  2:788840/229‑1 (MQ=255)
                                                                                                                                                          cGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGac                                                                                    >  1:788843/1‑229 (MQ=255)
                                                                                                                                                              caCCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCCAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTg                                                                                                                  <  2:76462/194‑1 (MQ=255)
                                                                                                                                                              aatcAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGCTCAAGTGACCGACACGCTcgac                                                                           >  2:374911/4‑234 (MQ=255)
                                                                                                                                                              aatcAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTca                                                                                                                                        >  2:821488/4‑173 (MQ=11)
                                                                                                                                                                   tgctcaggCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCGACAAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCtga                                                                                                                                                       >  1:106624/7‑152 (MQ=14)
                                                                                                                                                                   gCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGcc                                                                                                                                                                          >  2:412182/1‑134 (MQ=255)
                                                                                                                                                                            tCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTcg                                                                             >  1:442575/1‑218 (MQ=255)
                                                                                                                                                                                 aCCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGggcgg                                                                                                             >  2:150132/1‑181 (MQ=255)
                                                                                                                                                                                   cTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACcg                                                     <  1:821491/235‑1 (MQ=16)
                                                                                                                                                                                        cGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCCACCTGGGCATc                                                                                                                                                                  >  2:144071/1‑121 (MQ=255)
                                                                                                                                                                                               cTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGCCGGTCACACCTTCGAGAACCTGCAGCTGGGCGGCGACGCGTCg                                                                                                   >  2:1017371/1‑177 (MQ=255)
                                                                                                                                                                                                   gCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGGGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACC                                    >  1:170175‑M1/1‑236 (MQ=255)
                                                                                                                                                                                                    caacaaCGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCg                                                               <  2:999817/208‑1 (MQ=255)
                                                                                                                                                                                                    caacaaCGCTCAGGTCACCATCAAGGCCGGCGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGac                                                                                                                                                   >  1:586396/1‑123 (MQ=255)
                                                                                                                                                                                                     aacaacGCTCAGGTCACCATCAAGGCCGGTGAGACCAGt                                                                                                                                                                                                                                       <  1:915435/39‑1 (MQ=255)
                                                                                                                                                                                                       caacGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCG                                   <  2:945081‑M1/233‑1 (MQ=255)
                                                                                                                                                                                                        accgcTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAG                                <  1:111743‑M1/233‑1 (MQ=255)
                                                                                                                                                                                                             tCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTg                          <  1:740784/236‑1 (MQ=255)
                                                                                                                                                                                                             tCAGGTCACCATCAAGGCCGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCAACCAGAACGCTGGCCAGATCAGCCTGGGCATCACGTCGGCTGTGGATGTAGACGGTGACACCTTCGAGAACCTGCa                                                                                                                      <  1:960397/144‑1 (MQ=14)
                                                                                                                                                                                                               aGGTCACCATCAAGGCCGGTGATACCAGAGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGATGTAGACGGTGACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTcgacga                                                                         <  1:521662/187‑1 (MQ=14)
                                                                                                                                                                                                                    aCCATCAAGGGCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGa                                                                                                                                                    >  2:202091/1‑107 (MQ=255)
                                                                                                                                                                                                                       aTCAAGGCCGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCg                                                                                                                                <  1:914043/124‑1 (MQ=255)
                                                                                                                                                                                                                       aTCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCg                                                                                                                                >  2:914040/1‑124 (MQ=255)
                                                                                                                                                                                                                         cAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCCGCGGTGGATGCAGACGGTCACACCTTCGCGAACCTGCAACTg                                                                                                                  >  2:133123/1‑136 (MQ=14)
                                                                                                                                                                                                                              ccGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCAt          >  1:748041/1‑235 (MQ=255)
                                                                                                                                                                                                                              ccGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCA                                      <  1:1011711‑M1/207‑1 (MQ=255)
                                                                                                                                                                                                                              ccGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCA                                      >  2:1011708‑M1/1‑207 (MQ=255)
                                                                                                                                                                                                                                ggTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGcc                                                                                                                                                                          >  2:133147/1‑73 (MQ=40)
                                                                                                                                                                                                                                ggTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCAc       >  1:815830/1‑236 (MQ=255)
                                                                                                                                                                                                                                   gagaCCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGTAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGCCACCGAAGGCGGTGAGATCACCTACAcc            >  1:393369/1‑228 (MQ=255)
                                                                                                                                                                                                                                   gagaCCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCCACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGGGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGCGAGATCACCTacac             >  1:121969/1‑227 (MQ=255)
                                                                                                                                                                                                                                   gagaCCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGCGCATCACGTCCTCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTTGTGGCGAAGc                                                           >  1:907585/1‑181 (MQ=255)
                                                                                                                                                                                                                                     gaCCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCTGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGa  >  1:507490/1‑236 (MQ=255)
                                                                                                                                                                                                                                         aGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTga                         <  1:517828/209‑1 (MQ=40)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                     
AGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGA  >  AE015451/201882‑202346

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: