Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 197,926 C→T 10.5% G1144G (GGC→GGT PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451197,9260CT10.5% 195.3 / 8.1 74G1144G (GGC→GGTPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (30/36);  new base T (4/4);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

GCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGCCAGGCGCTGGACCGTCCGTTGACGGTAACCCTGTCTAACGGCGACACCGTGACGATC  >  AE015451/197755‑198160
                                                                                                                                                                           |                                                                                                                                                                                                                                          
gCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGa                                                                                                                                                                               >  2:677231/1‑233 (MQ=255)
           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)
                                         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)
                                                     aggCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTg                                                                                                                                                                                      <  2:690764/171‑1 (MQ=25)
                                                            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)
                                                                              cctGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGAc                                                                                                                                                            >  2:598886/3‑174 (MQ=11)
                                                                                tGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTcc                                                                                                                                                                                 >  1:243586/1‑151 (MQ=38)
                                                                                    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                                                                                                                                                                          >  1:971331/1‑148 (MQ=255)
                                                                                          ccAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCg                                                                                                                                                                          <  2:971330/148‑1 (MQ=255)
                                                                                               ccATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCg                                                                                                                                                                                >  1:757230/1‑137 (MQ=39)
                                                                                               ccATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCg                                                                                                                                                                                <  2:757229/137‑1 (MQ=255)
                                                                                                 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)
                                                                                                                 gcgaaaCCCAGGGCAGCGTCGATTTCCACACCCCGGCCAATGACGTCGACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACAGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACcag                                                                                      <  1:197904/206‑1 (MQ=255)
                                                                                                                   aaaaCCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAGCGCCACGggc                                                                                                                                                                                                        >  1:533384/1‑93 (MQ=18)
                                                                                                                   aaaaCCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAGCGCCACGggc                                                                                                                                                                                                        <  2:533383/93‑1 (MQ=21)
                                                                                                                     caCCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTaa                                                             <  1:313191/229‑1 (MQ=255)
                                                                                                                      aCCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTAc                                                                        <  1:938758/218‑1 (MQ=255)
                                                                                                                            ggCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACGGGCGGCAACTTCGAACATCTGACCCCGAATCCGACCCCAGCCTCGa                                                                                                                                                             >  2:585045/1‑127 (MQ=11)
                                                                                                                              cAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTaa                                                             <  2:997691/221‑1 (MQ=255)
                                                                                                                              cAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTaa                                                             >  1:997692/1‑221 (MQ=255)
                                                                                                                                  gTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGCCAGGCGCTGGAccgtcc                                          <  1:888240/236‑1 (MQ=255)
                                                                                                                                  gTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGCCAGGCGCTGGAccgtcc                                          <  1:973750/236‑1 (MQ=255)
                                                                                                                                   tCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGCACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCAcc                                                                                                                             >  2:108348/1‑152 (MQ=255)
                                                                                                                                    cGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGa                                                                   >  1:199384/1‑209 (MQ=255)
                                                                                                                                     gATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAg                                                            >  2:888239/1‑215 (MQ=255)
                                                                                                                                        gtCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAAt                                                                                                                                                  <  2:154182/125‑1 (MQ=255)
                                                                                                                                         tCCAGACCCCGGCGAATGACGTCTACAACAACGGCTCAACTGTCAGCGTCACCATCGAGAACGCCACGGGCGGCAACTTCGAGCATCTGAccc                                                                                                                                                                                  >  2:631832/1‑91 (MQ=11)
                                                                                                                                         tCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACtgtg                                                                    >  2:586076/1‑203 (MQ=255)
                                                                                                                                         tCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGCCAGGCGCt                                                   >  2:544607/1‑220 (MQ=255)
                                                                                                                                             taCCCCACCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGa                                                                                                                                                                         <  2:155914/97‑1 (MQ=255)
                                                                                                                                               ccccGGCGAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAAt                                                                                                                                                  <  2:314438/119‑1 (MQ=255)
                                                                                                                                                 ccGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGtc                                                                                                                                                        <  1:392580/111‑1 (MQ=255)
                                                                                                                                                       aaTGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCATCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAAtgt                                                                                                   >  2:395093/1‑158 (MQ=255)
                                                                                                                                                            cgacTACAACAACTGTTCGACCGTCAGCGTCACCATCGAGAGCGCCACAGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTaa                                                             <  2:969452/188‑1 (MQ=255)
                                                                                                                                                            cGTCTACAACACCGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTTACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATa                                                                                                                                     <  2:60768/119‑1 (MQ=255)
                                                                                                                                                             gTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTaa                                                             <  2:828167/190‑1 (MQ=255)
                                                                                                                                                                  ccacaacGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCGGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGCCAGGCGCTg                                                  <  2:73886/194‑1 (MQ=255)
                                                                                                                                                                    acagcGGTTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCAc                                                                                                                              <  2:386400/118‑1 (MQ=255)
                                                                                                                                                                    acaacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGCAATTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGAcca                                                                                                                                                          >  1:647328/1‑88 (MQ=14)
                                                                                                                                                                     gacaggCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACa                                                                                                                                              <  2:291984/97‑1 (MQ=255)
                                                                                                                                                                      aacGGCTCGACCGTCAGCGTCACCATCGCGAGCGCCACGGGCGGCATCTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGt                                                                                                                  <  2:85359/128‑1 (MQ=255)
                                                                                                                                                                      aacGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTaa                                                             <  2:365254/181‑1 (MQ=255)
                                                                                                                                                                        cGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGCCAGGCGCTGGAc                                               >  2:302633/1‑193 (MQ=255)
                                                                                                                                                                           cTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCg                                                                                                                                                                                              >  2:357637/1‑47 (MQ=38)
                                                                                                                                                                           cTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGCCAGGCGCTGGACCGTCCGTTGACGGTAACCCTGTCTAACGGCGACACCGTGACGATc  <  1:839791/235‑1 (MQ=255)
                                                                                                                                                                           |                                                                                                                                                                                                                                          
GCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGCCAGGCGCTGGACCGTCCGTTGACGGTAACCCTGTCTAACGGCGACACCGTGACGATC  >  AE015451/197755‑198160

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