Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 203,547 A→C 12.1% E3018A (GAG→GCG)  PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451203,5470AC12.1% 196.4 / 8.0 74E3018A (GAG→GCG) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base A (29/35);  new base C (2/8);  total (31/43)
Fisher's exact test for biased strand distribution p-value = 1.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTG  >  AE015451/203325‑203776
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                     
agaACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGGTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGc                                                                                                                                                                                                                           <  1:1098854/235‑1 (MQ=17)
                                aaCCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACcgacga                                                                                                                                                                                          >  1:795651/1‑236 (MQ=37)
                                                                     gCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACcg                                                                                                                                                                                              >  2:445770/1‑195 (MQ=37)
                                                                               ccAGGTCACCATCAAGGCTGGCGAAACCTGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGTCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACTTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAggcggc                                                                                                                                            <  1:1074157/235‑1 (MQ=21)
                                                                                cAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACGCCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAggcggc                                                                                                                                            <  1:281263/234‑1 (MQ=18)
                                                                                      aCCATCGAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCAc                                                                                                                                    <  1:674439/236‑1 (MQ=18)
                                                                                        cATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACtt                                                                                                                                  <  1:1038758/236‑1 (MQ=35)
                                                                                            aaGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACGACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGCTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACcga                                                                                                                                                                                             >  2:608496/1‑173 (MQ=25)
                                                                                            aaGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTaca                                                                                                                               <  1:566706/235‑1 (MQ=37)
                                                                                            aaGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGACGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTaca                                                                                                                               <  1:343752/235‑1 (MQ=25)
                                                                                                ccggCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCCTCGAGAACCTGGAGCCGGGCGGTGCCGCGAAGGGGGATGTCACCGAcacccc                                                                                                                                                                                                >  2:149147/3‑164 (MQ=255)
                                                                                                cTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTacac                                                                                                                              >  2:685333/1‑232 (MQ=35)
                                                                                                cTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCAt                                                                                                                           >  1:443173/1‑235 (MQ=34)
                                                                                                cTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATc                                                                                                                          <  1:300054/236‑1 (MQ=34)
                                                                                                     gAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGGGCGACCTTCGAGAACCTGGAGCTGGGCGGGGCCGCGCAGGTggtt                                                                                                                                                                                                              >  2:207072/1‑145 (MQ=255)
                                                                                                       aaCCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGa                                                                                                                   >  1:903229/1‑236 (MQ=37)
                                                                                                                       acCTACGCTGCGCAGGGTGATGACGTTTGCAACGACGCTGTCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTCTATGTCACCGACACCACCGACGAAGTCGTg                                                                                                                                                                                   <  2:599476/156‑1 (MQ=255)
                                                                                                                       acCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCt                                                                                                   <  1:589830/236‑1 (MQ=38)
                                                                                                                         ccACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATGACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACACAGATGGTCTGc                                                                                                 <  1:1062234/236‑1 (MQ=255)
                                                                                                                         ccACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGc                                                                                                 <  1:309480/236‑1 (MQ=37)
                                                                                                                          cACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCAc                                                                                                                                                                                                         <  1:218871/131‑1 (MQ=34)
                                                                                                                              cTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGTTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCGACCCCTTCgg                                                                                                                                                          >  1:445027/1‑174 (MQ=16)
                                                                                                                               tGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGTTGACCGCGACCCCTTCGGTGACCGAAggcg                                                                                                                                              >  1:673492/1‑185 (MQ=255)
                                                                                                                                 cgcAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATcaa                                                                                          <  1:842972/235‑1 (MQ=40)
                                                                                                                                 cgcAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCAACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATc                                                                                            <  2:608345/233‑1 (MQ=17)
                                                                                                                                 cgcAGGGTGATGACGGTTACAACGACGCTGGCCAGATCAGCCTGGGCATCGCCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACcg                                                                                                                                                                                              <  2:144247/135‑1 (MQ=21)
                                                                                                                                    aGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCgag                                                                                                                                         >  2:441838/1‑185 (MQ=33)
                                                                                                                                    aGGGTGATGACGTGTACAACGACGCGGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCg                                                                                                                                                                      >  1:599477/1‑156 (MQ=18)
                                                                                                                                     gggTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAAcc                                                                                     <  1:1092638/236‑1 (MQ=40)
                                                                                                                                      ggTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAAcc                                                                                     <  1:708474/235‑1 (MQ=38)
                                                                                                                                      ggTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGGGGCCa                                                                                                                                                                               >  2:74667/1‑145 (MQ=21)
                                                                                                                                         gatgaCGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGCGGAAGCCACTGGTGCGACCCTCGAGAACCTGGAGCTGGGCGGTGCcgc                                                                                                                                                                                                                        >  2:251516/1‑101 (MQ=255)
                                                                                                                                         gatgaCGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGAGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGGGCCGCGAAgg                                                                                                                                                                                                                   >  1:676967/1‑106 (MQ=255)
                                                                                                                                          atgaCGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGAGGTGGACGCAACTGGTGCGACCTTCGAGAACCTGGAGTTGGACGCTGCGGCGAAGGTGGATGTCACCGACACCACCGACGAAGTAGTGACCAAgat                                                                                                                                                                           >  1:839667/1‑143 (MQ=11)
                                                                                                                                             aCGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGc                                                                              <  1:879198/235‑1 (MQ=38)
                                                                                                                                              cGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGAt                                                                                                                                                                                                              >  2:1107680/1‑106 (MQ=35)
                                                                                                                                              cGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCgag                                                                                                                                         >  1:871426/1‑175 (MQ=33)
                                                                                                                                              cGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTg                                                                                                  >  2:1013779/1‑214 (MQ=38)
                                                                                                                                                   atcaggacgCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATcaa                                                                                          <  1:1096218/212‑1 (MQ=255)
                                                                                                                                                     aaCGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCAcc                                                                                                                       <  2:17540/186‑1 (MQ=17)
                                                                                                                                                                ccAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCAACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGc                                                                              <  1:788519/216‑1 (MQ=25)
                                                                                                                                                                      cAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCg                                                                                               <  1:559990/193‑1 (MQ=39)
                                                                                                                                                                      cAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGc                                                     <  1:1034498/235‑1 (MQ=38)
                                                                                                                                                                       aGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGa                                                                                              >  2:788518/1‑193 (MQ=25)
                                                                                                                                                                       aGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGc                                                     <  1:459707/234‑1 (MQ=40)
                                                                                                                                                                           tGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTaca                                                                                                                               <  1:12990/156‑1 (MQ=35)
                                                                                                                                                                                aTCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTg                                                                               >  1:858725/1‑199 (MQ=40)
                                                                                                                                                                                 tCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTaca                                                                                                                               >  2:919895/1‑150 (MQ=38)
                                                                                                                                                                                 tCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGGGGATGTCACCGACACCACCGACGAAGGCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCgagag                                                                                                                                       >  2:230604/1‑141 (MQ=17)
                                                                                                                                                                                     cTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGAt                                                                                             >  1:283502/1‑180 (MQ=39)
                                                                                                                                                                                      tCGGCGGTGGACGCCACTGGTGCAACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGt                                    >  1:447376/1‑236 (MQ=17)
                                                                                                                                                                                                aCGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACaaa                                                                                                           >  2:493096/1‑155 (MQ=37)
                                                                                                                                                                                                  gCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAAc                          >  1:1066662/1‑234 (MQ=37)
                                                                                                                                                                                                   ccACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTg                                   <  1:995143/224‑1 (MQ=21)
                                                                                                                                                                                                    tgcTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCAc                                                                                                                                    <  2:974201/124‑1 (MQ=17)
                                                                                                                                                                                                       tGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATcaaca                                                                                        >  2:925207/1‑167 (MQ=40)
                                                                                                                                                                                                       tGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATcaaca                                                                                        <  1:925208/167‑1 (MQ=255)
                                                                                                                                                                                                        ggTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATcaaca                                                                                        <  2:963631/166‑1 (MQ=255)
                                                                                                                                                                                                          tGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCAcc                    <  2:795650/232‑1 (MQ=39)
                                                                                                                                                                                                           gcaaCCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGt                  <  2:447375/230‑1 (MQ=37)
                                                                                                                                                                                                            cGACCTTCGTGAACCTTGAGCTGGGCGGTGCCGCGAAGGTGGATGTCATCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCg                                                                                               <  2:1048370/155‑1 (MQ=255)
                                                                                                                                                                                                            cGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCAcc                                                                                                                       <  2:1093300/131‑1 (MQ=35)
                                                                                                                                                                                                            cGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCAt                                                                            >  1:747637/1‑174 (MQ=39)
                                                                                                                                                                                                            cGACCTCCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGc                                                                                                 <  2:1090841/153‑1 (MQ=21)
                                                                                                                                                                                                            cGACCATCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCCTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGc                                                                                                 <  2:627524/153‑1 (MQ=255)
                                                                                                                                                                                                             gACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTTCATGTCACCGACACCACCGACGAAGTCGTGGCCAAGCTGACCGCCACCCCGTCGGTCACCGAAGGTGGCGAGATCACTTACACCATCACCCTGACCAACAAAGACGGCCTGCCg                                                                                               <  2:549202/154‑1 (MQ=1)
                                                                                                                                                                                                                 ttCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGc                                                                                                 <  2:958737/148‑1 (MQ=25)
                                                                                                                                                                                                                 ttCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCAc                                                                                                                                    <  2:541364/113‑1 (MQ=32)
                                                                                                                                                                                                                 ttCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCAc                                                                                                                                    >  1:541365/1‑113 (MQ=32)
                                                                                                                                                                                                                 ttCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATTACCCAGACCAACAAAg                                                                                                          >  2:647877/1‑139 (MQ=255)
                                                                                                                                                                                                                 ttCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGCCCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCaa                                                   <  2:1045379/194‑1 (MQ=38)
                                                                                                                                                                                                                       aaCCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCg    <  1:574277/235‑1 (MQ=25)
                                                                                                                                                                                                                         ccTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTg  <  1:932470/235‑1 (MQ=40)
                                                                                                                                                                                                                           gggAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGc              <  1:221361/220‑1 (MQ=37)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                     
AGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTG  >  AE015451/203325‑203776

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: