Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 197,866 G→C 22.8% A1124A (GCG→GCC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451197,8660GC22.8% 121.4 / 31.9 61A1124A (GCG→GCCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (25/22);  new base C (7/7);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01

GTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCA  >  AE015451/197639‑198078
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                    
gTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAAc                                                                                                                                                                                                              >  1:498322/1‑236 (MQ=255)
                  acaGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAAccc                                                                                                                                                                                                            >  2:865525/1‑220 (MQ=255)
                                                    cGACCCCGGCTCAGACCACGATCACCGACTCGGTCGATGCCACCACCGCGACCCTGACGGCCACGCCGAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCgg                                                                                                                                                                                                                     <  1:1130393/177‑1 (MQ=11)
                                                                aGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTacaa                                                                                                                                                                 >  1:291457/1‑217 (MQ=255)
                                                                         tCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGcc                                                                                                                                                                               >  1:497138/1‑194 (MQ=255)
                                                                         tCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTAACGTTGAAACGGGTAAAACCCAGGGAAGCGTCGATTTCCAGAct                                                                                                                                                                                     >  2:437779/1‑187 (MQ=255)
                                                                               aCTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGAACTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCgg                                                                                                                                                                                                  >  1:9063/1‑168 (MQ=255)
                                                                                      cGATACCACTACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACACTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCAAACGGCCAAACCATTACCGTTGAAGCgg                                                                                                                                                                                                                     <  2:532337/143‑1 (MQ=255)
                                                                                                  cgccaCCCTGCCCGCCACGCCGAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATCACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGa                                                                                                                                                                                       <  2:196582/157‑1 (MQ=11)
                                                                                                     taCCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCa                                                                                                                                                                              <  2:291456/166‑1 (MQ=255)
                                                                                                                    gCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGa                                                                                             >  2:677231/1‑233 (MQ=255)
                                                                                                                       gaGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAg                                                                                                                                                                                        <  2:146515/138‑1 (MQ=255)
                                                                                                                            tcaCCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGAc                                                                                                                                                                         <  2:798186/147‑1 (MQ=37)
                                                                                                                               ccGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAg                                                                                 <  2:498321/234‑1 (MQ=255)
                                                                                                                                            gATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTc                                                                                                                                                        >  2:409525/1‑150 (MQ=40)
                                                                                                                                              tCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAAt                                                                <  1:677232/236‑1 (MQ=255)
                                                                                                                                                      cccGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGtcatca                                                                  <  1:93676/225‑1 (MQ=255)
                                                                                                                                                          tgatccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAAc                                                                                                                >  2:938757/5‑176 (MQ=255)
                                                                                                                                                           gatccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAg                                                                                                                                               >  1:613734/4‑144 (MQ=21)
                                                                                                                                                            aCCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATAc                                                  <  1:439501/236‑1 (MQ=255)
                                                                                                                                                             cccTGAGCAATCCTGCCCAGACGCCGGTAACGGTGACCTGGTCCAAGGGCCAAACCATCACCGTAGAAACgga                                                                                                                                                                                                                    >  2:559376/1‑72 (MQ=11)
                                                                                                                                                             cccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCg                                                                                                                                                                                                   >  1:730909/1‑90 (MQ=255)
                                                                                                                                                             cccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCg                                                                                                                                                                                                   <  2:730908/90‑1 (MQ=255)
                                                                                                                                                             cccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTg                                                                        >  2:439500/1‑213 (MQ=255)
                                                                                                                                                             cccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATAcc                                                 >  1:325283/1‑236 (MQ=255)
                                                                                                                                                             cccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATAc                                                  >  1:365255/1‑235 (MQ=255)
                                                                                                                                                             cccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATAc                                                  <  1:865526/235‑1 (MQ=255)
                                                                                                                                                             cccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCg                                                      >  2:750173/1‑231 (MQ=14)
                                                                                                                                                              ccTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATAc                                                  <  2:325282/234‑1 (MQ=255)
                                                                                                                                                               cTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGCCCCCAGCCTCGACTGTCATCaa                                                                 >  1:140925/1‑218 (MQ=255)
                                                                                                                                                                     aaTCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGAc                                                                                                                                                     >  2:311255/1‑128 (MQ=40)
                                                                                                                                                                      aTCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTac                                                                                                                                                                   <  1:152809/113‑1 (MQ=17)
                                                                                                                                                                         aggCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTg                                                                                                    <  2:690764/171‑1 (MQ=25)
                                                                                                                                                                          tGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTaaa                                                                                                                                                                                                                <  2:836230/64‑1 (MQ=255)
                                                                                                                                                                          tGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGGTGAAGCGGGTAAAACCCCGGGCAGCGTCGATTCCCAGACCCCGGCCAATGACGTCTacaaccc                                                                                                                                                              >  2:316969/1‑112 (MQ=255)
                                                                                                                                                                                gACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCa                              <  1:318188/236‑1 (MQ=255)
                                                                                                                                                                                    ccGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGa                                                                           >  2:1056705/1‑187 (MQ=255)
                                                                                                                                                                                    ccGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTAGAACAACGGCTCGACCGTCAGCGTCACCATCGAAAACGGCACGGGCGGCAACTTCGAACATCTGACTCCCAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTa                                    >  1:634175/1‑226 (MQ=18)
                                                                                                                                                                                         aaCCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCAcc                                                                                                                                        >  2:478328/1‑121 (MQ=40)
                                                                                                                                                                                           ccGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACggg                                                                                                                       <  2:1020251/136‑1 (MQ=21)
                                                                                                                                                                                            cGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCgg                                                                                                                                                                                 <  2:803751/77‑1 (MQ=0)
                                                                                                                                                                                                  cctGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGAc                                                                          >  2:598886/3‑174 (MQ=11)
                                                                                                                                                                                                    tGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTcc                                                                                               >  1:243586/1‑151 (MQ=38)
                                                                                                                                                                                                        cAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCg                                                                                                                                                                                  >  2:829552/1‑64 (MQ=35)
                                                                                                                                                                                                        cAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGAcc                                                                                                                                                    >  2:328712/1‑94 (MQ=33)
                                                                                                                                                                                                        cAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGAcc      <  1:614703/236‑1 (MQ=255)
                                                                                                                                                                                                        cAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACTGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGa                                                                           >  2:973749/1‑167 (MQ=255)
                                                                                                                                                                                                         aaCGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTc                                                                                                             <  2:969105/132‑1 (MQ=255)
                                                                                                                                                                                                          aCGGCCAAACCATCACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACcag    <  1:1085523/236‑1 (MQ=255)
                                                                                                                                                                                                            ggCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACcagca  <  1:10381/236‑1 (MQ=255)
                                                                                                                                                                                                              ccAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGa           >  2:880310/1‑225 (MQ=255)
                                                                                                                                                                                                              ccAAACCATTACCGTTGAAGCCGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGaa          <  2:865097/226‑1 (MQ=255)
                                                                                                                                                                                                              ccAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCg                                                                                        <  2:971330/148‑1 (MQ=255)
                                                                                                                                                                                                              ccAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCg                                                                                        >  1:971331/1‑148 (MQ=255)
                                                                                                                                                                                                                  aCCATTACCGGTGAAGCGGGGAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAACGACGTCTACAACAACg                                                                                                                                                            >  2:385867/1‑76 (MQ=11)
                                                                                                                                                                                                                   ccATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCg                                                                                              <  2:757229/137‑1 (MQ=255)
                                                                                                                                                                                                                   ccATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCg                                                                                              >  1:757230/1‑137 (MQ=39)
                                                                                                                                                                                                                     atcactgTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATAc                                                  <  1:598887/173‑1 (MQ=255)
                                                                                                                                                                                                                     atcacggTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATa                                                   <  2:617476/172‑1 (MQ=255)
                                                                                                                                                                                                                     atcaCCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAgcaatgg                       <  1:750174/203‑1 (MQ=255)
                                                                                                                                                                                                                          cGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTAt                                   >  1:339853/1‑189 (MQ=255)
                                                                                                                                                                                                                           gTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCTCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGTCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATTCCGTCACCGTAAGTATTGt                                <  2:1067660/191‑1 (MQ=255)
                                                                                                                                                                                                                               aaGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGATCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAACCTCGACTGTCATCAATg                                                               >  2:700668/1‑156 (MQ=255)
                                                                                                                                                                                                                               aaGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGa                                                                           <  2:727783/144‑1 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                    
GTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCA  >  AE015451/197639‑198078

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