Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 206,332 C→T 9.2% T3946T (ACC→ACT PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451206,3320CT 9.2% 185.8 / 6.7 63T3946T (ACC→ACTPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (29/28);  new base T (2/4);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTT  >  AE015451/206131‑206524
                                                                                                                                                                                                         |                                                                                                                                                                                                
ccACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGc                                                                                                                                                                <  1:433553/236‑1 (MQ=255)
      tGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGc                                                                                                                                                                >  1:130420/1‑230 (MQ=255)
        cgcAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGt                                                                                                                                                                    >  2:313943/1‑224 (MQ=255)
         gcAAGGCGATGACGTCTACACCGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGc                                                                                                                                                                                     <  1:163635/206‑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)
                                                   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)
                                                                   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                                                                                                                                             <  2:758152/176‑1 (MQ=38)
                                                                               ggTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGc                                                                                                                                             >  1:758154/1‑176 (MQ=39)
                                                                                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)
                                                                                                             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)
                                                                                                                                                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:537141/180‑1 (MQ=17)
                                                                                                                                                          gCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACt                                                              >  2:157582/1‑180 (MQ=14)
                                                                                                                                                          gCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACt                                                              >  2:537139/1‑180 (MQ=14)
                                                                                                                                                          gCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACt                                                              <  1:157584/180‑1 (MQ=14)
                                                                                                                                                             aaGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCg                                                                                               <  2:709797/144‑1 (MQ=34)
                                                                                                                                                                cTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGtt  >  2:446773/1‑234 (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)
                                                                                                                                                                                                         |                                                                                                                                                                                                
CCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTCGCACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTT  >  AE015451/206131‑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: