Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 197,834 T→C 20.6% L1114L (TTG→CTG)  PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451197,8340TC20.6% 136.3 / 20.7 61L1114L (TTG→CTG) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (23/24);  new base C (7/6);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCA  >  AE015451/197611‑198050
                                                                                                                                                                                                                               |                                                                                                                                                                                                                        
agacTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAACTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACg                                                                                                                                                                                                                 <  1:684160/230‑1 (MQ=255)
                            gTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAAc                                                                                                                                                                                  >  1:498322/1‑236 (MQ=255)
                                              acaGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCACACCAtt                                                                                                                                                                                                      >  2:839790/1‑198 (MQ=255)
                                              acaGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAAccc                                                                                                                                                                                >  2:865525/1‑220 (MQ=255)
                                                               gCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGACGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGCCCAACGGCCAAACCACTAcc                                                                                                                                                                                                   >  2:30803/1‑184 (MQ=255)
                                                                        cccGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACGCCACCACCGCGACCCTGACCGCGTCCCCGAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGa                                                                                                                                                                                              <  2:613423/180‑1 (MQ=11)
                                                                              tccGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACTGCGACCCTGACCGCGTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAAc                                                                                                                                                                                                                  <  2:1016273/153‑1 (MQ=255)
                                                                                cGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTg                                                                                                                                                                                               >  1:769370/1‑171 (MQ=255)
                                                                                cGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTg                                                                                                                                                                                               <  2:769369/171‑1 (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)
                                                                                                      cACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGGGACCCTGTCCAACGGCCAAACCa                                                                                                                                                                                                        >  2:342751/1‑140 (MQ=255)
                                                                                                           aCTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGAACTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCgg                                                                                                                                                                      >  1:9063/1‑168 (MQ=255)
                                                                                                                  cGATACCACTACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACACTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCAAACGGCCAAACCATTACCGTTGAAGCgg                                                                                                                                                                                         <  2:532337/143‑1 (MQ=255)
                                                                                                                  cGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCGTGTCCAACGGCCAAACCATTACCGtgg                                                                                                                                                                                               >  2:162546/1‑135 (MQ=255)
                                                                                                                   gACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACACTGACCCTGTc                                                                                                                                                                                                                      >  1:52742/1‑113 (MQ=14)
                                                                                                                              cgccaCCCTGCCCGCCACGCCGAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATCACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGa                                                                                                                                                           <  2:196582/157‑1 (MQ=11)
                                                                                                                                 taCCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCa                                                                                                                                                  <  2:291456/166‑1 (MQ=255)
                                                                                                                                                gCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGa                                                                 >  2:677231/1‑233 (MQ=255)
                                                                                                                                                  cgtccGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTc                                                                                                                                                                                                                      <  2:482503/78‑1 (MQ=255)
                                                                                                                                                  cgaccGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGa                                                                                                                                                                                              <  2:1016795/102‑1 (MQ=255)
                                                                                                                                                   gaGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAg                                                                                                                                                            <  2:146515/138‑1 (MQ=255)
                                                                                                                                                       gtcaCCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTg                                                                                                                                                                                               <  2:547145/97‑1 (MQ=18)
                                                                                                                                                        tcaCCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGAc                                                                                                                                             <  2:798186/147‑1 (MQ=37)
                                                                                                                                                           ccGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAg                                                     <  2:498321/234‑1 (MQ=255)
                                                                                                                                                                        gATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTc                                                                                                                            >  2:409525/1‑150 (MQ=40)
                                                                                                                                                                         aTCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGt                                                                                                                                                                                                                       >  1:173850/1‑58 (MQ=35)
                                                                                                                                                                         aTCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGt                                                                                                                                                                                                                       <  2:173850/58‑1 (MQ=35)
                                                                                                                                                                          tCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAAt                                    <  1:677232/236‑1 (MQ=255)
                                                                                                                                                                                  cccGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGtcatca                                      <  1:93676/225‑1 (MQ=255)
                                                                                                                                                                                    cGTGACCGAGAGAGATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAAc                                                                                                                                                                                                          <  2:426837/60‑1 (MQ=14)
                                                                                                                                                                                      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)
                                                                                                                                                                                         cccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAAc                                                                                                                                                                                                                  <  2:924505/47‑1 (MQ=38)
                                                                                                                                                                                         cccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGccg                                                                                                                                                                                         >  1:215979/1‑70 (MQ=17)
                                                                                                                                                                                         cccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCg                                                                                                                                                                       <  2:730908/90‑1 (MQ=255)
                                                                                                                                                                                         cccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCg                                                                                                                                                                       >  1:730909/1‑90 (MQ=255)
                                                                                                                                                                                         cccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTg                                            >  2:439500/1‑213 (MQ=255)
                                                                                                                                                                                         cccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATAcc                     >  1:325283/1‑236 (MQ=255)
                                                                                                                                                                                         cccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATAc                      <  1:865526/235‑1 (MQ=255)
                                                                                                                                                                                         cccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATAc                      >  1:365255/1‑235 (MQ=255)
                                                                                                                                                                                         cccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCg                          >  2:750173/1‑231 (MQ=14)
                                                                                                                                                                                          ccTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATAc                      <  2:325282/234‑1 (MQ=255)
                                                                                                                                                                                           cTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGCCCCCAGCCTCGACTGTCATCaa                                     >  1:140925/1‑218 (MQ=255)
                                                                                                                                                                                              aGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAAc                                                                                                                                                                                                                  <  1:1069790/42‑1 (MQ=38)
                                                                                                                                                                                                 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)
                                                                                                                                                                                                           aGACGCCGGTAACCGTGACCTCGTCCAACGCCCACACCact                                                                                                                                                                                                      >  2:81762/1‑39 (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)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                        
CGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCA  >  AE015451/197611‑198050

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