Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 206,311 G→C 7.5% V3939V (GTG→GTC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451206,3110GC 7.5% 205.1 / 2.8 72V3939V (GTG→GTCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (29/36);  new base C (4/2);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTT  >  AE015451/206099‑206524
                                                                                                                                                                                                                    |                                                                                                                                                                                                                     
aTCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCa                                                                                                                                                                                                    >  1:557261/1‑232 (MQ=255)
                    cgcgcCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGgcggcg                                                                                                                                                                                                          <  2:452247/206‑1 (MQ=255)
                                ccACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGc                                                                                                                                                                <  1:433553/236‑1 (MQ=255)
                                      tGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGCCCGAAggcg                                                                                                                                                                                                             >  1:495524/1‑185 (MQ=255)
                                      tGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGc                                                                                                                                                                >  1:130420/1‑230 (MQ=255)
                                        cgcAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGt                                                                                                                                                                    >  2:313943/1‑224 (MQ=255)
                                         gcAAGGCGATGACGTCTACACCGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGc                                                                                                                                                                                     <  1:163635/206‑1 (MQ=255)
                                          cAAGGCGATGACGTCTACCACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCga                                                                                                                                                                                                         <  2:778722/185‑1 (MQ=255)
                                            aGGCGATGACGGCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTGCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAAcc                                                                                                                                                    <  1:159185/236‑1 (MQ=255)
                                                  tGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGt                                                                                                                                               <  1:453746/235‑1 (MQ=255)
                                                              cgacgCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGc                                                                                                                                                                <  2:601988/206‑1 (MQ=255)
                                                                 cgCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTc                                                                                                                                   >  2:729124/1‑232 (MQ=255)
                                                                   cTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACcaaca                                                                                                                                                                            >  2:633026/1‑189 (MQ=255)
                                                                         aGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGa                                                                                                                        >  1:735624/1‑235 (MQ=255)
                                                                          gATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACg                                                                                                                      <  1:333343/236‑1 (MQ=255)
                                                                                ccTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAggcggc                                                                                                                                                                                                           <  2:807252/145‑1 (MQ=255)
                                                                                   gggCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGCGAGCACCTACACCAGCAc                                                                                                                                                                                       >  2:167431/1‑162 (MQ=255)
                                                                                   gggCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCaaa                                                                                                                 >  2:766906/1‑232 (MQ=255)
                                                                                        tCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGCCGGTGACCGAAGGCGGCGAGATCACCACCACCATCACGCt                                                                                                                                                                                    >  2:28341/1‑160 (MQ=255)
                                                                                           cGTCGGCTGTGGATGTAGCCGGTCGCACCTTCGAGAACCTGCAGCTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGc                                                                                                                                             <  1:228433/196‑1 (MQ=255)
                                                                                             tCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATcaaca                                                                                                                                                       >  1:615421/1‑184 (MQ=255)
                                                                                             tCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATcaaca                                                                                                                                                       <  2:615419/184‑1 (MQ=255)
                                                                                             tCTGCTGTGGAGGTAGACGGTCGCACCGTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAggcggc                                                                                                                                                                                                           <  2:38906/132‑1 (MQ=255)
                                                                                                   gtgGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGc                                                                                             <  1:781964/236‑1 (MQ=255)
                                                                                                         gTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAAc                                                                                         <  2:572607/234‑1 (MQ=40)
                                                                                                         gTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACg                                                                                        >  1:541668/1‑235 (MQ=255)
                                                                                                         gTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCg                                                                                                                                                              >  2:12672/1‑165 (MQ=255)
                                                                                                          caGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGt                                                                                      <  1:451541/235‑1 (MQ=39)
                                                                                                            gACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTAc                                                                                    <  1:729126/236‑1 (MQ=255)
                                                                                                               ggTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGc                                                                                                                                             >  1:758154/1‑176 (MQ=39)
                                                                                                               ggTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGc                                                                                                                                             <  2:758152/176‑1 (MQ=38)
                                                                                                                gtaccACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAAcc                                                                                                                                                    <  1:556056/164‑1 (MQ=37)
                                                                                                                gTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCaacaac                                                                                                                                                     <  1:247375/167‑1 (MQ=39)
                                                                                                                  cGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGt                                                                                 <  2:735622/233‑1 (MQ=39)
                                                                                                                    cACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTaca                                                                                                                                                                                              >  2:769013/1‑122 (MQ=39)
                                                                                                                    cACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACg                                                                                                                                                                                      >  1:261322/1‑130 (MQ=39)
                                                                                                                      ccTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACCCCACCGACGAAGTTGTGGCCAAGCTCACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGcc                                                                                                <  1:140200/214‑1 (MQ=21)
                                                                                                                      ccTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGc                                                                                                                    <  1:648119/194‑1 (MQ=18)
                                                                                                                       cTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTa                                                                            >  2:764767/1‑233 (MQ=38)
                                                                                                                        ttCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTg                                                                                                  <  2:600500/210‑1 (MQ=35)
                                                                                                                             gaACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCg                                                                   <  1:458634/236‑1 (MQ=39)
                                                                                                                              acccTGCAACGGGGCGGTACCGCTACGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGc                                                                                                                                             <  2:158533/159‑1 (MQ=255)
                                                                                                                                       cTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCa                                                                                                                                            >  2:745610/1‑153 (MQ=35)
                                                                                                                                       cTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCAc                                                                                                                                 >  2:534628/1‑164 (MQ=37)
                                                                                                                                       cTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGCGCGACgg                                                                                                                     >  1:715654/1‑176 (MQ=25)
                                                                                                                                            cggcAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATc                                                                                                                                                                                                     >  1:648651/1‑91 (MQ=37)
                                                                                                                                             ggcAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCtt                                                                                                                                    >  2:40045/1‑155 (MQ=35)
                                                                                                                                             ggcAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGACaa                                                   >  1:706689/1‑236 (MQ=16)
                                                                                                                                             ggcAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCAGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTaat                                                                                         >  2:501839/1‑197 (MQ=19)
                                                                                                                                                  cGCTTCGGTTCAAGTCATCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACcaca                                                                                                                                                 <  2:636797/137‑1 (MQ=37)
                                                                                                                                                    cTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGTGGCGAGATCACTTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCaacaac                                                                                                                                                     <  2:557259/131‑1 (MQ=255)
                                                                                                                                                                tCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAAc                                                  <  2:639678/218‑1 (MQ=40)
                                                                                                                                                                 cACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCCACCGTCACTGCCCCg                                                        >  2:766213/1‑211 (MQ=21)
                                                                                                                                                                         ccaccGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAAc                                                                                         <  2:55382/170‑1 (MQ=33)
                                                                                                                                                                         ccaccGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGt                                          <  2:353293/217‑1 (MQ=39)
                                                                                                                                                                           accGACGAAGTGGTGGCCAAGCTGCCCGCGACACCGTCGGTGACCGAAggcggc                                                                                                                                                                                                           <  2:279860/54‑1 (MQ=11)
                                                                                                                                                                                cgaAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGa                   >  2:504680/1‑233 (MQ=255)
                                                                                                                                                                                 gaAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCAcc                                                                                                     >  2:614417/1‑150 (MQ=17)
                                                                                                                                                                                  aaGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCg                                                                   <  1:187421/183‑1 (MQ=18)
                                                                                                                                                                                     tggtggCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATc           <  1:534671/236‑1 (MQ=40)
                                                                                                                                                                                       gtggCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTg                    <  1:68634/225‑1 (MQ=39)
                                                                                                                                                                                          gCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACt                                                              <  1:157584/180‑1 (MQ=14)
                                                                                                                                                                                          gCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACt                                                              <  1:537141/180‑1 (MQ=17)
                                                                                                                                                                                          gCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACt                                                              >  2:537139/1‑180 (MQ=14)
                                                                                                                                                                                          gCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACt                                                              >  2:157582/1‑180 (MQ=14)
                                                                                                                                                                                             aaGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCg                                                                                               <  2:709797/144‑1 (MQ=34)
                                                                                                                                                                                                cTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGtt  >  2:446773/1‑234 (MQ=40)
                                                                                                                                                                                                  gACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCAc                                                                                                                                                                                                   <  2:542768/39‑1 (MQ=40)
                                                                                                                                                                                                  gACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGCACGGCTACCGTGACTGCCCCGGa                                                      >  2:232461/1‑180 (MQ=25)
                                                                                                                                                                                                  gACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCgtgaag                                                                                                                    >  2:657683/1‑113 (MQ=255)
                                                                                                                                                                                                      gcgACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCAc                                                                                                      >  1:245607/1‑128 (MQ=35)
                                                                                                                                                                                                                cGGTAACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAAcc                                                                                                                                                    >  2:501016/1‑72 (MQ=21)
                                                                                                                                                                                                                    |                                                                                                                                                                                                                     
ATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTT  >  AE015451/206099‑206524

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: