Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 213,652 C→T 8.4% G6386G (GGC→GGT PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451213,6520CT 8.4% 177.7 / 3.0 61G6386G (GGC→GGTPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (27/27);  new base T (2/3);  total (30/31)
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.88e-01

AGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGC  >  AE015451/213433‑213871
                                                                                                                                                                                                                           |                                                                                                                                                                                                                           
aGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCaacaac                                                                                                                                                                                                              <  1:260516/235‑1 (MQ=255)
          gTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGc                                                                                                                                                                                                      <  2:403455/233‑1 (MQ=255)
                       cTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCGTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGATGGTCTGCCTATcgacaagcatgcggc                                                                                                                                                                                                      >  1:660931/1‑211 (MQ=11)
                       cTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCGTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGATGGTCTGCCTATcgacaagcatgcggc                                                                                                                                                                                                      <  2:660930/220‑10 (MQ=11)
                            cAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGa                                                                                                                                                                                 >  1:320422/1‑236 (MQ=255)
                                  aGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACCGTGCATTGACCTTCACGCt                                                                                                                                                                                       >  2:681507/1‑224 (MQ=255)
                                                 tCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGTCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAAc                                                                                                                                                                        <  2:462922/224‑1 (MQ=255)
                                                 tCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAAc                                                                                                                                                                        >  1:462923/1‑224 (MQ=255)
                                                      aGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGAccc                                                                                                                                                                                              >  2:135904/1‑196 (MQ=255)
                                                      aGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACAGCCTGCCGATCAACAACCACAGTGCATTGACCTTCa                                                                                                                                                                                           >  2:669106/1‑200 (MQ=255)
                                                           aCGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCt                                                                                                                                                                                                                           <  2:22415/163‑1 (MQ=255)
                                                           aCGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGt                                                                                                                                                                     >  1:178318/1‑217 (MQ=40)
                                                           aCGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAAc                                                                                                                                                  <  1:632367/236‑1 (MQ=255)
                                                                aCTGGTGCGATCTTCGAGAATCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATcaacaa                                                                                                                                                                                                               <  2:692615/170‑1 (MQ=255)
                                                                            ttCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGCCGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGCAGGCGGCGAGATTACCTACACCATCACGCTGACCAACACAGACGGCCTGCCGACCAACAACcacact                                                                                                                                                                                                        >  2:11230/1‑163 (MQ=255)
                                                                            ttCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGc                                                                                                                                                                                                                         >  2:565595/1‑148 (MQ=255)
                                                                            tgcGAGAACCTGCAACTGGGCGGCAACGCTTGGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGAGAGATGACCTACACCATCACGCTGACCAACAAAGACGGGCTGCCGATCAACAGCCACAGTGCATTGACCTTCACGCTGAGCGACGGc                                                                                                                                                                             <  2:9628/190‑1 (MQ=255)
                                                                               gagaACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCAcc                                                                                                                                                              <  1:686922/204‑1 (MQ=38)
                                                                                   aCCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCg                                                                                                                            <  1:687930/234‑1 (MQ=255)
                                                                                                 ggcAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACa                                                                                                             >  1:436108/1‑235 (MQ=255)
                                                                                                  gcAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAAcc                                                                                                                                                                                                             >  2:720414/1‑138 (MQ=255)
                                                                                                  gcAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAAcc                                                                                                                                                                                                             <  1:720415/138‑1 (MQ=255)
                                                                                                      cGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCaa                                                                                                                                                                           >  1:346329/1‑168 (MQ=255)
                                                                                                                    tGACCGGCACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCg                                                                                                                            <  1:405605/201‑1 (MQ=38)
                                                                                                                    ggACCGACACCGCCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCGCAGTGCATTGACCTTCACGCTGAGCGACGGc                                                                                                                                                                             <  2:182640/151‑1 (MQ=255)
                                                                                                                       ccGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCg                                                                                                                            <  1:731757/198‑1 (MQ=37)
                                                                                                                             ccaccGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCg                                                                                                                            >  2:661395/1‑192 (MQ=37)
                                                                                                                             ccaccGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCg                                                                                                                            <  1:661396/192‑1 (MQ=37)
                                                                                                                             ccaccGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAAc                                                                                                           <  2:54997/209‑1 (MQ=40)
                                                                                                                             ccaccGACGAAGTAGTGGCCAAGGTGACCGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACGGTCTGCCGATTAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACg                                                                                        <  1:697419/228‑1 (MQ=11)
                                                                                                                                                    ttGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGa                                                         >  1:64844/1‑236 (MQ=40)
                                                                                                                                                      gACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCgg                                                                                                                >  2:538124/1‑179 (MQ=40)
                                                                                                                                                      gACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCgg                                                                                                                <  1:538125/179‑1 (MQ=255)
                                                                                                                                                      gACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGACTGCCGATcgacaagc                                                                                                                                                                                                             >  2:167238/1‑84 (MQ=255)
                                                                                                                                                              ccccGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCACCGGCACCGTGGGTACTGCCCCCGTGACTGCTCCGGACac                                                                                                            >  2:32975/1‑174 (MQ=255)
                                                                                                                                                                ccGTCGGTGACCGAAGGCGGCGAGATTACCTACGCCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCAGTGACCTGCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTc                                                                                                  <  2:15435/183‑1 (MQ=255)
                                                                                                                                                                ccGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGc                                                                                                                                                                                                                         <  2:712854/64‑1 (MQ=35)
                                                                                                                                                                ccGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGc                                                                                                                                                                                                                         >  1:712855/1‑64 (MQ=35)
                                                                                                                                                                ccGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATc                                                                                                                                                                                                                    <  2:57382/69‑1 (MQ=35)
                                                                                                                                                                  ggcGGTGACCGAAGGCGGCGAGATTACCTACACCAGCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACg                                                                                        <  1:3073/189‑1 (MQ=255)
                                                                                                                                                                  gTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCACCAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGt                                                                                               <  2:127673/184‑1 (MQ=39)
                                                                                                                                                                    cGGTGACCGAAGGCGGCGAGATTACCGACACCATCCCGCTGACCAGCAAAGACGGCCTGCCGGTCACCAACCACAGTGCATTGACTGTCACGCTGAGCGACgg                                                                                                                                                                              >  2:346328/1‑103 (MQ=11)
                                                                                                                                                                                   gcgAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTg                          <  1:746550/236‑1 (MQ=40)
                                                                                                                                                                                      agaTTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCATGGGTACTGCCACCGTGACTGCCCCGGAcct                                                                                                            >  2:340104/1‑149 (MQ=9)
                                                                                                                                                                                      agaTTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGGGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCg                                  >  1:191095/1‑225 (MQ=255)
                                                                                                                                                                                                    aTCACGCTGACCAACATAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAACGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGc                                                                                                                                 <  2:705574/116‑1 (MQ=16)
                                                                                                                                                                                                    aTCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCTACCGtt                                                           >  2:409376/1‑186 (MQ=37)
                                                                                                                                                                                                    aTCACGCTGACAAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGc                                                                                                                                 >  1:705575/1‑116 (MQ=17)
                                                                                                                                                                                                        cccTGACCAACAAAGACGGTCTGCCGATTAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCg                                                                                                                 <  2:408023/126‑1 (MQ=255)
                                                                                                                                                                                                        cccTGACCAACAAAGACGGTCTGCCGATTAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCg                                                                                                                 >  1:408024/3‑128 (MQ=255)
                                                                                                                                                                                                        cGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGt                                                                                                   >  1:547003/1‑142 (MQ=38)
                                                                                                                                                                                                        cGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGt                                                                                                   <  2:547002/142‑1 (MQ=39)
                                                                                                                                                                                                        cGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGATGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTc     >  1:693641/1‑236 (MQ=255)
                                                                                                                                                                                                           tGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGc  <  1:600243/236‑1 (MQ=255)
                                                                                                                                                                                                              ccAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTa                                                                                                      >  2:142811/1‑133 (MQ=37)
                                                                                                                                                                                                                aacaaAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCg                                                                                                                 <  1:665315/120‑1 (MQ=40)
                                                                                                                                                                                                                aacaaAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCg                                                                                                                 >  2:665314/1‑120 (MQ=40)
                                                                                                                                                                                                                aacaaAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTg                                                                             >  2:410326/1‑156 (MQ=37)
                                                                                                                                                                                                                    aaGACGGCCGGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTg                          <  1:569509/203‑1 (MQ=255)
                                                                                                                                                                                                                          gCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTg                                                                             >  1:364551/1‑146 (MQ=39)
                                                                                                                                                                                                                          gCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTg                                                                             <  2:364550/146‑1 (MQ=39)
                                                                                                                                                                                                                           |                                                                                                                                                                                                                           
AGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGC  >  AE015451/213433‑213871

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