Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 197,632 T→C 19.6% T1046T (ACT→ACC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451197,6320TC19.6% 151.1 / 23.4 78T1046T (ACT→ACCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (38/24);  new base C (8/8);  total (46/32)
Fisher's exact test for biased strand distribution p-value = 5.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

GGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGG  >  AE015451/197422‑197867
                                                                                                                                                                                                                  |                                                                                                                                                                                                                                           
ggTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGGCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACGGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCAt                                                                                                                                                                                                                               >  2:680512/1‑225 (MQ=9)
                                   cGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCg                                                                                                                                                                                                                             >  2:435229/1‑192 (MQ=16)
                                            ccGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCa                                                                                                                                                                                                 <  2:1009190/211‑1 (MQ=17)
                                                                   acCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGCCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGcc                                                                                                                                                                                                                     <  2:1161988/168‑1 (MQ=38)
                                                                   acCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGCCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGcc                                                                                                                                                                                                                     >  1:1161991/1‑168 (MQ=38)
                                                                       tGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGCTCGACTGTCAGCGTCACCATCGAGAACGCCACGGGCGGCAACTTCGaa                                                                                                                                                                                                   >  1:198613/1‑181 (MQ=255)
                                                                        gacccTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGACTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAAcc                                                                                                                                                                                   >  1:181562/4‑197 (MQ=11)
                                                                        gATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGCTCGACTGTCAGCGTCACCATCGAGAACGCCACGGGCGGCAACTTCGAACACCTGACCCCg                                                                                                                                                                                       >  1:853969/1‑193 (MQ=255)
                                                                          tCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGcc                                                                                                                                                                            >  2:849520/1‑202 (MQ=40)
                                                                          tCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCAc                                                                                                                                                                >  2:989573/1‑214 (MQ=39)
                                                                          tCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATCACCGTTGAAGCCGGCAAGACCCAGGCCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGAc                                                                                                                                                                               >  1:715827/1‑199 (MQ=255)
                                                                          cccTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGacac                                                                                                                                           >  1:1118946/2‑235 (MQ=11)
                                                                            cTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCgagagc                                                                                                                                                                                                                        >  2:124793/1‑154 (MQ=2)
                                                                                gCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGCAAATTc                                                                                                                                                                                                      >  2:343004/1‑170 (MQ=11)
                                                                                  aacccTGCCCAGACGCCGGTGACCGTGTCCCTGTCCAACGGCCAAACGAGGACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAAt                                                                                                                                                                                                         >  1:476702/4‑165 (MQ=11)
                                                                                  aaTCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCg                                                                                                                                                                                                                             >  2:87545/1‑145 (MQ=37)
                                                                                  aaTCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCggca                                                                                                                                                                                                           >  2:427973/1‑161 (MQ=25)
                                                                                    tCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCa                                                                                                                                                                                                                                         <  2:660824/131‑1 (MQ=18)
                                                                                       tGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCaac                                                                                                                                                         >  2:976936/1‑206 (MQ=2)
                                                                                       tGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCa                                                                                                                                                                                                 >  2:295418/1‑168 (MQ=9)
                                                                                       tGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGCAATTTCGAACAACTGACCCCGAACCCGACCCCGGCCCCGACCACCATc                                                                                                                                                            >  1:21256/1‑205 (MQ=255)
                                                                                       tGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGCTCGACTGTCAGCGTCACCATCGAGaa                                                                                                                                                                                                                         >  1:115677/1‑144 (MQ=37)
                                                                                                gCCGGTCACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGcc                                                                                                                                                                                                                     <  2:380439/139‑1 (MQ=255)
                                                                                                gCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGa                                                                                                                                              <  2:1118943/210‑1 (MQ=40)
                                                                                                 ccGGTAACCATGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGAATTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCgac                                                                                                                                                                                                                           >  2:937121/1‑131 (MQ=12)
                                                                                                      gacagTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCg                                                                                                                                                                                 <  2:888905/165‑1 (MQ=14)
                                                                                                      gaCCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGAc                                                                                                                                                                                           >  1:1098045/2‑159 (MQ=255)
                                                                                                      gaCCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGAc                                                                                                                                                                                           <  2:1098042/158‑1 (MQ=255)
                                                                                                              cccTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGa                                                                                                                                                                                <  2:1245665/162‑1 (MQ=33)
                                                                                                               ccTGTCCAACGGCCAAACGATTACCGTTGATGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACATGCGGTAATTTCGAGCAACCGACCCAGAACCCGACCCCGGc                                                                                                                                                                         >  2:642253/1‑168 (MQ=11)
                                                                                                               ccTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACACGACGCCGGCACAGa                                                                                                                                                                    >  2:65853/1‑173 (MQ=25)
                                                                                                                cTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGCTCGACTGTCAGCGTCACCATCGAGAGCGCCACAGGCGGGaacgtcgg                                                                                                                                                                                                    >  1:174463/1‑134 (MQ=11)
                                                                                                                     cAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAAAGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTaaa                                                                                                                                                                                                         >  1:233859/1‑129 (MQ=11)
                                                                                                                       aCGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAggcggc                                                                                             <  1:620209/236‑1 (MQ=255)
                                                                                                                       aCGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAggcggc                                                                                             <  1:65853/236‑1 (MQ=255)
                                                                                                                         ggCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGt                                                                                           <  1:989576/236‑1 (MQ=255)
                                                                                                                            cataCCATTTGCGTTGGAGCCCGCAAGACCCAGGGCAGCGTCGATTTCCCGACCCCGGCGCAGGACGTCTACAACAACGGTTCGACTGt                                                                                                                                                                                                                                           <  2:668154/86‑1 (MQ=11)
                                                                                                                               aCCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCCAAGGCGGCGTGATCAc                                                                                     <  1:849523/236‑1 (MQ=255)
                                                                                                                                ccATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACg                                                                                                                                                                                                                       <  1:694150/105‑1 (MQ=32)
                                                                                                                                ccATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACg                                                                                                                                                                                                                       >  2:694147/1‑105 (MQ=32)
                                                                                                                                ccATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAGCaa                                                                                                                                                                                                >  2:874618/1‑128 (MQ=2)
                                                                                                                                ccATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCg                                                                                                                                                   >  2:470680/1‑173 (MQ=38)
                                                                                                                                ccATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCg                                                                                                                                                   <  1:470683/173‑1 (MQ=39)
                                                                                                                                  aTTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCacta                                                                                                                                      <  1:1186743/184‑1 (MQ=255)
                                                                                                                                    tACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGt                                                                                                                                                 <  1:231167/171‑1 (MQ=38)
                                                                                                                                      acGTTGAAGCCGGCAACACCCAGGGCAGCGGCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGa                                                                                                                                                                    <  2:24344/149‑1 (MQ=255)
                                                                                                                                       cGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCACCGACTCGATCGACACCACTACTGCGACCCTGACCGCCTCGCCAAACGTGACCGAAgccc                                                                                               >  1:1105825/1‑215 (MQ=255)
                                                                                                                                            aaGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCCCCTACACCGTGAcc                                                                        >  1:441722/1‑236 (MQ=255)
                                                                                                                                            aaGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGAcc                                                                        >  1:1230353/1‑236 (MQ=255)
                                                                                                                                            aaGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGa                                                                                                                                                                    >  2:773102/1‑144 (MQ=14)
                                                                                                                                            aaGCCGGCAAGACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGa                                                                                                                                                                    <  1:773105/144‑1 (MQ=17)
                                                                                                                                              gCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGCGAACGCCACAGGCg                                                                                                                                                                                                              >  1:218421/1‑100 (MQ=18)
                                                                                                                                                 ggCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTg                                                                                                                                   >  2:1252177/1‑172 (MQ=255)
                                                                                                                                                 ggCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGa                                                                                                                                                                    <  2:855851/139‑1 (MQ=12)
                                                                                                                                                    aaaaCCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACAGGTCGACTGTCAGCGTCCCCATCGAGCACGCCACAGGCggc                                                                                                                                                                                                            >  2:217701/4‑95 (MQ=11)
                                                                                                                                                    aaGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCAc                                                                                                                                                                <  1:618107/140‑1 (MQ=18)
                                                                                                                                                    aaGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCAc                                                                                                                                                                >  2:618104/1‑140 (MQ=17)
                                                                                                                                                    aaGACCCAGGGCAGCGTCGATTTCCAATCCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTCACCCCGAACCCAACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAggcggc                                                                                             <  2:768859/207‑1 (MQ=255)
                                                                                                                                                     aGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCg                                                                                                                                                                           <  2:95798/128‑1 (MQ=38)
                                                                                                                                                        cccAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCGTCGAGAACGCCACAGGCGGTAACTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGt                                                                                                                                                 <  2:351846/151‑1 (MQ=16)
                                                                                                                                                         ccAGGTCAGCGCCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAACTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGAcc                                                                                                                              <  2:1103775/169‑1 (MQ=255)
                                                                                                                                                               cAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCa                                                       <  2:1230350/234‑1 (MQ=255)
                                                                                                                                                                     cGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAAttt                                                                                                                                                                                                       >  1:287890/1‑84 (MQ=33)
                                                                                                                                                                     cGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAAATTCGAGCAACCGAccc                                                                                                                                                                                         >  2:230839/1‑98 (MQ=255)
                                                                                                                                                                                     gCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAggcg                                                                                               <  2:861248/172‑1 (MQ=255)
                                                                                                                                                                                      cGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCAGCGAGAACGCCACATGCGGCAATTTc                                                                                                                                                                                                      >  1:348246/1‑68 (MQ=11)
                                                                                                                                                                                       gAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGAACCCGGCCCAGACCACGATCaacc                                                                                                                                                        >  2:149905/1‑110 (MQ=11)
                                                                                                                                                                                        aaTGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCa                             >  1:1058421/1‑235 (MQ=255)
                                                                                                                                                                                                       aacGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGa                                                                                                                                                                                <  2:993147/73‑1 (MQ=39)
                                                                                                                                                                                                            ttCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGATCGATACCACTACTGCGAcca                                                                                                                             >  1:602164/1‑118 (MQ=255)
                                                                                                                                                                                                             tCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGa       <  1:1126978/236‑1 (MQ=255)
                                                                                                                                                                                                             tCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGCCCCCGACCCCGCCGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAAc                           <  1:292656/216‑1 (MQ=255)
                                                                                                                                                                                                              cGACTGTCAGCGTCCCCATCGAGAACGCCACAGGCGGCAAtt                                                                                                                                                                                                        >  2:263042/1‑42 (MQ=14)
                                                                                                                                                                                                              cGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGt                                                                                           <  2:742250/151‑1 (MQ=255)
                                                                                                                                                                                                              cGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGt                                                                                           >  1:742253/1‑151 (MQ=255)
                                                                                                                                                                                                               gACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAATTTCGAACAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAggcggc                                                                                             <  1:1155751/148‑1 (MQ=255)
                                                                                                                                                                                                                 cTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCCCAGACCACGATCACCGACTCGGGCGACACCACTACTGCGACCCTGACCGc                                                                                                                      >  1:146916/1‑121 (MQ=255)
                                                                                                                                                                                                                 cTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAACAACTGACCCCGAACCCGACGCCGGCCCAGACCACGATCACCGACTCGGTCGACACCACTACTGc                                                                                                                                  >  1:1065952/1‑109 (MQ=255)
                                                                                                                                                                                                                  tGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCgg  <  1:1252180/236‑1 (MQ=255)
                                                                                                                                                                                                                  |                                                                                                                                                                                                                                           
GGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGG  >  AE015451/197422‑197867

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