Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 196,027 T→A 9.5% T511T (ACT→ACA PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451196,0270TA 9.5% 182.1 / 5.8 68T511T (ACT→ACAPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (32/26);  new base A (2/4);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACA  >  AE015451/195802‑196229
                                                                                                                                                                                                                                 |                                                                                                                                                                                                          
ggCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCGCCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGt                                                                                                                                                                                                  <  1:435371/236‑1 (MQ=255)
                        cTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACAGTGAccc                                                                                                                                                                                                     >  1:42948/1‑209 (MQ=25)
                           gACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGt                                                                                                                                                                                                  <  1:1243894/209‑1 (MQ=2)
                            aCCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAAc                                                                                                                                                                                             <  2:658389/213‑1 (MQ=14)
                              cGTCAGCGTCCCCATCGAGAGCGCCACTGGCGGCAACTGCGAACAGCTGCCCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                    <  1:139631/236‑1 (MQ=255)
                                 cAGCGTCACCATCGAGAACGCCACAGGCGGCAACTTCGAACAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                    <  2:1080852/233‑1 (MQ=255)
                                                 agCGCCACTGGCGACAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCTTGACCCTGTctc                                                                                                                                                                                               >  1:791932/1‑188 (MQ=255)
                                                    gcCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAAc                                                                                                                                                                                             <  2:35461/189‑1 (MQ=25)
                                                      cACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACTGTGACCCTGTc                                                                                                                                                                                                 >  1:399208/1‑183 (MQ=18)
                                                      cACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGAtt                                                                                                                                            >  1:602176/1‑236 (MQ=255)
                                                         tGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAAc                                                                                                                                                                                     <  2:860071/192‑1 (MQ=38)
                                                           gcggcAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGGAAAACCCAGGGCAGCGTCGATTTGCAg                                                                                                                                       >  1:370754/1‑236 (MQ=255)
                                                              gcAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGcc                                                                                                                                                                                         >  1:224762/1‑183 (MQ=18)
                                                                             tGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGtt                                                                                                                                                                           <  1:103304/182‑1 (MQ=37)
                                                                               aCCCCGAACGCGACCCCGGCTCAGACCACGATCCACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                    <  2:314546/187‑1 (MQ=16)
                                                                                 cccGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAAc                                                                                                                                                                                     <  2:205013/168‑1 (MQ=18)
                                                                                              ccGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGAc                                                                                                    >  1:583933/1‑236 (MQ=37)
                                                                                                         cACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACtgtg                                                                                                                                                                                                         >  2:242637/1‑124 (MQ=17)
                                                                                                         cACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACtgtg                                                                                                                                                                                                         <  1:242641/124‑1 (MQ=38)
                                                                                                               cAACGACTCGGTCGATACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACTGTGACCCTGTCCAACGGCCAAACgc                                                                                                                                                                                   >  2:1200975/1‑139 (MQ=2)
                                                                                                               cAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGGGACCCTGTCCAACGGCCAAACGATTCCCGTTGa                                                                                                                                                                         >  2:390943/1‑150 (MQ=255)
                                                                                                                    aCTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGaa                                                                               >  1:102651/1‑235 (MQ=39)
                                                                                                                         gTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGcc                                                                                                                                                                                         >  1:450404/1‑124 (MQ=255)
                                                                                                                               cccaccaccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGa                                                                                                                                                                         <  2:357628/133‑1 (MQ=37)
                                                                                                                                 caccaccGCGACCCTGCCGGCGAGCCCGTCGGTCACAGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCgag                                                                                 <  1:71278/220‑1 (MQ=255)
                                                                                                                                       cgcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGc                                                                                             >  1:386603/1‑202 (MQ=40)
                                                                                                                                         cgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTc                                                                                                                                                                                                 <  1:705328/100‑1 (MQ=255)
                                                                                                                                            cccTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAg                                                                                                                                       <  2:602171/155‑1 (MQ=38)
                                                                                                                                                  cGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGcc                                                                           <  2:583928/209‑1 (MQ=255)
                                                                                                                                                        gCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGAACCTGTTCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGccgg                                                                                                                                                   >  2:106935/1‑128 (MQ=11)
                                                                                                                                                        gCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGaa                                           <  1:1046099/235‑1 (MQ=255)
                                                                                                                                                           cgaccgTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGa                                                                                                                                                                         <  2:1101418/101‑1 (MQ=255)
                                                                                                                                                                               gggATCACGTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                    <  1:1150746/89‑1 (MQ=255)
                                                                                                                                                                                tcagcACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGa                                      <  1:383873/212‑1 (MQ=255)
                                                                                                                                                                                    cACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAAc               <  1:1166170/235‑1 (MQ=255)
                                                                                                                                                                                            gcGTAACCCAGACAAATCCTGCCCAGACGCCGGTAACAGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                    <  2:974351/77‑1 (MQ=11)
                                                                                                                                                                                               tGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTc                                                                                                                                                                                                 <  1:1030783/46‑1 (MQ=38)
                                                                                                                                                                                               tGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTc                                                                                                                                                                                                 >  2:1030778/1‑46 (MQ=40)
                                                                                                                                                                                                gACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGa                                                                                                                                                                                                        >  1:1176785/1‑38 (MQ=40)
                                                                                                                                                                                                gACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGc                                                                                                                                                                      <  2:1258569/72‑1 (MQ=255)
                                                                                                                                                                                                gACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGaca  >  1:133866/1‑236 (MQ=255)
                                                                                                                                                                                                gACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGccg                                                                                                                                                                    >  1:488114/1‑72 (MQ=12)
                                                                                                                                                                                                gACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGCAAAACCCAGGGc                                                                                                                                                      >  1:790557/1‑88 (MQ=255)
                                                                                                                                                                                                 aacaagAGCCACCCTGCCCAGACGCCGGTAACAGTGACCCTGTCCAACGGCAAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGAt                   <  1:293414/213‑1 (MQ=255)
                                                                                                                                                                                                  cccTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGt                                                                                                                                                                                                  <  2:1036734/42‑1 (MQ=40)
                                                                                                                                                                                                  cccTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGa                                                                                                                                                                         >  2:906636/1‑67 (MQ=255)
                                                                                                                                                                                                  cccTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGa                                                                                                                                                                         <  1:906641/67‑1 (MQ=255)
                                                                                                                                                                                                  cccTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGAc                         >  2:923779/1‑211 (MQ=255)
                                                                                                                                                                                                  cccTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCCCCGACCACGATCCACGACTCGGTCGa    >  2:1166165/1‑232 (MQ=255)
                                                                                                                                                                                                  cccTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGa                                                                                                     >  2:786271/1‑135 (MQ=21)
                                                                                                                                                                                                    cTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGc                                                                                                                                                                      >  1:687066/1‑68 (MQ=255)
                                                                                                                                                                                                       aGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCCTctctc                                                                                      >  1:637333/1‑141 (MQ=255)
                                                                                                                                                                                                        gCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACgg                                                                                                          >  2:286757/1‑124 (MQ=38)
                                                                                                                                                                                                        gCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGt       >  1:1122647/1‑223 (MQ=19)
                                                                                                                                                                                                          aaTCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTAcc                                                                                                                                                                              <  2:288415/54‑1 (MQ=255)
                                                                                                                                                                                                          aaTCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCACTGACGTCTACAACAACGGTTCTACGGTCATCGTCACCATCGAGAACGCCACggg                                                                      >  2:85823/1‑158 (MQ=255)
                                                                                                                                                                                                          aaTCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAAtt                                                              >  2:961655/1‑166 (MQ=255)
                                                                                                                                                                                                          aaTCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTc                             >  2:527513/1‑199 (MQ=255)
                                                                                                                                                                                                          aaTCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGa    >  2:907254/1‑223 (MQ=255)
                                                                                                                                                                                                               tGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACg              <  2:1086886/209‑1 (MQ=255)
                                                                                                                                                                                                               tGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGt                                                                                                                                                                  <  1:238129/61‑1 (MQ=40)
                                                                                                                                                                                                               tGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGt                                                                                                                                                                  >  2:238125/1‑61 (MQ=40)
                                                                                                                                                                                                                      aCGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGCTTACCGTTGAAGCCGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACg              <  1:19955/202‑1 (MQ=255)
                                                                                                                                                                                                                        gCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGa                          <  2:326263/188‑1 (MQ=255)
                                                                                                                                                                                                                         ccGGTAACTGTGACCCTGTCCAACGGCCAAAc                                                                                                                                                                                     >  2:828572/1‑32 (MQ=40)
                                                                                                                                                                                                                         ccGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                    >  1:364162/1‑49 (MQ=255)
                                                                                                                                                                                                                         ccGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCg                                                           <  2:946917/154‑1 (MQ=255)
                                                                                                                                                                                                                         ccGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCg                                                           >  1:946922/1‑154 (MQ=255)
                                                                                                                                                                                                                         ccGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCAt                                                                                     >  1:652019/1‑128 (MQ=25)
                                                                                                                                                                                                                         ccGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGCAAAACCCAGGGCatc                                                                                                                                                   >  1:325802/1‑64 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                          
GGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACA  >  AE015451/195802‑196229

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