Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 201,007 T→C 6.5% G2171G (GGT→GGC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451201,0070TC 6.5% 225.0 / 2.4 71G2171G (GGT→GGCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (30/35);  new base C (3/2);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCT  >  AE015451/200789‑201237
                                                                                                                                                                                                                          |                                                                                                                                                                                                                                        
tcTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGGCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGt                                                                                                                                                                                                                                          <  1:82055/218‑1 (MQ=255)
 cTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGACGGCGGT‑‑GAGAGCa                                                                                                                                                                                                                                 >  1:323063/1‑225 (MQ=255)
        aaCCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTg                                                                                                                                                                                                                <  1:1256725/235‑1 (MQ=39)
            aGCGCGCCGTACACCCACGAAGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACcaa                                                                                                                                                                                                           <  1:793435/236‑1 (MQ=39)
              cgcgcCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCAc                                                                                                                                                                                                                                <  2:1110241/213‑1 (MQ=40)
                 gcCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGAt                                                                                                                                                                                                                                   >  2:821010/1‑207 (MQ=39)
                  ccGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTaca                                                                                                                                                                                                                           <  2:897129/214‑1 (MQ=39)
                      agacCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACCGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGC‑‑GAGATCACCTACACCATCACGCTGACcaa                                                                                                                                                                                                           >  1:643337/3‑226 (MQ=11)
                          ccACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTg                                                                                                                                                                                              <  1:1190641/235‑1 (MQ=255)
                           cACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTACGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCCAAGCTGACTGCCACTCCGTCGGTCACCGAAGGCGGC‑‑GAGATCACCTACACCATc                                                                                                                                                                                                                      <  2:643334/210‑1 (MQ=11)
                              gacgCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGt                                                                                                                                                                                                                                          <  2:896973/189‑1 (MQ=39)
                                 gcgcAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGC‑‑GAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATcaa                                                                                                                                                                                      <  1:1157023/236‑1 (MQ=255)
                                      aGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGt                                                                                                                                                                                                                                          <  2:833302/181‑1 (MQ=39)
                                             gACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGa                                                                                                                                                                          <  1:801501/236‑1 (MQ=40)
                                              aCGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTGGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATcaa                                                                                                                                                                                      <  2:549471/223‑1 (MQ=39)
                                                gTCTATCAGGATGCGGGCCAGACCAGCCTCGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCAc                                                                                                                                                                                                                                <  2:191004/179‑1 (MQ=255)
                                                     tCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTAGCCGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGt                                                                                                                                                                                                                                          <  2:1001568/166‑1 (MQ=39)
                                                             cGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCCGTCACCGAAGGCGGC‑‑GAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCaa                                                                                                                                                              >  1:527886/1‑232 (MQ=2)
                                                             cGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTg                                                                                                                                                          <  1:1009226/236‑1 (MQ=35)
                                                              gggCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTg                                                                                                                                                          >  2:942545/1‑235 (MQ=34)
                                                                  cAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGTAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCAc                                                                                                                                                                                                                                <  2:881898/161‑1 (MQ=18)
                                                                  cAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACt                                                                                                                                                                              >  2:799079/1‑211 (MQ=34)
                                                                       cAGCCTCGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTc                                                                                                                                                                             >  2:1206798/1‑207 (MQ=35)
                                                                          ccTCGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGcacc                                                                                                                                              <  1:922533/235‑1 (MQ=33)
                                                                            tCGGCATCACTTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAAcc                                                                                                                                                                                 <  2:528330/198‑1 (MQ=21)
                                                                                   cACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACt                                                                                                                                                                              <  1:256950/194‑1 (MQ=33)
                                                                                      gTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTg                                                                                                                                  <  1:825608/235‑1 (MQ=35)
                                                                                       tacGCTATGGATGTAGACGGTCACACCTTCGAGAACATGCAACTGGGCGGCAACGCTTCGGTGCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCt                                                                                                                                                                                               <  1:115679/171‑1 (MQ=16)
                                                                                       tCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCa                                                                                                                                                                                                                        >  1:219980/1‑148 (MQ=32)
                                                                                       tCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTcgtcgt                                                                                                                                   >  2:526163/1‑233 (MQ=33)
                                                                                       tCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTAGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTg                                                                                                                                  >  2:766817/1‑234 (MQ=21)
                                                                                                 aTGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGC‑‑GAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAg                                                                                                                       >  1:1223483/1‑235 (MQ=9)
                                                                                                          gtgacgccTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACTGAAGGCGGT‑‑GAGATCACCTaca                                                                                                                                                                                                                           <  2:732295/120‑1 (MQ=255)
                                                                                                           tCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTg                                                                                                                                                                      >  1:899707/1‑178 (MQ=33)
                                                                                                           tCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCt                                                                                                            >  1:547853/1‑236 (MQ=37)
                                                                                                             acacCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCt                                                                                                            >  2:736431/1‑234 (MQ=37)
                                                                                                              cacCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCAc                                                                                                                                                                                                                                >  1:754532/1‑117 (MQ=33)
                                                                                                              cacCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCAc                                                                                                                                                                                                                                <  2:754529/117‑1 (MQ=33)
                                                                                                              cacCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACCGGCTCgg                                                                                                         >  1:938820/1‑236 (MQ=18)
                                                                                                               tccTTCGAGAACCTGCAACTGGGCGGCCACGCGTAGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCAACGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGcacc                                                                                                                                              <  1:282066/197‑1 (MQ=16)
                                                                                                                 gttCGAGAACCTGGAGCTGGGCGGCTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATc                                                                                                                                                                                                                      <  1:1136022/123‑1 (MQ=11)
                                                                                                                  ttCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGcacc                                                                                                                                              <  1:431515/195‑1 (MQ=33)
                                                                                                                    cGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTg                                                                                                   <  1:766820/236‑1 (MQ=39)
                                                                                                                        aaCCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTa                                                                                                                                                                    >  1:1088293/1‑167 (MQ=33)
                                                                                                                        aaCCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACggg                                                                                                              >  1:1083319/1‑221 (MQ=35)
                                                                                                                               atcTGGGCGGCAACGCTTCGGTTCAAGTGCCCGACGCCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCa                                                                                        <  1:369737/234‑1 (MQ=255)
                                                                                                                                aCTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAg                                                                                       >  1:1061980/1‑236 (MQ=40)
                                                                                                                                       ggcAACGCGTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCa                                                                                        <  2:651492/228‑1 (MQ=255)
                                                                                                                                        gcAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGG‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTc                                                                                                                                                                             >  2:213121/1‑142 (MQ=21)
                                                                                                                                           aCGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTa                                                                             <  1:65484/235‑1 (MQ=40)
                                                                                                                                             gCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGCTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGTTCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAg                                                                                       <  1:1207438/223‑1 (MQ=255)
                                                                                                                                                          tcaCCGACACCCGCGACGAGGTGGTGGCCAAGCTGACCGCCACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTgtcg                                                                                                                                       <  2:28113/160‑1 (MQ=255)
                                                                                                                                                          tGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGTGAGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCGGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACggg                                                                                                              >  2:1212695/1‑189 (MQ=255)
                                                                                                                                                               gACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCggt                                                        <  1:942548/236‑1 (MQ=255)
                                                                                                                                                                 cacCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCCACAGCACCACCGGCTCGGCCACTGCAACCGCAc                                                                                          >  1:1079491/1‑200 (MQ=17)
                                                                                                                                                                         acgaGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAAc                                                                                               >  2:995535/1‑187 (MQ=39)
                                                                                                                                                                         acgaGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCa                                                                                        <  2:1042438/194‑1 (MQ=255)
                                                                                                                                                                         acgaGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGc                                               >  1:287096/1‑235 (MQ=255)
                                                                                                                                                                                 gtggCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAAc                                                                                               >  2:1046273/1‑179 (MQ=40)
                                                                                                                                                                                 gtggCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTcc                                                                              >  1:1164526/1‑195 (MQ=255)
                                                                                                                                                                                  tggCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTc                                                                                                                                                                >  1:1133988/1‑113 (MQ=35)
                                                                                                                                                                                   ggCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATTACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCTGCACTAATGCGCCGGTGGTCAATGCCATCGACGcc                                     >  1:383812/1‑235 (MQ=255)
                                                                                                                                                                                    gCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGt                                   <  1:963568/236‑1 (MQ=255)
                                                                                                                                                                                     cGACGCTGACCGCCACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTg                                                                                                                                  <  1:110749/140‑1 (MQ=255)
                                                                                                                                                                                          cTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCaa                                                                                                                          <  2:1097164/143‑1 (MQ=37)
                                                                                                                                                                                           tGACCGCGACCCCTTCAGTCACCGAAGGCGGt                                                                                                                                                                                                                                          <  2:232480/32‑1 (MQ=40)
                                                                                                                                                                                            gACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCaa                                                                                                >  2:829064/1‑167 (MQ=255)
                                                                                                                                                                                                 cgACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAGCCGcca                                                                                          >  2:635780/1‑166 (MQ=255)
                                                                                                                                                                                                 cgACCCATTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCAcc                                                                                                                                                                                                                   <  1:447346/47‑1 (MQ=255)
                                                                                                                                                                                                          cAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATGAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAAc                                                                                  >  1:689691/1‑167 (MQ=255)
                                                                                                                                                                                                            gTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACg                                       >  2:394938/1‑208 (MQ=255)
                                                                                                                                                                                                             tCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTa                                                                             <  2:1233618/169‑1 (MQ=37)
                                                                                                                                                                                                               aCCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTg                                                                                                                                                                                              <  2:255959/54‑1 (MQ=35)
                                                                                                                                                                                                               aCCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGtt           <  2:547850/233‑1 (MQ=255)
                                                                                                                                                                                                                   aaGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGGGGaa              >  1:278064/1‑226 (MQ=255)
                                                                                                                                                                                                                     ggCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAaacct  >  1:449476/1‑236 (MQ=255)
                                                                                                                                                                                                                          |                                                                                                                                                                                                                                        
GCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGT‑‑GAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCT  >  AE015451/200789‑201237

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